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Summary. — To introduce the concept of impact parameter into the 
analysis of high energy events seems possible and also useful. Showing 
the inadequacy of Fermi’s theory with impact parameter, we propose 
another way to infer the impact parameter for each event in the emul- 
sion from the experimentally observable quantities. The impact para- 
meters introduced in this way have admissible order of magnitude and 
are shown to have a somewhat definite functional relation with the mul- 
tiplicity of mesons whose plot versus primary energy has been known 
to suffer large fluctuation. Consistency with the uncertainty principle 
is also discussed. 


1. — Introduction. 


Since HEISENBERG pointed out the possibility of the explosive production 
of meson in high energy nuclear impacts, (1), many people have attempted to 
reveal, both theoretically and experimentally, characteristic features of this 
phenomenon, the investigation of which is expected to bring us into a glittering 
expanse of virgin country. 

Various opinions have been stated by different authors. HEISENBERG (') 
has emphasized the possible existence of second kind interactions and the 


(1) W. HEISENBERG: Zeits. f. Phys., 101, 35 (1936); 113, 61 (1939); 126, 569 (1949); 
133, 65 (1952). 
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classical behavior of highly excited states; FERMI (?*), the statistical nature 
of many meson systems; LANDAU (*), the inadequacy of particle picture in the 
early stage of Fermis process; LEWIS, OPPENHEIMER and WOUTHUYSEN ‘ey 
the instantaneity of the production process; and so on. 

These theories have indeed constituted convenient references for experi- 
mentalists, whose attempts have thus consisted in deciding which one of those 
features prevails in nature. Unfortunately, the rarity of events and therefore 
the poorness of statistics hindered the decision. As is very natural in the 
traditional application of quantum mechanics, only one parameter, the energy 
of primary nucleons in the center-of-mass system (C.M.S.) 2y,Mc?, has hitherto 
been used to specify the initial condition in all the theories but those of 
FERMI (*). The analysis of experiments has also been carried through it. For 
example, the multiplicity of meson production n; (*) has been analysed through 
the relation »,= f(y.), which is predicted differently theory by theory and 
therefore considered to serve as a crucial point. However, this way of ana- 
lysis has been not so successful because of the poorness of statistics and of the 
large fluctuation of the multiplicity itself. Taking this situation for granted, 
it seems desirable to find another way to push our investigation forward. We 
have to adapt our method of investigation to the new field of physics under 
consideration, the extremely high energy region. 

Recently, angular distributions of secondary particles have been found to 
vary considerably event by event, even when we restrict ourselves to those with 
similar primary energy. In the usual scattering experiments, we at first 
classify events with respect to their primary energy and then add up the 
number of events to give the angular distribution which is to be compared 
with theory. In this paper, however, we give the term « angular distribution » 
a meaning different from that usual in ordinary scattering experiments; our 
angular distribution is the distribution of many secondary particles, one for 
each event. To compare the distributions with the theory hitherto proposed, 
we must average over the events with the same or at least a similar primary 
energy. 

In the multiple meson processes, there are many quantities which are ob- 
servable in the emulsion, say, multiplicity of mesons, angular distribution, inelas- 
ticity ete.; they are fluctuating event by event. Instead of averaging, we may 


(*) E. FERMI: Prog. Theor. Phys., 5, 4 (1950). See also reference (18). 

3) E. FERMI: Phys. Rev., 81, 115 (1951). 

(4) L. D. Lanpav:' Usp. Akad. Nauk SSSR. 17, 51 (1953); S. Z. BELEN/KIJ and 
L. D. LANDAU: Usp. Fiz. Nauk, 56, 309 (1955). 
(5) H. W. Lewis, J. R. OPPENHEIMER and S. A. WOUTHU YSEN : Phys. Rev., 78, 
127 (1948). 

(*) In this paper the number of charged swift particles is used as a representative 
of the multiplicity. 
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also inquire about the correlation between any two of these fluctuating quan- 
tities. When there are many quantities correlated with each other, we may 
combine them with some parameters (parametric representation). In addition 
these parameters themselves will have some physical meaning. 

Thus, I should like to point out, in this paper, that the experimental facts 
suggest the possibility to introduce one more parameter other than the primary 
energy, the impact parameter b of nuclear collision, and then to analyze the 
multiplicity, for example, through n,= g(y., b). This point of view will enable 
us to draw more information from experiments. 

There may indeed be many questions as to the introduction of such a para- 
meter, to the discussion of which we will return later. One remark to be 
stressed here is that the usefulness of the impact parameter, a rather geometrical 
concept, is a characteristic feature of high energy phenomena. 

To conclude this section, we summarize the experimental material to be 
used in this paper. Besides the number of charged swift secondaries »,, the 
inelasticity AK and the Lorentz factor of colliding nucleons in the C.M.S. y., 
we use the index of peakedness of the angular distribution m. This index m 
is, according to T. F. HoANG (5) and the Bristol group (7), defined through 
the angular distribution (C.M.S.) of secondary particles as 


d n;(0) 


CO LUE 
Ol COS) 


(1.1) 


Although some authors (5) call our attention to the fact that the angular 
distribution is not monotonic even in one of the separate cones, forward and 
backward, but has a maximum in each cone, we shall not, at present, enter 
into such a discussion. The Bristol group (7) showed that m can be conveniently 
calculated by the formula: 


(1.2) {cosec 0»,,= 1.25Vm (for large m), 


the average being taken over the charged swift secondaries. Hereafter we 
shall regard Eq. (1.2) as the equation defining m (*), keeping in mind its 
meaning expressed in Eq. (1.1). 


(6) T. F. Hoane: Journ. Phys. et Rad., 14, 395 (1953). 

(7) B. Epwarps, J. Losty, D. H. Perkins and J. ReyNoLDS: Phil. Mag., 3, 
237 (1958). The data of this paper are supplemented by a private communication from 
Dr. J. IWADARE, to whom the present author is very grateful. The m for the event 
p20 is taken to be 2000. 

(8) K. Niu: Lecture at the Annual Meeting of the Physical Society of Japan, (Oct. 
1957). Uchusen-Kenkyu (Mimeographed circular in Japanese), 3, 85 (1958). The paper 
will appear in Nuovo Cimento; G. Cocconi: Phys. Rev., 111, 1699 (1958). 

(*) Use is made of this definition also for small m. 
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2. — The Fermi theory with impact parameter. 


FERMI developed a theory to describe the angular distribution of each 
event. It is the only theory hitherto proposed which takes the impact para- 
meter into account. We shall briefly summarize his results and compare them 
with experiments, showing some discrepancies. In the following, the equation 
cited as, for example, (F.20) refers to the Eq. (20) in FERMI'S paper (*). 

Fermi’s theory gives for the angular distribution of secondary particles 
when observed in C.M.S., the formula (*) 

dn, 


(F.21) dE 9(B) falon) ; 


where g(B) is a certain function of f, f corresponding to the hot spot temperature 
p= c/JkT, and 


(F.22) falc) = [2/o2(1 — «?)] — (1/02) In [1+)/1— 0]. 


n = cos 0, and o is related to the impact parameter b by (F.20), b/R—3j,(0)/2f,(0) 
(R= luc, u—=T-meson mass), as to the right hand side of which the reader 
may consult the Table I of his paper. 

The numbers of charged secondaries are given by (F.15) and (F.17) after 
some modifications of numerical coefficients to take the charge independence 
into account (a > 3a, b + 3b, and taking two thirds of n for n,): 


(2.1) ns =3.4[2K(y.—1)]*-f(o)/[falo)]* , 
where 

fe) = A+ e*)/o°] In (1 + 0)/(1 — 6) — 2/0? , 
and 


fle) = [1/e] In + 0)/1— 0) + 2/(1 — e). 


The angular distribution (F.21) agrees qualitatively with experiments. As 
a first step towards the quantitative comparison we shall take up the value m 
defined in Eq. (1.2), the agreement of which with experiment is a necessary 


() Some changes of notation: y, r in Fermi’s paper are changed into B and d 
respectively. 
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condition for the theory to be correct. According to Eq. (1.2) and (F.20), 


(2 
i f.(0 cos 0) dé 
(2.2) 1.25Vm = — 


1/2 
PAC cos 0) sin 0 do 


uv 


That m is independent of y, is a characteristic of Fermi’s theory. 
Thus comparing Eq. (2.1) and (2.2), we see that, in Fermi’s theory 
n,/[2K(y,— 1)} is connected uniquely with m through 0, and therefore through 


SG | 
ia 
VD | 
| 
“i 14 20n'-8 
P Sg 23 =| Hel 
ape 1 


ep25-7 


Fermi theory 


Fig. 1. — The relation between multiplicity of charged particles n, and m, the latter 
being defined by Eq. (1.2). The curve shows the relation due to Fermi’s theory. The 
points with question mark correspond to the events with apparent lack of forward- 
backward symmetry in C.M.S. The numbers beside the points show the number of 
gray tracks. For the events without gray tracks and for the events with m < 2 (gray 
tracks < 4), nothing is written. The events with m< 2 are separately listed. The refe- 
rences are designated: p25 etc.: Bristol group (7), D: DEBENEDETTI et al. (°), Dm: DE- 
NTU en) eH LOAN GAS) SH M EOPPER CNT (>). PICKUP en al. (LE), 
S: ScHEIN (13), T: TEUCHER et al. (14). The inelasticities for the events of Pickup et al. 
are taken to be 1, those of the Bristol group to be their K’. 


(°) A. DEBENEDETTI, C. M. GARELLI, L. TALLONE and M. ViconE: Nuovo Cimento, 
4, 1142 (1956). 

(10) M. DEMEUR, C. DitwortH and M. ScHONBERG: Nuovo Cimento, 9, 92 (1952). 
(1) V. D. HoPPER, S. Biswas and J. F. Darsy: Phys. Rev., 84, 457 (1951). 
(2) E. Pickup and L. Voyvopic: Phys. Rev., 84, 1190 (1951). 
(13) R. G. GLASER, D. M. HAsKin, M. ScHEIN and J. J. Lorp: Phys. Rev., 99, 
1555 (1955). 

(14) K. GorTsTEIN and M. TEUCHER: Zeits. f. Naturf., 8a, 120 (1953). 
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the impact parameter b. While a certain relation between them is apparent in 
the experimental plot (cf. Fig. 1) (*), the discrepancy with Fermi’s theory is 
so deep that the conventional procedure to adjust the hot spot volume 
4xR*/3 (to translate the theoretical curve upwards or downwards) does not work. 

That the experiments show values for n,|[2K(y.-— 1)? higher than those 
expected in the theory, seems to show the inadequacy of Eq. (F.21) and con- 
sequently of Eq. (2.2). It seems probable that, in contrast to Fermis theory, 
m depends also on y, in such a way that the higher value of m does not neces- 
sarily correspond to the more distant collision. 

Thus, it is shown that Fermi’s theory is not adequate for us to infer the 
impact parameter of each event from the angular distribution. In the next 
section, we shall provide another method for this purpose. 


3. A method to infer the impact parameter. 


The discussion in this Section is very crude, and the method presented is 
rather heuristic in nature. 


(i) To begin with, we shall consider the physical meaning of m intro- 
duced in Section 1. Since m is originally defined to be the power of cos 0 
in the angular distribution of mesons, it may have some relationship to the 
angular momentum of each meson; the following reasoning will confirm this 
to some extent and lead us to the conclusion that «the angular momentum 
of each meson is probably perpendicular to the direction of incidence of a 
nucleon having initiated the process (this direction is hereafter referred to as 
z-axis), and the most frequent value (the mode, in the terminology of sta- 
tistics) of the angular momentum among secondary mesons of the jet in 
question will be mA/2 ». 

As described in the introduction, our angular distributions are different 
from those in the usual scattering experiments which are the average taken 
over impact parameters. When we restrict ourselves to the event with definite 
impact parameter, we can, in principle, calculate the wave function of secon- 
dary particles using the wave packet for the primary nucleons. 


a) Consider a case where the angular part of the wave function for 
the secondary meson-system is given by 


(3.1) ‘(03 Pi ; 0.923 ‘003 0,n) = DI P ‘Pi (cos 0,) exp [— imp; | 5 
P 


‘Pia (cos 0;) exp [— D; Sch “Pi (cos 0,.) exp [— iun®;,] s 


(*) As for the Bristol data (7), their K’ is taken for our K. 
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where > P is the symmetrizing operator, P*(cos 0), the associated Legendre func- 
= 
tion. 4’s and w’s are assumed to be different from each other. Then the 


angular distribution of mesons is: 


dn 


9 ne 
(3.2) 10 


x [ae ose do, AE 
oc >| Phi(cos 0)|?, 


since the interference term vanishes on account of our assumption. 
Alternatively, consider the other extreme case 


MERA) A= p= =, =) Im (3.1), then 
dn , Ane 
(3.3) 10% [Pz (coso) 


In fact, the wave function of secondary mesons may be more complicated 
than those assumed above; the angular momentum state for each meson may 
not be pure, and so on. We cannot at present enter into such a discussion. 
In connection with this, the observation of possible azimuthal anisotropy may 
be of interest (1°). 


b) Experiments show us that the angular distribution is peaked for- 
wards and backwards, and therefore that we may assume u;,= 0 (cf. Eq. (3.2), 
(3.3)). The angular momenta of mesons are thus perpendicular to the 2-axis. 


c) Owing to the fact that, in experiments, the angular distribution 
is expressed by the term cos” 0 alone, m/2 may be interpreted as the 4 in 
Eq. (3.3) or the À around which the density of À, in Eq. (3.2) is.greatest. This 
is the reason why we have called m7/2 the most probable angular momentum 
among secondary mesons of the jet in question. 


d) One additional remark. Although P#(cos 0) contains also terms 
cos 0, k<A, it is the highest power term cos” 0 that is responsible for the peaks, 
forward and backward. And the shape of the angular distribution near 0 = 0 
or 7 is important in the determination of the index m, because of the weight 
factor cosec 0 in the defining equation (1.2). Thus the neglect of cos" 6, 
k< À, and the choice A= m/2 is justified (*). 


(5) K. Niu: Soryusiron-Kenkyu (Mimeographied circular in Japanese), 7, 836 (1955); 
Z. KoBa and S. TAKAGI: Nuovo Cimento, 10, 755 (1958). 

(*) Mr. K. KoBAYAKAWA, Osaka University, remarked that we have to pick up 
the term cos* 0 with the greatest coefficient in the expansion |P,(cos 0)|?, then to 
take m= k. The change in guess into A= «am (x— some numerical constant other 


than 4) is not so serious to our present discussion. 
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ii) Now the knowledge of the total angular momentum of secondary 
mesons J, allows a guess on the impact parameter 6 through the conser- 


vation law: 


(3.4) MODE (M — nucleon mass). 


The left-hand side is the angular momentum of two nucleons in the initial 
stage, that is still perpendicular to the ¢-axis. In the right-hand side, the omis- 
sion of the angular momenta carried off by particles other than mesons may 
be admitted (*), taking into account the crudeness of our argument. 


iii) Let us consider how to obtain the total angular momentum J,,. 
No detailed mechanism of the production process being as yet known, we 
cannot but rely upon a probability-theoretical argument (*). According to the 
theory of random walk (15), the probability distribution for the resultant R 
of N (>1) vectors x; whose magnitudes and directions are randomly distri- 
buted over a plane (2-dimensional random walk), 


R? 


1 
3.5 ST dog, 
(3.5) W(R)dR = —— esp | aa ote 


That this distribution is characterized only by o?= <a>,, is the merit of our 
assumption N>1. 
The most probable value of R is obtained from 


d 
— [2aR-W(R)]=0, 


dk 
as 
(3.6) Re VINo? | 
Now identifying 
(3.7) sp ae ra Se I 


we obtain for the most probable resultant angular momentum, 


1 3n, (mh\2 = 
(3.8) Jr VE i (0) E ve nee 


2 


() Cf. K. Kopayakawa’s comment, which will be published in the near future. 
Vo K. Husimi: Kakuritsu-ron oyobi Tôkei-ron (Kawade-Syobò, 1940), p. 215. 
It is not difficult to take into account the Markoffian correlation, into the discussion of 


which we have not entered. See loc. cit. p. 260; S. CHANDRASEKHAR: Rev. Mod. Phys. 
15, 1 (1943). i 
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iv) Thus eq. (3.4) gives, as the most probable value (*) 


V3 mVn, h 
. mV nN, a (M: nucleon mass). 
4 Ve Me 


(3.9) b 
The impact parameters obtained in this way are found to range from (1/50) (h/Me) 
to 10(%/Mc), admissible order of magnitude, showing that Eq. (3.9) is suffi- 
ciently trustworthy. 

Other facts supporting Eq. (3.9) are the systematic functional relations 
found between b and some other quantities, for example, multiplicity of me- 
sons (see next section). 

As to the statistical distribution of the impact parameter over all events 
observed in the emulsion, b does not seem to obey Prob (r<b<yr-+ dr) a 
oc 2ardr. It is not yet clear whether this results from the poorness of statistics 
or more profound nature of the process in question. 

When we assume a law n,o y?b * partially invoking the statistical theory 
of multiple production, then 


(3.10) m oc bite - 


This shows, in contrast to Fermis theory, that m is not only a function of b, 
but also an increasing function of y,: this agrees with the guess at the end of 
Section 2. 


4. — The impact parameter and other observable quantities. 


In this Section some graphs are presented to show the possible existence 
of systematie functional relations between the quantities observable in emul- 
sion experiments and the impact parameter introduced in the last section. 

Besides the impact parameter, there is another parameter that character - 
izes the initial condition, the energy of nucleons y,; it is therefore desirable to 
investigate the relation, for example, in the form ns= g(y.; d). In the fol- 
lowing, we draw graphs with impact parameter b as abscissa and discriminate 
the energy y. by some marks. The dependence on y, will be found weaker 
than that on b. 

Since the argument in the last Section was probability-theoretic in nature, 
so should be the relation of any quantity with b. The fluctuations seen in the 
graphical representation n,=g(y.,b) are definitely smaller than in the usual 


(*) As to the events with m= 0 see Sect. 5. 
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representation n,= f(y.). This fact shows the importance of the impact para- 
meter for the analysis of experiments. 
The experiments referred to are the same as in Fig. 1. 


i) Multiplicity and impact parameter. See Fig. 2. The question marks 
beside the points tell the lack of forward-backward symmetry in the corres- 
ponding events. The numbers just after the hyphen mean the number of 


= 
= 
= 
i (ui | 
0 25/20 ApII NS 
P p5 i Tea (pis excluded ) 
© 
4 A p20n1-8 APS 
py x - 
È A x p25-7 CR eP10 
17 ® pl6- 
è i 8 |, SE) pl aP24 
ox e 
e” @x 0) pîa 
b x 5 e: 
p15 
<a è? p20 
DMHHEPT pig 
m<2 | m=2 | 
i — —_T_ 
005 01 


1 5 10 
b (unit:/Mc : M=nucleon mass; 


Fig. 2. — The relation between the multiplicity of charged particles n, and the impact 
parameter calculated through Eq. (3.9). The experimental materials are the same as 
those of Fig. 1. 


heavy tracks accompanying the events. For example p25-7 is accompanied 
with 7 heavy tracks and is inadequate for our analysis. The events with 
m< 2 are separately listed (*). 

The multiplicity (of charged particles) n, decreases rather systematically 
with increasing impact parameter d; this is very natural. The deviation 
of those points with question mark from such a trend may be considered to 
support, though negatively, our method of analysis. From the graph, roughly, 


TOS Oa aa 


However, this relation should not be taken seriously because of many uncertain 


(*) Under the restriction to integer value, m (<2) must be zero on account of the 
forward-backward symmetry of the angular distribution. 
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factors in the argument of Section 8. It is the only aim of this paper to 
indicate the possible rôle taken by the impact parameter. 


il) Inelasticity and impact parameter. — The Bristol group (7) introduced 
two kinds of inelasticity, A’ and AK”. K' is defined on the assumption that all 
the secondary particles are 7-mesons with transverse momenta 0.5 GeV/c. K” 
is calculated on the same basis but excluding the particles (about 20%) emitted 
at the smallest forward and backward angles as «x » particles different from T's. 


Afa 100 £Ye 
© 50 <Y <100 
e 30 <7ç < 50 
® 10 <}ÿe < 30 
005 01 SO À HS À 
K' 
1 
05 
04 
03 
02 E 
Fig. 4 
fC. pede RE ee 
si 005 01 GS 5 10 


b (unit: Mac » M=nucleon mass) 


Fig. 3 and 4. — The relation between the inelasticity and the impact parameter calcu- 
lated through Eq. (3.9). The experimental materials are the same as those of Fig. 1. 
As to the definitions of A’ and K" see the text. 


We identify the inelasticities used by authors other than the Bristol group 
with K'. (See Fig. 3 and 4). We see, in the region of small impact parameter, 
the fact that the inelasticity decreases with the decrease of the impact pa- 
rameter; this trend is surprising. However, the fact that this trend is more 
prominent for K” makes us surmise the increase in the proportion of «x» 
particles (*) as the impact parameter decreases. If this is in fact the case, 
the investigation on high energy collisions will yield some knowledge about 
nucleon structure. 


(*) This terminology was introduced by the Bristol group (°). 
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5. — The impact parameter and the uncertainty principle. 


To introduce the concept of impact parameter, we have to examine its 
compatibility with the uncertainty principle. We shall, meanwhile, consider 
the problem in the laboratory system (hereafter abbreviated as L.S.), the z-axis 
being taken parallel to the direction of incidence of the primary nucleon. 
Denoting the uncertainty in the position and momentum of the incident nuc- 
leon Ag, ...,; Ap,, ..., respectively, the uncertainty relations are 


h 
(5.1) AGIND = =e CU: 


D 

To make the impact parameter b meaningful, the condition Ax < b must be 
satisfied. When the momentum p of the incident nucleon is large enough, 
Ap may also be large and Aw small, allowing us to use the concept of impact 
parameter; this is the conventional argument. In addition to the largeness 
of Ap, however, there must be an appropriate relation among the Fourier 
amplitudes to realize a small wave packet. (For example, in the non-rela- 
tivistie theory, to obtain a minimum packet for which Ax-Ap, = h/2, 


p(x) SI exp [ip,~|dp, ’ 


we must take the amplitude f(p,) to be a function of gaussian type.) As to the 
realization of the latter relation we know nothing. An observation of position 
causes the wave packet to contract. But observation with an atom (ionization) 
is insufficient to localize the nucleon in such a way as to be able to speak of its 
impact parameter. The idea which arises next is to observe, in the emulsion, 
many ionized atoms. Does a line suitably drawn to pass them enable us to 
find the nucleon path more accurately than observing one atom? This does 
not seem to be the case even when we note that the nucleon in question is 
extremely relativistic and that the spreading of the wave packet is very slow 
because of the well known retardation of a moving clock, since the momentum 
transfer to an atom is very small and randomly oriented for different atoms 
in emulsion. 

Then, how can we understand the regularity found in Section 4? The 
present author can’t answer this question. In this connection, we remember 
that the direction of motion can be determined to 10-4 radian in the emulsion 
experiments; this fact implies that Ap,/p, of a particle traversing emulsion is 
smaller than 10-4, and that 

TLC EN? 


h 
(5.2) NO ss}. eee 
2AP: 2p 4y? Me” 
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Thus, for y,— 100, the impact parameter (*) may be definite to the extent 


h 


5.3 ~ : 
ae 4 Me 


This fact will afford a basis for the introduction of the impact parameter. 

In connection with this, we notice the following facts: there is only one 
event found with m+ 0 for y.< 30, and few events with m=0 for y. > 30 (**). 
See Fig. 2. The reader may be perplexed with the fact that the events with 
m= 0 (which would have b~ 0 if Eq. (3.9) were applied) have rather small 
mutiplicity. The reason is that the events with m = 0 are also of low energy. 
These facts seem to show the inadequacy of the impact-parameter concept 
for low energy events. 

The above arguments show only the possibility of introducing the impact 
parameter and it should be decided only by further experiments whether or 
not the nucleon-wave packet is so small as to make the parameter useful. 


6. — Conclusions and discussions. 


i) The importance of the impact parameter for the description of expe- 
rimental results in the high energy region is pointed out. In this connection, 
Fermi’s theory with impact parameter is examined. The relation between n, 
and m, resulting from this theory, is shown not to be in accord with experiments ; 
and so this theory cannot be used to infer the impact parameters for each 
event. An alternative method, heuristic in nature, for this purpose is proposed. 
In spite of the crudeness of reasoning, the impact parameter, calculated in 
the latter way from experimentally observable quantities, is satisfactory in 
various respects: (i) admissible order of magnitude 


edo (M: nucleon mass), 


and, (ii) rather definite functional relations with multiplicity and probably 
with inelasticity. 


(*) For the sake of eliminating the uncertain situation as to the target nucleon, 
the reader may consider the problem in C.M.S., remembering the Lorentz trans- 
formation of angle cms. = 2:y.Ous.- 

(**) See footnote on page 754. 
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The concept of impact parameter is possibly compatible with the uncer- 
tainty principle for the high energy region under consideration. However no 
positive basis that the nucleon wave packet must be small enough could be 
found. 


ii) Now that the importance and utility of the impact parameter for 
high energy events is established, there arise many possibilities for theory and 
experiments, the examples of which will be described below. 

For experiments, it is interesting not only to confirm the relations men- 
tioned in Section 4, but also to investigate the composition of secondary par- 
ticles as a function of the impact parameter. Especially it is interesting to inves- 
tigate whether there arises a knick-like behavior in the functional relation 
with b. Isn’t there, for example, a sudden increase in the nucleon pair pro- 
duction at some small values of 6? This knowledge should contribute much 
to our understanding of the nucleon-structure. 

Theoretically, it becomes possible to treat, as a first step, only events with 
large impact parameter; a rather easier problem being separated from the com- 
plicated whole. This method is the same as that which is advocated by Prof. M. 
TAKETANI (17) and is adopted by Japanese workers on nuclear force problems. 
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like to note here, with grateful recognition, that Prof. Y. FUJIMOTO and J. 
NISHIMURA have called his attention to the angular distribution problem of 
high energy jets, and that K. KoBAYAKAWA gave him critical comments. Last 
but not least, the author expresses his gratitude to the Iwanami Fûjukaj 
Fellowship for the financial aid. 


(1) M. TAKETANI, S. NAKAMURA and M. SASAKI: Prog. Theor. Phys., 6, 581 (1951). 
See also, Suppl. of Prog. Theor. Phys., 3 (1956). 


Note added in proof. 


(i) Although we have not been able to find any positive basis that a small 
nucleon wave packet is realized in nature, we may still make some arguments in 
favour of the concept of impact parameter, from another point of wiew. 

When the incident nucleon is duly expressed by plane waves rather than by a 
wave packet, we may ascribe impact parameter b = 1%/p to the partial waves with 
angular momenta ~ l# (p being the momentum of the nucleon). Doubtful as it may 
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appear at first sight whether we can ascribe an impact parameter to each event in 
the emulsion (many partial waves may contribute to the cross section), this doubt 
will be resolved when we can assume the phase relations to bo very complicated 
among the matrix elements of those partial waves. Since interferences are destructed, 
the cross section will be practically identical to that which would be calculated for 
the « Gemisch » of nucleon states, each of these states belonging to a rather small range 
of angular momentum. If this is the case, a member of the « Gemisch » corresponds 
to each event in the emulsion, thus allowing the concept of impact parameter. This 
is again a characteristic feature of high energy phenomena. 

Other points of view will also be possible. 

These are essentially the same as the arguments given by Prof. Z. KoBA on the 
Annual Meeting of the Physical Society of Japan (Oct. 1958). Also, the author owes 
much to the members of the symposium on Interpretation of Quantum Mechanics 
and Quantum Field Theories (Feb. 1959). To all of them, the author expresses his 
sincere gratitude. 


(ii) We have to note, with great pleasure, that Koba and Takagi’s introduction 
of «classifying variables » is similar in spirit to our introduction of the impact parameter. 
Their argument is adapted for the framework of field theory. See, Z. KoBA and 
S. TAKAGI: Fortschritte der Physik, in press. 


RIASSUNTO (*) 


L’introduzione del concetto di parametro d’impatto nell’analisi degli eventi di alta 
energia appare possibile ed opportuna. Dimostrando l’insufficienza della teoria di Fermi 
col parametro d’impatto, proponiamo un’altra via per ricavare il parametro d’impatto 
per ogni evento in emulsione in base alle grandezze sperimentalmente osservabili. 
I parametri d'impatto introdotti in tal modo hanno ordine di grandezza ammissibile 
e, come si dimostra, una abbastanza definita relazione funzionale con la molteplicità 
dei mesoni il cui diagramma in funzione dell’energia primaria presenta, come è noto, 
grande fluttuazione. Si discute anche la compatibilità col principio d’indeterminazione. 


(*) Traduzione a cura della Redazione. 
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Extended Invariance Properties of Quantum Fields. 


I. — Lagrangian formalism. 


M. E. MAYER (*) 


Joint Institute for Nuclear Research - Dubna 
Laboratory of Theoretical Physics 


(ricevuto il 15 Novembre 1958) 


Summary. — This is the first paper from a series in which the consequences 
of invariance of Quantum Field Theory under an f-parameter Lie group 
are investigated. Part I deals with the consequences of invariance of a 
general Lagrangian under such a group. It is shown that the invariance 
leads to the conservation of f «current density operators», the corres 
ponding charge operators being a representation of the Lie group on tl : 
state vector space. Invariance under extended transformations, i.e. tra? s- 
formations in which the parameters are functions of the space-time po nt 
leads to the necessity of introducing f non-commuting vectors fields. The 
general structure of the Lagrangian corresponding to these fields is est ab- 
lished and their possible physical counterparts are discussed. Part II 
will contain a similar investigation in the framework of an S-matrix 
theory, and Part III will be devoted to a detailed classification of pos- 
sible transformation groups and the corresponding calssification of par- 
ticle fields and their interactions. 


(‘) On leave of absence from the Faculty for Mathematics and Physics of the 
C.I. Parhon University, Bucharest and the Institute for Atomic Physics of the Aca- 
demy of the Rumanian People’s Republic, Bucharest. This work was started at these 
Institutions in 1956. Present address: Physico-Mathematical Faculty, C.I. Parhon 
University, Bucharest. 
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1. — Introduction. 


During the last years considerable progress has been made in the study 
of elementary particle interactions and the selection rules governing them, 
and there have been several, more or léss successful attempts at theoretical 
interpretations or classifications of the corresponding interactions, the best- 
known of which is the GELL-MANN-NISHIJIMA interpretation in terms of the 
strangeness quantum number, and the various modifications of this scheme. 

Although it has been clear that any conservation law in quantum field 
theory must be the consequence of some invariance property of the fields under 
investigation, there have been relatively few attempts to study in quite a 
general form the consequences of the invariance of a theory (we will specify 
below what exactly is meant by this term) and to investigate the possibility 
of obtaining all possible fields and their interactions on the basis of the funda- 
mental postulates of quantum field theory and the structure of the admissible 
transformation groups, which, as will be shown below, are Lie groups which 
turn out to possess the important (from the mathematical point of view) pro- 
perty of being compact and which are well studied in the mathematical lite- 
rature. 

Such an investigation has been started by the author in 1956 and some of 
the results obtained have been quoted in an unpublished dissertation (1). 
Essentially the same results have been published at about the same time by 
UTIYAMA () and are mostly a generalization of Schwinger’s analysis of 
constant phase transformations and their extension to gauge transformations, 
which leads to the introduction of the electromagnetic field (**). These results, 
together with some general properties of the groups which leave invariant the 
Lagrangian form of quantum field theory, form the contents of Sections 2 
and 3 of this paper. It is shown that the « internal degrees of freedom » (as 
contrasted to the space-time properties) of the fields are such that the «in- 
ternal space » is a finite dimensional Euclidian space so that the corresponding 
transformations are orthogonal lnear transformations. The invariance of the 
(quantum) Lagrangian with respect to such transformations leads to the con- 


(1) M. E. MAYER: Some Functional Methods in Quantum Field Theory (in rumanian). 
Unpublished Dissertation, Parhon University, Bucharest (1956). 


(2) R. UTIYAMA: Phys. Rev., 101, 1597 (1956). 

(3) J. SCHWINGER: Phys. Rev., 91, 713 (1953); Proc. Nat. Acad. Sci. U.S.A., 44, 
223 (1958). 

(4) M. E. MAYER: Cimpuri Cuantice si Particule Elementare (in rumanian), (Bucha- 
rest, 1959). 
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servation of f charges, where f is the number of parameters char: cterizing the 
given group, and the commutation relations of the corresponding operators 
coincide with those of the infinitesimal operators of the group, so that the 
charge operators are a representation of these infinitesimal operators in the 
state-vector space of the system. | 

In Section 3 the notion of «extended invariance » is introduced. This 
means that invariance of the Lagrangian is required not only for transformations 
in which the parameters are constant in space-time, but also under point- 
dependent transformations. Physically such a requirement can be motivated 
by imposing the condition that if a theory is invariant under a group of internal 
symmetry transformations, this invariance must not depend upon the point 
of location of the observer, so that two observers in distinct space-time points 
should be able to carry out independently arbitrary transformations from the 
group under consideration, without destroying the invariance of the observable 
quantities. In order to ensure the invariance of the Lagrangian formalism under 
such extended transformations it becomes necessary to introduce supplemen- 
tary operator fields (in the same manner as the Maxwell field can be intro- 
duced to reestablish the invariance of a charged field Lagrangian, when ex- 
tending constant phase transformations to gauge transformations). The co- 
variant vector fields introduced in this manner are subjected under our group 
to generalized gauge transformations of the second kind (3.7). The possible 
general form of the «free Lagrangian» corresponding to these vector fields 
is also investigated and it turns out, that it leads to complicated non-linear 
field equations, which have not been investigated in the general case. 

A similar analysis, but in the framework of an S-matrix theory is the 
subject of Part II of this series and Part III will contain a detailed analysis 
of the structure of admissible transformation groups with the final scope of 
obtaining a classification of «elementary » particles and the structure of the 
corresponding Lagrangian (this was in part carried out by SCHWINGER (3)) 
or S-matrices. This investigation is only a preliminary step to the construction 
of a quantum field theory in which physically meaningless quantities, such as 
interacting field operators etc., no longer occur and the basic quantities of 
the theory are closely connected with group representations and invariant 
transition amplitudes (in this connection attention must be called to attempts 
in this direction by SEGAL (5), which are however in a much too abstract form 
as yet). 


C) ly He SEGAL TE ro Nat. Acad. Sct. US wa sda 103 (1955); 42, 670 (1956). 
and other mathematical papers to which references can be found there. 
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2. — Lie Groups, infinitesimal transformations, structural costants and current 
densities. 


We take as our starting point Schwinger’s general form of the Lagrangian 
density (3) 


(2.1) LY) == a1 4 Ady — Cy Aly) — Ay) , 


where A” are a set of antihermitian matrices and y sets of hermitian field ope- 
rators. The matrix indices of 4" refer only to the space-time transformation 
properties of the operators yx. 

In order to be able to describe various physical properties of the particles 
associated with the fields 7 the latter must possess some internal degrees of 
freedom (i.e. they must consist of sets of components suitably labelled) which 
can be subjected to various transformation groups. The following natural 
requirements are to be imposed on such «internal transformations »: i) L must 
be invariant under the transformations of the group G, or at least invariant 
up to the addition of a divergence; ii) all physically significant relations, such 
as commutation relations, expectation values of S-matrix elements etc., should 
also be left invariant. From these requirements it follows immediately that 
the internal space, 1.6. the system of components of y, must have some sort 
of scalar product defined in it, which shall assure the positive definite character 
of e.g. the anticommutator of fermion operators or the energy density of bosons 
(cf. also (5)). 

Thus the internal space must be a Riemannian manifold, but if we do not 
want to get entangled with non-linear field equations it is reasonable to sup- 
pose it a finite-dimensional Euclidian space. The corresponding transformations 
of the internal space are consequently linear orthogonal transformations = ro- 
tations. 

Consider the f-parameter group of linear transformations of the n-dimen- 
sional internal space 


(2.2) Ge a DI M,,(41, SOLE, As) Xs ’ 


sel 


The infinitesimal operators of this group can be written in the form 


f 
(2.3) MOOD a OA dd ATE 


a=1 


(6) J. SCHWINGER: Ann. of Phys., 2, 407 (1957). Cf. also Proceedings of the 7th An- 
nual Rochester Conference (1957). 


en 
[el 
wu 
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or in obvious matrix notation 
(2.4) ty = (I — i0/,1")4 


In order to preserve the Hermitian character of the field operators it is neces- 
sary that the matrices 7" be imaginary. 

The A” commute with the 7 as they refer to different degrees of freedom, 
so that the invariance requirement on the Lagrangian leads to the Hermitian 
property of the infinitesimal matrices 


i 


(225) TS =.= Ta 5 (io) — Ta af 


Let the structural constants of the group, Ci’ be defined by the ¢ relations 


f 
(2.6) [LP Te Te Tee SOLE 


K=1 
with the well known properties 
Yip __ ji 
Ci, SALA Ci, ’ 


x (CCE yin na OF Ox ale Oe (09) = 0 È 


m m 
m 


A variation of the Lagrangian with respect to À, can be put into the form (*) 
vi 
(2.8) Ô,L = — DAI UT" y) 8,02, (&), 


and defines the current densities 


which satisfy the conservation laws 


fa ag" (@) 0, 


(2.16) > 
Où = à o"47 (a) = const. 


Q 


() Here 04,(x) are arbitrary functions of the co-ordinates, vanishing on the two 
limiting spacelike surfaces which intervene in the action principle; cf. for example (4), § 21. 
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Repeating the usual argument (e.g. (4), $ 22) it can be seen that 


(2.11) [x O°] = Ty, 
and (2.6) implies 
(2.12) (9, Q*] = > 020", 


(4 
that is, the various charges do not commute in general. 
On the other hand, the transformation (2.2) can be expressed as a unitary 
transformation of the field operators 


vj 


(2.13) 'y=My=U+*yU, U0=UU+=1I, 


with U acting on the state vectors (e.g. in the occupation number represen- 
tation U can be expressed in terms of the creation and annihilation operators). 
For infinitesimal transformations (2.3) U can be expressed by means of Her- 
mitian operators 4 in the form 


(2.14) U=I+iY Q°64,. 


It is easy to verify that the operators Y, satisfy the commutation relations 
(2.11) and (2.12) and are thus identical with the operators (2.10). Thus for 
finite transformations we obtain 


(2.15) Diexpii SQ, "exp 


x | do. : 


a 

Now let us consider briefly the structure of the possible transformation 
groups. As shown in Chapter XI of (7) the Lie algebra of an euclidian space 
is compact and decomposes into the direct sum of its center (the commutative 
subgroup) and a set of simple compact non-commutative algebras which are 
ideals in the original algebra. The vector space of the Lie algebra is unitary 
with the Hermitian scalar product 


(2.16) (a, 6) =— > 00,0 = 9 ,a°b". 


Thus it is sufficient to consider separately commutative and purely non-com- 
mutative algebras. 


(7) L. PONTRIAGIN: Topological Groups (Princeton, 1939); second russian edition 
(1953). 


12 
È 
= 
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The commutative groups, which can be considered as rotations around an 
axis in the internal space, or alternatively as constant phase transformations 
have been thoroughly studied in (*) and (4) and will not be considered here. 
We only note that the « charges » associated with various phase or gauge trans- 
formations commute and can be used for classification of the states of a system 
(this problem will be treated in detail in Part III of this series of papers). 
Examples of non-commutative groups have been investigated by different 
authors (8-11), and some special examples will be investigated below (Part III) (*). 


3. — Extended group invariance and associated vector fields. 


As a generalization of the extension of the group of constant phase trans- 
formations that leads to the introduction of the electromagnetic field with its 
own gauge transformation of the second kind (cf. (4), Sect. 22, or SCHWINGER’S 
original work (*)) we will now consider «extended » transformation groups in 
which the parameters 7, are arbitrary functions of the space-time co-ordinates. 
This extension is based upon the independence of the internal degrees of free- 
dom on the space-time transformation properties and means that we may 
apply to each space-time point its own internal transformation (**). 

Consider the extended transformations 


(onl) ASEM AIA) 
or their infinitesimal form 


(3.2) OE DDA GIO), 


where we have made use of the antisymmetry (2.5) of the matrices 7. The 
variation of the Lagrangian has the form 


(3.3) è,L=—iYyA4T% 


a 


4 


d4,(0)= > fx) 02 (a), 


(8) C. N. Yane and R. L. Mirus: Phys. Rev., 96, 191 (1954). 
(*) A. Sanam and J. C. PoLKINGHORNE: Nuovo Cimento, 2, 685 (1955). 
10) B. D’ESPAGNAT and J. PRENTKI: Nucl. Phys., 1, 33 (1956). 
(7) B. D’EspaGnAT, J. PRENTKI and A. SALAM: Nuclear Physics, 5, 447 (1958). 
and A Tentative Scheme for Weak Interaction, CERN preprint by the first two authors. 
() Cf. also S. BLuDMAN: Nuovo Cimento, 9, 433 (1958). [Note added in proof]. 
(”) Essentially the same result was obtained by UtIYAMA (2) by means of a less 
elementary proof. For the special case of the 3-dimensional iso-rotation group the 
same result has been obtained by YANG and Mirus (7). 
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which is no longer zero, as À(x) does not vanish on the boundary surfaces. 

However we can reestablish the invariance of the Lagrangian by intro- 
ducing supplementary fields UV”, and a supplementary term in the Lagrangian, 
L, such that 


(3.4) Ô(L + Ly) =0, Lie Dit ? 


It is clearly insufficient to subject the field U” to the gauge transformations 
of the second kind known from electrodynamics 


(3.5) Ue > UH — ©} (x) 


as in this case 


(8.6)e OL (AU) pane OA (+6 S x4,[1°, T']UL2,(0) = 
= — Vj 02, (@) — > F(x) UA (2). 


The first term cancels with (3.3) and the second cancels only if we subject 
the vector fields U“ to the generalized gauge transformations of the second kind 


(3.7) ‘UE (e) = U4 (x) — 0A, (a) + Y 08 UEA, (x) . 


In order to have a complete description of the U”-field we must investigate 
the possible form of the corresponding «free Lagrangian », imposing the na- 
tural conditions that it be of the general form (2.1) or reducible to that form 
by means of introducing suitable supplementary variables (since we are looking 
for first order equations), and that the Lagrangian be invariant under Lorentz 
transformations and extended internal transformations (here invariant means 
as usual, invariance up to possible appearance of expressions in the form of 
a four-dimensional divergence). 

As well known from electrodynamics, or the theory of the Proca field, the 
supplementary variables needed to obtain first order equations for the vector 
field, are antisymmetrical tensors of the second rank, so that besides U” we 
introduce the f antisymmetrical tensor fields Z“” for which we suppose the 
following transformation law under the extended transformation 


vo 
D 
— 


OPA Den (Oa ey Yk uy 
TNT ES OZR A, , 
kt 


(there are no terms involving the derivatives of 1;, because of the antisym- 
metry). This means that the internal space components of the vector U and 
the tensor Z transform cogrediently under our group (leaving aside the deri- 
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vative of À,, which for our present purpose can be considered as à (constant »). 

Using the well-known facts from the theory of Lie groups, that the struc- 
tural constants are «tensors » under the transformations of the group, and 
that, for fixed index a they form a matrix representation of the generators of 
the Lie group (one can easily verify that the identity (2.7) is just the repre- 
sentation of (2.6) in an f-dimensional Euclidian space with the scalar product 
defined by (2.16), we are able to write the equations (3.2), (3.7) and (3.8) in 
the form 


Wa = DATA 


b 


(3.9) (Da = SE OSA) TY PI, 


d 


he Pe > (6 A i > 0? A WARE: 


b 


and the «metric tensor » 


(3.10) CES OO COTS ODO 


can be used for lowering or raising internal indices and for forming invariants. 

Taking into account all the mentioned facts, we see that the only inte- 
resting bilinear Lorentz invariants, linear in the first derivatives of U* and 
Zw are 


= DE (0, 0, — DI Gite ot! 0? — z o EZ 
a ab 
(3.11) T,= 3 21,20, I, = > US Zi, 
RADIO Un) Ze 


a.b.c 


These forms are not invariant under the extended group, because of the deri- 
vatives of 2, appearing in (3.9) (they are however invariant, for constant 2,) 
but they can be combined so as to yield an invariant, the various terms with 
derivatives cancelling one another, or by using antisymmetry and neg- 
lecting terms of the form of a divergence. By these arguments we are led 
almost uniquely to the free Lagrangian for the U-Z- field, of the form 


(8.12) L(U, Z) =— 33 {Z}, (0,0; — 2,09} — 083, ZH} — 


st Sipe VASTO } +4 


a? iw 


to} 


ab [TTL TTY MONT 
Lins Cr (Ut he UE): ’ 
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where summation over repeated indices is understood. The first two terms 
in the Lagrangian can be evidently transformed into one another, up to a 
divergence. The form (3.12) reduces to the one used by SCHWINGER (°) for 
pO: 

The resulting field equations (taking into account the coupling with the 
y-field through L, (3.4)) are non-linear generalization of the Maxwell equations 
of the form first proposed by YANG and MILLS (8). 


È x wv “T7 ‘ss T yab TT Tr T 
NES zo = qu SU“ + 3 CU“ U'=— VUE) 
and 

‘ Mi QUIZ _—_ 1 A [TAC ATEN 
(3.14) A NAT 4 TT 


The non-linear character of these field equations leads to great difficulties 
in analysing their quantum properties, which will not be investigated here. 
We only want to call attention to the fact, that due to the vanishing of the 
rest mass of the U-Z field (a term UU? would lead to non-invariance of 
the U-Lagrangian under the extended group) and the existence of the gauge- 
transformation, not all the components of the U-field will be dynamical va- 
riables (the same is well known for the electromagnetic field, which, as yet, 
is the only U-field known to havel'a physical counterpart). 

As already mentioned in the introduction, the only special cases of U-fields 
treated in the literature so far, excepting of course the electromagnetic field, 
are the Yang-Mills B-fields, associated with extended isotopic invariance (”) 
and the gravitational field or more exactly, the field of the Riemann affine 
connection, which as shown by UTIYAMA, plays the role of the U-field for 
«extended Lorentz transformations » i.e. Lorentz transformations with point- 
dependent coefficients (?) (*). UTIYAMA also called attention to the fact, that in 
order to obtain the invariant interaction term between the field y and the 
new U-field, one has simply to replace the ordinary derivative 0,4 by the 
« covariant derivative » 


(3.15) Vi i= XIX Uh 


For the extended Lorentz group this coincides with the usual notion of co- 
variant derivative, and in the case of electrodynamics it reduced to the well 
known recipe of introducing gauge invariant interaction terms. 


(*) Urryama (7) has proved that the free Lagrangian of the field can depend only 
on the combination of U and its derivatives in the r.h.s. of Eq. (3.13). It is evident 
that our Lagrangian (3.12) differs only by a divergence from the square of this expres- 
sion. The latter would however lead to second order equations. 

(**) Cf. Note added in proof on p. 751. 
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As far as the physical existence of the U-fields is concerned, it must be 
stressed that one cannot expect to find real particles corresponding to each 
of them, as not all of the considered invariance properties are strictly universal. 
For example isotopic rotation invariance is true only as long as electromagnetic 
effects are being negligible, and so one may doubt the real existence of the 
Yang-Mills field. On the other hand the conservation laws of baryonic and 
leptonic charge and the associated invariance properties of the fields are pos- 
sibly rigorous laws, and should therefore lead to observable U-fields, coupled 
respectively to the baryons and leptons and possibly to the electromagnetic 
field. These two fields may in some respects resemble the Z-field, which ac- 
cording to SCHWINGER (5) is responsible for the weak interactions. These 
questions will be investigated in Part ITI. 

From the point of view of renormalizability, the U-fields lead in general 
to non-renormalizable non-linear theories, which are difficult to investigate 
with the usual methods. However, for commutative groups, the difficulties 
will not be greater than in quantum electrodynamics. This question will also 
be investigated separately in Part III. 


* oF 


The author wishes to express his gratitude to Professors D. I. BLOHINGEY 
and N. N. BoGoLsuBoy, Directors of the Joint Institute and Theoretical Labo- 
ratory, respectively, for the kind hospitality extended to him. It is a pleasure 
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RIASSUNTO (*) 


Il presente è il primo articolo di una serie in cui si esaminano le conseguenze del- 
l’invarianza della teoria quantistica dei campi rispetto ad un gruppo di parametri Î 
di Lie. La parte I tratta le conseguenze dell’invarianza di un lagrangiano generale 
rispetto a un tale gruppo. Si dimostra che l’invarianza conduce alla conservazione di f 
«operatori di densità di corrente », i corrispondenti operatori di carica essendo una 
rappresentazione del gruppo di Lie nello spazio dei vettori di stato. L’invarianza rispetto 
alle trasformazioni estese, cioè trasformazioni in cui i parametri sono funzioni dei 
punti dello spazio-tempo, porta alla necessità di introdurre f campi vettoriali non 
commutanti. Si determina la struttura generale del lagrangiano corrispondente a tali 
campi e si discute la loro entità fisica. La parte II conterrà una simile indagine nel 
quadro di una teoria della matrice S, e la parte III sarà dedicata a una classificazione 
dettagliata di possibili gruppi di trasformazioni e la corrispondente classificazione dei 
campi di particelle e delle loro interazioni. 


(*) Traduzione a cura della Redazione. 
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On K--Hydrogen Interactions. 
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G. DASCOLA, 8. MORA 


Istituto di Fisica dell Universita - Parma 


(ricevuto il 27 Novembre 1958) 


Summary. — Experimental data are given on the interaction of K -mesons 
with hydrogen nuclei. These emulsion data do now agree with the new 
bubble chamber data. They indicate an essentially isotropic angular 
distribution, both for scattering and for absorption events, a slight decrease 
of the elastic scattering cross-section with energy and a larger one of 
the absorption cross-section. 


The purpose of the present experiment was to collect further data on 
K=-Hydrogen interactions in nuclear emulsion, applying fixed criteria for their 
selection and analysis and uniformly correcting data to reduce to a minimum 
systematic errors. The interactions studied are elastic scatterings (K7+p — 
— K-+p) and absorptions (K~+p > &*+77). 


(*) At present at University of California Radiation Laboratory, Berkeley. 
(*) Now at the Brookhaven National Laboratory. 
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Four emulsion stacks, each composed of 110 pellicles 3.5 in. x 5.5 in. x 600 um 
were exposed at the Berkeley Bevatron to the separated Barkas-Heckman 
K~ beam. The K~ mesons of this beam had a mean energy of 80 MeV. 

The scanning procedure consisted in identifying K~ tracks at 3mm from 
the front edge where they had a grain density of about 2.5 times minimum 
and then in following them through the stack until an interaction in flight, 
a decay or stop in flight, or until the particle comes to rest. A total of 
300 K=-tracks has been followed. 

From those events in which the K~ interacted in flight, the K~-H elastic 


a 


DT 


scatterings or absorptions were recognized by checking carefully angular cor- 
relations and consistency of ranges, with the energy balance. 

In order to minimize observation biases the following selection criteria 
have been rigidly applied: 


a) For scattering events (K” +p > K~+p): 


1) All events in which the range of the proton is < 5 um irrespectively 
of the scattering angle have been rejected. The angular cut-off of elastic 
hydrogen scatterings is therefore a function of the energy of the K-particle. 
The reason for such a cut-off was to reject a type of event for which the scan- 
ning efficiency and the accuracy of angular measurements is low. 


2) Events occurring less than 10 um from either surface of the un- 


processed emulsion have also been rejected. 


3) If the energy of the KT at the point of scattering was smaller than 
10 MeV, only those events in which the K~ ended in a capture star with at 
least one heavily ionizing prong were considered. This restriction was intro- 
duced in order to eliminate the considerable number of K~ two prongs capture 
stars at rest which, for a particular range and direction of the outgoing par- 
ticles may well simulate a very low energy K-p scattering. 


b) For absorption events (K-+p->Z*+x?7). 


4) Only events in which the laboratory projected angle between the 
x direction and the incident K~ is > 10° have been considered. In some 
vases the thin pion track was overlooked during the scanning and hence the 
event was classified as « single scattering » of a K-particle. Therefore all the 
events classified as « single scattering » by the scanner have been re-examined 
and a careful search for a possible z-track has been made in the angular region 
suggested by kinematical considerations. 


5) Absorptions occurring at a K7 energy of less than 10 MeV have 
been rejected owing to the difficulty of discriminating them from absorption 
stars at rest having one thin and one heavily ionizing prong and accidentally 
fitting the kinematics of an absorption in flight by hydrogen. 
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6) Events occurring at a distance of < 50 um from the surface of 
the unprocessed emulsions have been rejected. 


For selection criteria 2) and 6) the. correction is straightforward; for 3) 
the correction coefficient is given by the fraction of the K stars leading to 
capture stars with NV, >1. This fraction has been taken equal to 22% from 
rete (2): : 

The correction for the selection criterion 4) has been applied assuming the 
angular distribution in c.m.s. of the absorption events to be isotropic. The 
observed frequencies of absorption have finally been increased by 10% in 
order to take into account the loss of fast pions. 

To the data collected during this experiment we have also added those 
previously obtained by the GOttingen group (1), after having applied the se- 
lection and correction criteria given above. We are then left with a total of 
135 scattering and 47 absorption events. 

In Table I we have listed the track lengths followed in different energy 
intervals and the corrected number of events of the two types. 

The observed number of events are given in brackets. 


James IP 
Energy | Track Number of events Cross-section (mb) 
interval et ———t ee e dl di 
- | followed | | | 
| (MeV) Imola Elastic Inelastic Elastic Inelastic 
MOZIO 294 421 (3) | cut-offregion 58.5-L33.8 = 


10—20| 542 | 10.34 (10) 8.65 (6.5) | 59.4418.8 


56.1411.6 Lane 
2030 | 80.6 | 13.97 (13.6) | 11.35 © (855) | 5394147 


30—40| 107.2 |22.76 (22) | 6.65 (5) | 66.0-+14.1 


40—50| 134.0 |20.68 (20) | 13.30 (10) | 48.04-10.7 


1 50—60)) 160.0 | 212012055) | 11:97 (9) 9) 41-226 9.1 


| | | 
60—70| 177.0 |, 24.81  (24)-| 7.98! (6) | 43.64 9.9 


70-80| 143.5 | 20.69 (20) | 2.00 (1.5) | 44.9+10.0 


| | 
80—90| 14.57 | 2.07 (2) | 0.67 (0.5) | 442+31.3 
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(1) W. ALLes, N. N. Biswas, M. CECCARELLI and J. Crussarp: Nuovo Cimento, 
6, 571 (1957). 
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In Fig. 1 the c.m. angle and energy of the K7 for each elastic scattering 


on hydrogen is represented by 

to the 5um cut-off mentioned 

above. 
Fig. 2. 


results 


shows the correspond- 
ing for K° absorption. 
The plotted those 
between the directions of motion 
of the K° and the 
c.m.s. The sign of the X is shown; 
open correspond to the 
events leading to & 


angles are 


x in the 
circles 
»’s of undeter- 
mined sign. 

Fig. 3 and 4 show the cross- 
section versus energy for absorp- 
tion and for scattering events. 

In the same figures are also 
shown the curves calculated by 
DALITZ (*) fitting the low energy 


(2) KR. EN DATI Report of 
the VIII Annual Conference on 
High Energy Physics, Geneva (1958), 
D'190: 


one point. 
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bubble chamber point (*) by means of the zero range approximation treat- 
ment proposed by JACKSON et al. (4) for a purely S state interaction. 
EMERSO and er 

the c.m.s. angular distri- n K+P>+K+P À 
butions of elastic scatter- 207 10-30 MeV 20 
ing for three K~ energy 
intervals and in Fig. 8 

the angular distribution 

of absorption events of 

all energies are repre- 

sented. These results in- Fig. 5. 
dicate: 


= 
(©) 


1) A slight decrease 


of the scattering cross- 
section and a stronger 
than 1/v decrease of the 
absorption cross-section 
with increasing K° e- 


Number of events 


= 


nergy; in good agree- 
ment with the calculated 
curves. 


nora 


2) An essential isotropy of the scattering and absorption angular distri- 
butions; above 60 MeV one may notice a certain forward peaking in the scat- 
tering angular distribution. 


These features seem to be also present in the collected data of other emul- 
sion groups, and in those of the Berkeley hydrogen bubble chamber (2). 


(3) Report of the VIII Annual Conference on High Energy Physics, Geneva (1958), 


D'MITONLTS 
(4) J. D. JACKSON, D. G. RAVENHALL and H. W. Wvrp: Nuovo Cimento, 9, 834 
(1958). 


RIASSUNTO 


Si presentano risultati sperimentali sull’interazione di mesoni K negativi con 
nuclei di idrogeno. Questi dati, ottenuti con emulsioni nucleari, sono in accordo con 
i nuovi dati di camere a bolle. Essi indicano una distribuzione angolare essenzialmente 
isotropa sia per eventi di diffusione elastica che di assorbimento, una leggera dimi- 
nuzione coll’energia della sezione d’urto di diffusione ed una più forte diminuzione della. 
sezione d’urto d’assorbimento. 
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The Electric Dipole Moment of Elementary Particles. 


R. L. GARWIN 
IBM Watson Lab. and Dept. of Physics, Columbia University - New York 


L. M. LEDERMAN (*) 


CERN - Geneva 


(ricevuto il 12 Dicembre 1958) 


Summary. — The detection of electric dipole moments (EDM) of ele- 
mentary particles is discussed and a new limit to the size ot the electron 
EDM is presented (+ e:1.5:10-15em). This result is used in a discussion 
of recent observations of positronium annihilation. 


1. — Introduction. 


It is still of considerable interest in the testing of time-reversal invariance 
and parity conservation to attempt to detect the electric dipole moment (EDM) 
of an elementary particle. The existence of an EDM implies non-invariance 
under time reversal and non-conservation of parity (*°) and its absence under 


(*) J. S. Guggenheim Fellow, 1958-59. Permanent address: Department of Physics, 
Columbia University, New York. 
(1) E. M. Purcern and N. F. Ramsey: Phys. Rev., 78, 807 (1950); J. SMITH, 
E. M. PurceLL and N. F. Ramsey: Phys. Rev., 108, 120 (1957). 
(2) L. Lanpau: Zu. Eksper. Teor. Fiz., 32, 405 (1957); Nucl. Phys., 3, 405 (1957). 
(3) T. D. LEE and C. N. YANG: Elementary Particles and Weak Interactions, BNL 443 
) 


CEO LW LOST) peal a 

(4) N. BF. RAMSEN Phys. Rev., 109, 225 (1958). 

(5) D. BerLEy, R. L. Garwin, G. GipaL and L. M. LEDERMAN: Phys. Rev. Lett., 
1, 144 (1958). 
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a severe test is a sensitive indication of invariance of the particle under space 
reflection or time reversal (5). 

SMITH, PURCELL and RAMSEY (1) in an atomic beam resonance experiment 
have set an upper limit to the electric dipole moment of the neutron. Their 
result is D,< 107% em where D, multiplied by the electric charge gives the 
dipole moment. Since we know that parity is not conserved, we may regard 
this measurement as showing that |F|?-/@|?<3-10-!, where G is the proba- 
bility amplitude for parity mixing (which may be large), and F is the relative 
amplitude non-invariant under time reversal (*). 

Motivated by the desire to find some difference between the muon and the 
electron other than the rest mass, a search has been made for a muon electric 
dipole moment with the following result (°): 


(1) D 10m. 


A discussion (7) of the effects of a possible electron electric dipole moment 
in the spectroscopy of the hydrogen atom led to the results (this applies equally 
well to the proton): 


(1) D, < 10-* em. 


Recent reports (7) of the observation of effect proportional to o-E in posi- 
tronium have led us to a further examination of this question. We call at- 
tention in this note to three points: 


i) Many proposed experiments designed to detect the electric dipole 
moment of charged particles are based on the incorrect notion that it is, in 
practice, possible to apply an electric field to such a particle without the con- 
sequent free acceleration implied by Newton’s second law. 


ii) There exist in the literature, data on the behaviour of free electrons 
which enable a more sensitive limit than (2) to be found for the electron EDM. 
This type of experiment may be improved by an additional order of magnitude. 


iii) The limit thus set makes it extremely unlikely that intrinsic electric 
moments are responsible for the «o-H » effect observed in positronium. 


(6) I. ZeL'povicn: Journ. Exp. Theor. Phys., 6, 1148 (1958). Translation. This 
illuminating paper discusses the connection between electric moments and asym- 
metric decay. 

(*) Alternatively, |#|?<3-10-¥, where # is that amplitude which introduces both 
parity and time-reversal mixing. 

(7) G. FEINBERG: Phys. Rev. (to be published). 

(8) F. OBENSHAN and L. PAGE: Phys. Rev. (to be published). 
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2. — Analysis of some incorrect experiments. 


Many experiments have been performed or proposed for the measurement 
of the EDM of charged particles. These generally run into a difficulty which 
does not exist for neutral particles (1). The problem caused by the charge 
is simply that to measure an EDM one must apply an electric field at the par- 
ticle in order to generate an interaction energy for the parameter being ob- 
served. A charged particle cannot exist unaccelerated in an electric field. 
However, charged particles in matter have time-averaged motions which are 
either rest or uniform (drift) velocity with respect to the applied electric field. 
Such arrangements have essentially no sensitivity to EDM. The essential 
reason is that all forces acting on the particle in matter are electrical. For 
example, the mobility of ions in liquids is well known to be of the order ~ some 
em/s per kV/em of applied field. Charged particles, e.g. ions in the liquid, 
move at constant velocity in zero average field since the electric drag exerted 
by the medium exactly equals the applied field. Thus, for example, proposals 
to observe a shift in the resonant frequency in nuclear magnetic resonance 
carried out in the presence of a strong electric field, i.e. observing the para- 
meter d in (u:H+do:Æ), fail because the nucleus shifts to the new minimum 
in the total electric potential which it experiences. Quantum mechanical con- 
siderations do not alter the conclusion that the charged particle must be acce- 
lerated by the applied field in order to detect the EDM (*). 


3. — EDM of the free electron. 


It happens that an experiment has been performed for quite a different 
purpose which has considerable sensitivity to the electron EDM. This is the 
anomalous g-factor experiment of PIDD, LOUISELL and CRANE (*). In this 
experiment electrons polarized normal to a uniform magnetic field are caused 
to execute a large number (~ 107) of orbital revolutions in a magnetic field. 
The vertical plane precession induced by the supposed EDM builds up only 
for a single half-period of the anomalous moment precession before decreasing 
again, since thereafter the spin points in the opposite sense in the orbit plane. 
However, no significant reduction in polarization is observed thereby setting 


(*) W. PAUL however notes that a second order effect may be observed in a strong 


electric field gradient. In this case the resultant electric field is proportional to the 
EDM times the field gradient. 


(*) R. W. Pipp: Bull. Am. Phys. Soc., 7, 364 (1958) and private communication. 
See also H. R. CRANE, R. W. Prpp and W. H. LouISELL: Phys. Rev., 91, 475 (1953). 
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a limit on the electric dipole moment equal to that of the anomalous magnetic 


moment (v/c&1): 
CAE 
Mes 
° = 2x \me) ? 


De 10 em. 


A quite different calculation can be made using the results of this experiment. 
An electric dipole moment of magnitude e times the anomalous magnetic mo- 
ment would change the anomalous precession period by a factor ~1-+ e2/2. 
This can be seen by treating the EDM precession in the electric field (v/c) x Hy 
as if a rotating magnetic field of magnitude (xe/2x)(v/c) H, were applied and 
EDM= 0. In the CS rotating with the cyclotron frequency we observe two 
stationary orthogonal fields: (ex/2x)(v/c)H, and H, — (w/y) = (x/2x)H,. Since 
this experiment agrees with measurements of the anomalous moment of the 
bound electron (7) and with theory to within at least 10~° of the anomalous 
moment, we find immediately that 


(3) D, < 1.5-10-5 em. 


A further increase in sensitivity can be obtained by à direct analysis by Mott 
scattering of the out-of-plane polarization for electrons which have been stored 
in the magnetic field for + and 3 periods of the anomalous precession. A measu- 
rement to 1° would yield a sensitivity D, ~10-* cm. 


4. — o È effect in positronium. 


In this experiment, an effect is observed when positronium is formed by 
polarized positrons in a 15 kV/cm electric field (§). The effect (a change in 
the detected rate of 2-quantum annihilation) has the symmetry o:-H. We 
only comment here on the impossibility of this effect being due to intrinsic 
EDM of the elementary particles as limited by (3) (we invoke the CPT theorem 
here only to similarly restrict the EDM of the positron). The argument is 
independent of the details of the experiment and is made as follows. The 
maximum interaction energy produced by the electric field is, using (3) 


eD,E < (5:10-%)(10-5)(50) = 3-10-88 erg. 


The corresponding angular frequency of the free oscillations between eigen- 
states of the perturbed system is thus 


E 
gag. 
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Since this corresponds to a period considerably longer than the lifetime of 
triplet positronium (~ 10-78) the maximum amplitude of such a presumed 
intrinsic EDM effect is 


< (10-7)(4-104) = 4-10-* 


which is far from the 5% effect observed. This is certainly an absolute max- 
imum, since detailed calculation of matrix elements, symmetry conside- 
rations, ete., can only reduce the effect. Thus, if the effect should turn out 
to actually have the symmetry o-E (and not bea p-E effect masquerading as 
c-E because of the o-p correlation in weak interactions (5)) it would have 
to be attributed to the interactions involved in positronium formation and 
decay. 


One of us (LML) would like to thank CERN for hospitality, and I. 
HALPERN and A. PETERMAN for useful discussions. 


RIASSUNTO (*) 


Si discute la rivelazione dei momenti di dipolo elettrico (EDM) delle particelle 
elementari e si ottiene un nuovo limite della grandezza dell EDM dell’elettrone 
(< e:1.5:10- em). Si usa tale risultato nella discussione di recenti osservazioni del- 
Vannichilamento del positonio. 


(*) Traduzione a cura della Redazione. 
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(ricevuto il 12 Dicembre 1958) 


Summary. — The limitation of the «dipole effect» (electron-positron 
screening) as a method for determining the energy of electron pairs in 
nuclear emulsion is discussed. It is then shown that the multiple scattering, 
which is mainly responsible for the separation of the trajectories at a 
great distance from the origin (~ 1mm) can be evaluated from the 
measurement of the « geometrical effect». This phenomenon is studied 
experimentally on electron pairs of energy between 0.6 and 1 GeV, for 
trajectory separations up to 0.6 um. 


1. — Introduction: dipole and geometrical effects. 


PERKINS first observed (1) that the grain density on the initial track of 
high energy electron pairs is less than two times the value of a plateau track 
and increases slowly to reach this value. He attributed this phenomenon to 
a screening effect of the two electrons (Dipole effect). The calculations by 
several theoreticians (*>) of the ionization of such a dipole were found to be 
in agreement with the experimental observations (9-8). The grain density 


(*) On leave of absence from «Laboratoire de recherches nucléaires E.P.U.L. », 
Lausanne. 
(1) D. H. Perkins: Phil. Mag., 46, 1146 (1955). 
(2) A. E. CupaKov: Compt. Rend. Acad. Sci. U.R.S.S., 19, 651 (1955). 
G. YEKUTIELI: Nuovo Cimento, 5, 1381 (1957). 
I. Miro and H. Ezawa: Progr. Theor. Phys., 18, 437 (1957). 
G. H. BURKHARDT: Nuovo Cimento, 9, 375 (1958). 
W. WoLter and M. Mrrsowicz: Nuovo Cimento, 4, 648 (1956). 
J. IDAWARE: Phil. Mag., 3, 680 (1958). 
A. A. VARFOLOMEEV: Private communication. 
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starting from the value zero at the origin rises to about 1.8 times the plateau 
density at a point where the separation between the two trajectories is about 
3:10 um. From this point on the two electrons can be considered as inde- 
pendent for the ionization processes. 

However, it was shown (#1) that another phenomenon appears then, 
following the dipole effect. The blob density on the electron pair track, at 
the point where the dipole effect stops, is only 1.5 times the plateau density 
and increnses to two times, even though the ionization of the two electrons 
is now twice the plateau. The blob density is much more sensitive to this new 
phenomenon than the grain density. This change of the track structure was 
called geometrical effect, and is related to the fact that two electrons can cross 
the same grains. A part of the energy released in the grains is then lost for 
the photographie process. The measurement of the geometrical effect on 7 
cosmic electron pairs of energies between 10 and 100 GeV gave results in agree- 
ment with the prediction of a three dimensional model of track formation (°°). 


2. — Measurement of the pair aperture through the dipole effect. 


It has been proposed to use the dipole effect as a method for the deter- 
mination of the angular aperture of high energy electron pairs. The Borsel- 
lino relation (1!) can then give a most probable value for the energy of the 
photon. However, it was pointed out in a previous work (1!) that an angular 
measurement cannot be performed at a distance from the origin of the pair 
greater than 1 mm because after such a length the separation between the two 
electron trajectories is due mostly to multiple scattering. This means that 
the use of the dipole effect for angle determination will be useless for pairs 
above an energy which can be caiculated in the following way. 

We can admit that an initial electron pair track can be recorded, provided 
its grain density is greater than 0.3 times the plateau density. According to 
the CupAKOV dipole theory (2) this corresponds to a trajectory separation of 
about 2-10 um. In order that the ionization variation be measured with 
enough accuracy, it is necessary to record at least 200 grains on the track 
segment starting from the point where the separation is 2:10 ym, and ending 
1 mm after the (unknown) origin. In the case of normal G-5 emulsions, with 
about 20 blobs/100 um for a plateau track, the integration of the Cudakov 
relation shows that this requirement is not fulfilled for aperture angles of 


3 


(9) R. WeiLz, M. GaAïILLOUD and P. RosseLET: Nuovo Cimento, 6, 413 (1957) 
6, 1430 (1957). 

(10) R. WEILL: Helv. Phys. Acta, 31, 641 (1958). 

(1!) A. BORSELLINO: Phys. Rev., 89, 1023 (1953). 
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less than 10 radians. This means that the angle measurement by means of 
the dipole effect is impossible for photons of more than 10% GeV. The hyper- 
sensitized Bogomolov emulsions with 105 grains/100 um at plateau (12) should 
give the possibility of reaching 106 GeV.) 


3. — The use of the geometrical effect for the measurement of the multiple 
scattering. 


It would then be natural to use the dipole effect for the measurement of 
the multiple scattering (see Sect. 4). However, this can be done more success: 
fully with the help of the geometrical effect, which occurs for such cases in 
the region where the separation is due mostly to multiple scattering. The 
measurement of this effect will then provide a measurement of the multiple 
scattering. Suppose that the variation of the track structure (blob density 
and/or mean gap length) as a function of the separation of the trajectories 
is known and represented by a «calibration curve ». We can then measure 
the length L over which the blob density changes by a given amount for a 
pair of unknown energy. The calibration curve will give us the variation AY 
of the lateral separation of the trajectories over the length Z. From these 
two data, L and AY, we have to estimate the energy of the photon K=#,+E_. 
The relative projected lateral displacement due to multiple scattering can 
be considered as determined by a Gaussian distribution. The probability of 
having a separation between AY and AY+dY will be 


K, AY?K? 
[BB ea ee oa 
V/2nb Lt 26208 


where 
Chih Tee eect Fs (AE 
b is the multiple scattering constant. The principle of « maximum likelihood » 


tells us that the best value X* we can attribute to K, in our measurement is 
that for which: 


dlog P … O 
char 
This gives: 
es AI 
TK —= AY . 


(12) C. BocomoLov: Congrès International de Photographie Corpusculaire, Montreal 
(August 1958). 
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4. — The uncertainties and the limits of the method. 


The first uncertainty in the determination of A> arises from the error in 
the blob number, caused by statistical fluctuations. This corresponds in some 
way to the reading error (noise) which occurs in normal multiple scattering 
measurements. It can be reduced to less than 5% for pairs of more than 20 GeV 
for which the available length is greater than 1 mm. 

The «statistical inaccuracy » on KY will be defined by the two intervals 
A'K, and A’K, such that the likelihood function P drops by a factor of two 
for K,= K* {4 A'K, or K*—A'K,. For one cell the accuracy will be cha- 
racterized by K* i: For photons of more than 10? GeV the rapid deve- 
lopment of the initiated shower (appearing mostly as pseudo or real tridents 
on the initial track of the pair) will leave only about 1 em for the measure- 
ments. But on such a length the geometrical effect can be studied with a 
maximum of 3 cells (not more, in order to keep the blob density fluctuations 
low enough). The statistical accuracy will then be given by ÆŸ* es for two 
cells the accuracy will be given by HE The statistical accuracy must 
of course be combined with the statistical error on blob counts, and the number 
of cells used (1, 2 or 3) should be adjusted to the available length in à way 
which minimizes the total error. 

The proposed method is free of the stage and grain noise inherent in the 


normal multiple scattering measurements. 


Limit of validity of the energy determination on electron pairs 


Multiple 
scattering 
measurement 
direct |by geometrical 


Relative multiple 


scaltering measurement by dipole effect 


Angle 
measurement 


geom_|dipole 


1 103 10° 
Fig. 1. — Limits of application of the different methods of electron pair energy measu- 
rement. 


This method, however, is not applicable to pairs of more than 104 GeV. 
For this energy the geometrical effect will begin 1 em after the origin. But 
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then the dipole effect can be used for the measurement of the multiple scat- 
tering. However, it is easy to calculate that in this case, admitting the neces- 
sity of recording at least 200 grains, the dipole method could not be useful 
for angles of less than 10-® radians (X,= about 10° GeV). The limit will be 
about 1071 radians with Bogomolov emulsions. This value seems to be the 
highest electron pair energy measurement directly accessible with nuclear 
emulsions. We give in Fig. 1 the limits of application of the different methods 
of electron pair energy measurement. 


5. — Measurement of the geometrical effect (calibration curve). 


We have measured the geometrical effect on 76 electron pairs recorded 
in Ilford G-5 emulsions exposed to the photon beam of the Cornell synchrotron. 
The measurements were performed on electron pairs for which the initial un- 
resolvable track was at least 90 um in length (pairs of about 600 to 1000 MeV). 
For such an energy the dipole effect is effective for only about 1 um. The 
tracks were divided into cells of 30 um and the blob count was done for each 
cell. The projected separation of the two electron trajectories at the middle 
of each cell was deduced from the measured projected opening angle of each 
pair. (In this photon energy range the influence of the multiple scattering is 
negligible (#°).) We have averaged 
the individual measurements, group- 
ing the data in classes according to 
the trajectory separation. (The width 
of a class is 0.1 um.) Fig. 2 shows 
the variation of the blob density of 
electron pairs relative to the plateau 
track, as a function of Y, the pro- 0 01 02 2603 04 05 
jected separation of the trajectories. 


Fig. 2. — g* (blob density relative to pla- 


The errors quoted are the r.m.s. de- teau) as a function of the projected sepa- 
viations. We have also reported on ration Y um of the trajectories. V is the 
this graph the point obtained by average of Varfolomeev’s measurements 

ae i irs (8 a A 
averaging the measurements done on 9 cosmic pairs (8). Errors are r.m.s. 


ue deviations. 
by VARFOLOMEEV (°) at the beginning 


of the geometrical effect on 9 pairs 
whose energy range was from 10? to 3-10? GeV. These data were obtained 
in NIKFI R emulsions, whose characteristics, before and after processing, are 
similar to those of the G-5 used in our work. 

The geometrical effect starts with a blob density relative to the plateau 
of 1.55. This is in fairly good agreement with the value predicted by our three 
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dimensional model (°°) (1.55 to 1.60). Calibration against plateau tracks was 
done in all the plates (same batch 
and processing) in all the regions 
where the measurements were per- 
formed. All the measurements were 
taken in a layer situated 50 um from 
the surfaces of the emulsions. 

We give also (as an indication) 
the variation of the mean gap length 
as a function of Y, as deduced 


Fig. 3. G* (mean gap length relative to 
plateau) as a function of the projected 


separation Y um of the trajectories. Errors from the measurements of 50 pairs 
are r.m.s. deviations. (Fig. 3). 


6. — Influence of the processing on the calibration curve. 


The emulsions used for this work were Ilford G-5, with the following 
y) too} 
characteristics : 


Plateau track: Blob density: g = 20/100 um 
Mean gap length: G = (4.35 + 0.08) um 
Mean track diameter: D = (0.71 + 0.02) um. 


Let us call d, the mean grain diameter in the unprocessed emulsion. We 
define, according to O’CEALLAIGH (1%), a factor f which can characterize the 
processing by: 

D = doi + f). 

The three dimensional model enables us to calculate the influence of the 

processing on the calibration curve. For two different processings characterized 


by fi, and f,, the ratio of the relative blob densities will be according to the 
relation (5) of reference (1°): 


gig pal, CS) 


The definitions of the symbols are given in reference (1). Using the fact 
that P.a <2:Pia, and the experimental value of P,, = 0.09 we obtain: 


Gilg: = exp [— 0.09(f1 — fa)] - 


(13) C. O’CEALLaIGH: Suppl. Nuovo Cimento, 12, 412 (1954). 
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For a given standard processing, f varies less than 5%. This means that 
the calibration curve will change by an amount less than 1%. For extreme 
cases, for which f goes from 1.36 to 1.79 the calibration curve would be changed 
by less than 3%. The three dimensional model predicts that the calibration 
curve for the mean gap length should be independent of the processing. 


‘ 


* Ok OK 


We wish to thank Professors G. Cocconi and J. OREAR for many helpful 
discussions and Mrs. L. VAN Nest for her help in the measurements. 


RUTASISTUENELZIONO) 


Si discutono i limiti dell’efficacia dell’« effetto di dipolo » (schermo elettrone-positone) 
come mezzo per determinare in emulsione nucleare Venergia di coppie di elettroni. 
Si mostra quindi come lo scattering multiplo, che è la causa principale della separa- 
zione delle traiettorie a grande distanza dall’origine (~ 1mm) possa essere calcolato 
dalla misura dell’effetto «geometrico ». Si studia sperimentalmente il fenomeno su 
coppie di elettroni di energia fra 0.6 e 1 GeV per separazioni di traiettorie fino a 0.6 um. 


(*) Traduzione a cura della Redazione. 
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On the Decay of Xenon 155. 


S. Joa, H. G. DEVARE, M. NARAYANA Rao and R. M. SINGRU 


Tata Institute of Fundamental Research - Bombay 


(ricevuto il 17 Dicembre 1958) 


Summary. — The y-rays from 1#%Xe have been studied in a scintillation 
spectrometer. The y-rays have the energy of 80 keV, 160 keV, 300 keV 
and 380 keV with the relative intensity 1, 1.4:10-?, 8.4-10-4 and 4.3-10-4. 
By studying the 80 keV y-rays in coincidence, the relative intensity of 
the cascade and the cross-over transition from the 160 keV state has 
been found. Taking this into account, the log ft values for the f-tran- 
sition to the 80 keV, 160 keV and 380 keV states have been calculated 
to be 5.6, 7 and 5.7 respectively. 


1. — Introduction. 


13*Xe has been known to decay from its ground state (spin 3+) with a half- 
life of 5.3 days emitting a single group of f--rays of maximum energy of 
346 keV (1). This decay leads to the first excited state of 130$ at 80 keV (42). 
The isomeric state of this isotope (spin a) decays by the emission of a 233 keV 
y-ray. Recent studies of !*3Ba (*°), the ground state of which is 3*, have shown 
that the electron capture decay of this isotope excites the levels in 13*Cs at 
80 keV, 160 keV, 381 keV and 437 keV. The levels at 80 keV, 160 keV and 


(1958). 


(5) S. D. Korcxt et al.: Bull. Inst. Nuclear Sciences « Boris Kidrich », 8, No. 144, 
(March 1958). 


() M. G. Stuart, D. C. Lu and F. M. CLIKEMAN: Bull. Am. Phys. Soc., 3, 208 
(1958). 
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381 keV are reached also in the Coulomb excitation of 1#%Cs. The spins and 
parities assigned to these states of Cs are 3*, 5* or 3*, 3* and 3* respectively, 
the ground state of 15*Cs being 3*. If the 160 keV state is 3*, one cannot 
understand the absence of a predominant electron capture decay of !**Ba to 
this state; on the other hand, whether this state is 3* or 3*, there should be 
a considerable 8 -branching to this state in the @~-decay of Xe. In order 
to investigate this point, the study of the y-rays of Xe was undertaken (*). 


2. — Experimental investigations. 


Them Xe" source 
for this investigation 
was obtained from HEART y-ray spectrum in a 
Harwell. The radio- scintillation spectrometer 
active gas was sealed 
in a glass tube 5cm 
long and 1 em in dia- 
meter. The Harwell 104 
catalogue of isotopes 
states that, as the gas 
is collected from the 
iodine fraction of the 
fission products, ap- 
proximately 0.5% of 


230 keV 


3 
xenon activity is due 10 
to 181Xe” (half-life 12 - 
days, y-ray energy = 
163 keV and 80 keV). S 
Bs 300 keV 
T= 
3 
LS) 
Fig. 1. — The spectrum ne 
of the y-rays for 1#Xe | 
and 131 Xe”. a) was taken 380 keV 


three days after the ar- 
rival of the source in 
Bombay. b) was taken 4 
days later. The 230 keV 
y-ray peak can be seen to 
have come down to 10 
nearly one fourth its pre- 
vious value. 


Energy (keV) 


(*) This study was also suggested in a private communication by Dr. B. CRASE- 
MANN of the University of Oregon, U.S.A. 


for) 
oO 
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The strength of the source was 200 millicuries to begin with. The tube was 
enclosed in a brass tube, one end of which was closed with a platinum foil 
(555 mg/cm?) in order to cut down the strong intensity 80 keV y-ray and the 
K X-ray of Cs. The spectrum of the y-rays was taken in a scintillation spectro- 
meter with the NaI(T1) crystal, 14 in. in diameter and 1 in. high. The activity 
was followed for three weeks. The spectrum of the y-rays is reproduced in 
Fig. 1. One can see that y-rays of energy 80 keV, 160 keV, 230 keV, 300 keV 
and 380 keV appear distinctly. It was found that the 230 keV y-ray decayed 
with a half-life of about 2.3 days and could, therefore, be attributed to #?Xe”, 
which is formed (7) in the decay of I. The y-rays of energy 80 keV, 160 keV, 
300 keV and 380 keV decayed with a half-life of about 5.3 days. It was con- 
cluded that at least a very large fraction of those y-rays originated from the 
same source and from !83Xe. 

When the activity had died down considerably, a spectrum of the y-rays 
of Xe was taken without the platinum foil covering the source. From this 
spectrum the intensities of the y-rays were calculated. Corrections were made 
for the efficiency of detection, and the differences with energies of the ratio 
of the photo-peak area and total area. The relative intensity of the 80 keV, 
160 keV, 300 keV and 380 keV y-rays was found to be about 1, 1.4-10-?, 
8.4-10-4 and 4.3-10-* respectively. 


233 keV 


Eg =48 keV 
133 
Va 


_ 


137keV- IE Si 
381 keV 3/2 


Eg=347 keV 
~ 98 Vo 
log f{=5:6 
À 160keV 5/,* 

80keV 5/2 * 


0 keV Yo 


Fig. 2. — The decay scheme of 138Xe and 1#Ba. 


It is known from the studies of **Ba that the 160 keV state decays by the 
emission of a y-ray of energy 160 keV and also by the cascade emission of two 
80 keV y-rays (***). In order to estimate the relative intensities of 8--branching 
it was necessary to find the intensity of the 80 keV y-ray emitted in cascade. 


(7) A. H, SNELL and F. PLEASONTON: Phys. Rev., 111, 1338 (1958). 
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For this purpose, a fast-slow coincidence unit with a resolving time of 1.2-10-7s 
was used. It was found that the number of 80 keV y-rays arising from the 
160 keV state was 0.0037 times that emitted from the 80keV state. 
Taking this into account, and neglecting the contribution from the 380 keV 
state, the £--branching ratios to the 80 keV state, 160 keV state, and 380 keV 
State were calculated to be 98%, 2% and 0.1% respectively. From the known 
maximum ~~ energy emitted in transition to the 80 keV state, the maximum 
8° energies involved in transition to the 160 keV state and the 380 keV state 
were calculated to be 287 keV and 48 keV respectively. From these data, the 
log ft values for the @ -transition to the 80 keV, 160 keV, and the 380 keV 
states were found to be 5.6, 7, and 5.7 respectively. This is illustrated in the 
decay scheme given in Fig. 2. 


3. — Diseussion. 


From the occurrence of the y-rays of energy 160 keV, 300 keV, and 380 keV 
in the decay of Xe in addition to the well-known y-ray of energy 80 keV, 
it is concluded that, as expected from the results of the studies of Ba, 
B"-decay of *Xe leads also to the levels in 13°Cs at 160 keV and at 380 keV. 
However, the intensity of the 6 --branching to the 160 keV state is rather low; 
consequently, the log ff value is much higher than what one would expect 
it to be if one assigns to the 160 keV state the spin and parity 3* or 3* as sug- 
gested from the study of the y-rays from *°Ba. It is worth noting in this 
connection that, as pointed out in the introduction, there is almost no electron 
capture decay of !**Ba (4*) leading to the 160 keV state of "Cs. One can 
see that this is consistent, so far as the decay of Ba is concerned, with the 
assignment of the spin and parity of only 3° to this state. The spin and parity 3° 
assigned to the 380 keV state are consistent with the log ft value found here 
in the B -decay of 1**Xe. 

The relative intensities of the 300 keV and the 380 keV y-rays, though 
subject to large errors in our measurements because of their low intensities, 
are 2:1 as found bi KOICKI et al. (*) in their studies of Ba. The cascade 
to cross-over transition from the 160 keV state has been found here to be 1:4. 
The intensity of the 220 keV y-ray emitted in the transition from the 380 keV 
state to the 160 keV state is rather low. It may be mentioned that there was 
a very low intensity y-ray of energy about 500 keV, the origin of which was 
not investigated because of its low intensity. Lastly, towards the end of these 
investigations, the results of SNELL et al. (9) became available. Although they 
confirm the existence of these y-rays, they found the intensity of the 300 keV 
y-ray so low that they could not make the half-life measurement. 


ih 
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Our grateful thanks are due to Dr. P. AMMIRAJU and Miss G. C. PRAMILA 
for their assistance in this work. 


RIASSUNTO (*) 


Per mezzo di uno spettrometro a scintillazione si sono studiati i raggi y del 1#Xe. 
I raggi y hanno energie di 80 keV, 160 keV, 300 keV e 380 keV con intensità relative 1, 
1.4-10-2, 8.4-10-4 e 4.3-10-4. Studiando i raggi y di 80 keV in coincidenza, si sono 
trovate l’intensità relativa della cascata e della transizione dallo stato 160 keV. Tenendo 
conto di queste determinazioni, si sono calcolati i log ft delle transizioni 6 agli stati 
80 keV, 160 keV e 380 keV trovando rispettivamente i valori 5.6, 7 e 5.7. 


(*) Traduzione a cura della Redazione. 
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v--Capture in Light Nuclei. 


E. G. BELTRAMETTI 
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Istituto Nazionale di Fisica Nucleare - Sottosezione di Genova 


L. A. RADICATI 


Istituto di Fisica dell Universita - Pisa 
Istituto Nazionale di Fisica Nucleare - Sezione di Pisa 


(ricevuto il 29 Dicembre 1958) 


Summary. — A shell model analysis of the u--capture by light nuclei 
without heavy particle emission, is given. Some experiments are sug- 
gested to measure the coupling constant and the relative contribution 
of the Gamow-Teller and Fermi interactions. A detailed study is presented 
of the u--capture in 3He, ‘Li and 160. 


1. — The purpose of this paper is to discuss the possibility of measuring 
the coupling constant g for u--meson capture and the relative importance of 
the Fermi and Gamow-Teller interactions. This problem has already been 
discussed for medium-heavy nuclei by LUYTEN and TOLHOEK (') who have 
used a simplified shell model description of the nucleus. However if one wants 
a more accurate evaluation of the coupling constant g and of the relative 
contribution x of the Gamow-Teller interaction one needs a situation in which 
the initial and final nuclear states are completely specified. 

This suggests to analyze a capture in a light nucleus, possibly without emis- 
sion of neutrons. The formation of the final nucleus in a bound state is not 
actually an unlikely process in a light nucleus. Indeed by considering the 
available volume in phase space one gets, for a nucleus of mass number A: 


wo _ 1.56 
VAI 


(1) A. H. ToLHoEK and J. R. LuytEen: Nucl. Phys., 3, 679 (1957). 
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where #, and w, are the probabilities per unit time for the two processes 


no (AZ) (Ay Ze) 
u-+ (4, Z) = (A—1, Z—1)+n-+y. 


In Section 2 we present a formula for the matrix element for the first pro- 
cess (without neutron emission) valid for transitions of any order. Indeed 
one expects that higher order transitions should be fairly important in p-- 
capture due to the high energy release. To calculate the nuclear matrix ele- 
ment we assume the validity of the shell model in its L-S coupling version. 
This seems a fairly good approximation for light nuclei and can at any rate 
be used as a starting point for a more refined intermediate coupling calculation. 
To simplify the treatment we shall assume that the proton which absorbs the 
wu -meson belongs to the unfilled shell. This seems a plausible assumption at 
least for those processes which do not lead to neutron emission. 

In Sections 3 and 4 this formula is used for discussing the capture rate 
by *Li, *He and !*0 which seem to be the simplest cases amongst the light 
nuclei. 

Other cases are briefly discussed in Section 5. 


2.— The operator causing the transition is a tensor operator of order j which 
can by written in the form 


en) A == y > CH hO, (0) ru) (i) 7 
î=1 y 

where H i he a i are tensor operators of order k, a and 1, operating respect- 
ively on the space, spin and isotopic-spin parts of the wave function. The mag- 
gnetic quantum numbers g, y vary from —j to +7 or —k to + res- 
pectively. k represents the order of the transition; a = 0 or 1 for Fermi and 
Gamow-Teller transitions respectively. 

By using standard methods we may write the reduced matrix element 
of the operator (1) in the form (2) 


(2) FIFO TD = VE] E°) Gr]: 
pan Th 
D Kala, ><a ja > KEL PT ITs IM) aes S's 
p 
i Ab dk’ 


*U(kKI,LL,; biL')U(a} SS,; #8")[(2j + 1)(2Z +1)(28 + (ese DIE 


(*) See, for a similar discussion in the case of y transitions: A. M. LANE and L. A. 
RapICATI: Proc. Phys. Soc., A 67, 167 (1954). 
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where |J) stands for 


[0-11 aLSTM,J) 


and represents a state with n —1 particles in the /-shell, one particle in the 
L-shell, coupled to orbital angular momentum L, spin S, isotopic spin 7, third 
component M,. J is the total angular momentum, and « stands for all the 
other necessary quantum numbers. In equation (2), ¢«|}«,> are the coeffi- 
cients of fractional parentage (*), AK is the product of a Clebsh-Gordon coef- 
ficient by a Racah function: 


K(abedef: 20) =D er CS U(ebef: ad) 


and 
U(abed; ef) = [(2e + 1)(2f + 1)}}W (abed; ef) 
abe 
d e f;is the Wigner nine-j symbol (*). The value of the one-particle re- 
ILE 


|A|4> is 1 or 3 for Fermi (a=0) and Gamow- 


duced matrix element <4 
Teller (a=1) transitions. 

To evaluate the relative importance of allowed and second forbidden tran- 
sitions we shall assume that the u--meson wave function is a constant Cover 
the nuclear volume. 

With obvious notations the one particle space operator may be written 
in the form 

ily = V4ngi* Cj,(qr) aS (0, ®)» 


where q is the neutrino momentum. Using an infinite potential well with radius 
2,54-10-1 em (appropriate to À — 6) one gets for transitions between states 
of the 1 p-shell, and for 4— 100 MeV/c: 


App) _ 


i GPA) — 


3. — The probability for formation of *He due to the capture of u--mesons 
by ‘Li (J=1-+), can be measured by detecting the @--activity of the target. 
If it should turn out that the transition to the excited state (J= 2+; 


(8) The p-shell fractional parentage coefficients have been tabulated H. A. JAHN 


and H. VAN WIERINGEN: Proc. Roy. Soc., A 209, 502 (1951). 
(4) See, for example, M. E. Rose: Hlementary Theory of Angular Momentum (New 
York, 1957), Ch. XI, or A. R. EpMmonDps: Angular Momentum in Quantum Mechanies 


(Princeton N.J., 1957), Ch. 6. 
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1.77 MeV) is unimportant, then the capture rate for the reaction 


u- + SLi > ‘He + v 


can be used to determine in the usual way the G.T. coupling constant. 
The states of “Li and He are, in L.S coupling, of the form. 


6Li ground state: 138,9 
‘He ground state: 3180) 


‘He first excited state: |°4D,>. 


One can therefore easily see that in L.S coupling there is no allowed transition 
to the first excited state and no forbidden transitions to the ground state. 

The relevant matrix elements, calculated under the assumptions that lead 
to equation (3), are, apart from the constant |¢|?, which will be consistentiy 
omitted : 


(a) Peed? 4:70 for the allowed transition to the ground 
state 
(b) |Mar|} = 138-102 for the second forbidden transition to the 


excited state. 


An intermediate coupling calculation (*) modifies (a) by only a few per- 
cents, (| M,,|? still vanishes), while (b) becomes even smaller (~ 4.9-10-*) due 
to interference effects, (| W,/?= 2.3-10-5). We see from the above discussion 
that we can neglect the probability for the formation of *He in its excited state. 

In order to determine the relative importance of the F and G.T. couplings 
one could use the reaction 


(5) u + 5He >v+ 3H. 


*H has no bound excited state and, as “He, undergoes B7-decay. Though this 
is presumably a very difficult experiment, the theoretical situation is very 
simple (**). 


() For the coupling parameter we have used the values a/K=1.35, K = —1.25 MeV 
and L/K=6 as suggested by InaLIS (Phys. Rev., 87, 915, (1952); Rev. Mod. Phys., 
25, 390 (1953)) and by FRENCH (University of Pittsburg: Technical Rep. n. IX (1958)), 
to whom we refer also for the notations. 

(”) The possibility of this experiment has been suggested to us by Professor M. GoLp- 
HABER (private communication) whom we wish to thank for some useful suggestions. 
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To evaluate the matrix element we have used for both *He and °H the 
space wave function suggested by IRVING (5) 


yp = exp[—eV7r, +72, +735], (1/e =1.089-10-! cm), 


where 7,3, 113, 7»g are the relative co-ordinates of the three nucleons. With 
this function the values of the two matrix elements are 


||? = 0.89, 


| Mon |? = 2.67. 


4. — The determination of g and x we have discussed depends upon the 
values of M,, and M, or, in others words, upon the choice of a particular nu- 
clear model. 

As in the analogous case of the ~-decay, one would like to have some 
other independent way for determining g and x in order to check the consis- 
tency of the assumptions made about 
the matrix elements. A case which 391 key = 
can be analyzed in some detail oc- 300kev-G 
curs in the y--capture by 10 lead-  313Kev-22 
ing to 16N. 16N has only four bound ground 2 
states whose spins and parities are 
reasonably well known (°) (see Fig. 1). 

All the transitions to the excited = = 4 | 
states contribute, through fast y-de- 60 

cay, to populate the ground state Fig. 1, 

which undergoes B--decay (t =7.35 8; 

E — 10.4 MeV). All these transitions are at least first forbidden. However, 
because of the higher charge of 10, the probability for formation of !*N is 
of the same order of magnitude as the probability for the formation of *He 
from SLi that has been discussed before. 

In Table I we have collected the values of the first forbidden matrix ele- 
ments for the different transitions, calculated with equation (2). For this eva- 
luation we have used the wave functions for *N calculated by ELLIOT and 
FLOWERS in intermediate coupling (7). The reduced single particle matrix 


(5) J. Irvine: Phil. Mag., 42, 338 (1951). 
(6) D. H. WILKINSON: Phys. Rev., 105, 686 (1957). 
(7) J. P. ELLIOTT and H. B. FLowERS: Proc. Roy. Soc., A 242, 57 (1957). 
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elements of the space operator H® (see equation (2)) have been evaluated 


E È : pei, 
using an infinite square well with radius a = 3.53 =10718 em. 


TABLE I. — 160 —16N first forbidden transitions. 


First order transition 
Final state 
16N k 
pel | Mon |? | Me? 
1— (391 keV) 5.18-10-2 1242-105 
O— (1138 keV) 5.54-10-2 0 
2— (ground) 11.3: 1072 0 
(3—) (300 keV) 0 0 


We have neglected higher order (3-rd order) transitions to the 2 — and 
(3 —) levels. A rough estimation indicates that they are unimportant. 

It is apparent from the above results that the measurement of the f-activity 
of the target can actually provide only the value of the G.T. coupling constant. 


5. — If we require that the capture processes should be observed by de- 
tecting the f-activity of the daughter nucleus, then it is difficult to find any 
other case that could be analysed in such a complete way as *He, ‘Li and 160. 
Usually there are too many accessible excited states whose spins and parities 
are unknown. 

The nuclei with a ground state of isotopic spin 7= è and third component 
M,=—2 could be reached from the ground state of nuclei with T=}, 
M,=—3}. They are B-active and have presumably very few excited states. 
Amongst the more promising candidates in this group we suggest 


°Li (Be +p- > *Li-+ v) and #0 (EN + u- > 0+ y), 


Another possibility would be to study the reactions which lead to the emis- 
sion of one neutron; for example: 


bil +p >°He4-n+,, 

*Be-- w= shi +n y. 
In both cases the final nucleus is f-active and has few excited states. 
This would probably provide a more feasible way for determining the Fermi 


coupling constant than the y--capture by *He that has been previously dis- 
cussed. A detailed study of these cases will be given in a forthcoming paper. 


wo 
a 
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Finally one could look for those transitions which lead to formation of the 
final nucleus in an excited state. These transitions could be observed by 
detecting the subsequent y-transition. 

Some examples have been studied in connection with the 190 + 16N tran- 
sition (see Table I), and a few other cases could probably be found. 


* KOOKS 


We wish to thank Professor M. GOLDHABER for a stimulating discussion, 
and Mr. G. FrorIo for performing some numerical calculations. 


RIASSUNTO 


Servendosi del modello a shell si analizza il processo di cattura dei mesoni u+, senza 
emissione di particelle pesanti, da parte di nuclei leggeri. Allo scopo di poter conoscere 
il contributo relativo dell’accoppiamento di Gamow-Teller e di Fermi, si suggeriscono 
alcune esperienze e si studia dettagliatamente la cattura u- in SLi, $He e 160. 
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Summary. — The equivalence of the pseudoscalar and pseudovector 
couplings of pseudoscalar pions to nucleons has been investigated here 
in a straightforward manner without making any use of a canonical 
transformation as has been done previously (1). PS-coupling is shown 
to be equivalent, as far as the S-matrix is concerned, to the usual p-wave 
and s-wave terms along with some non-local terms whose interpretation 
is not yet clear. 


There has been considerable work done to investigate the extent of the 
equivalence of pseudoscalar (PS) and pseudovector (PV) couplings of pseudo- 
scalar pions to nucleons (1). The usual way to do this has been to perform a 
canonical transformation on the Hamiltonian #,, of pseudoscalar coupling 
theory and approximating the transformed Hamiltonian in an appropriate way 
such as to get the PV-coupling term (?) i(G/2M)py“y,t,p(ep,ex") and the 
S wave scattering term (G2/2M)pyœ;. 

The purpose of this note is to show that a similar result can be obtained 
more directly without the use of a canonical transformation. In fact, starting 


(1) F. J. Dyson: Phys. Rev., 73, 929 (1948); K. M. Case: Phys. Rev., 76, 14 (1949); 
L. L. Fozpy: Phys. Rev., 84, 168 (1951); G. WENTZEL: Phys. Rev., 86, 802 (1952); 
L. L. FoLpy, R. K. OsBorx and J. M. BERGER: Phys. Rev., 87, 1061 (1952); J. V. LEPORE: 
Phys. Rev., 88, 750 (1952). 

(?) We shall follow the notation used in S. S. ScHwEBER, H. A. BETHE and 
F. DE HOoFFMANN: Mesons and fields (New York, 1955). Also # = ¢ = 1. 
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from the PS-PS Hamiltonian #,,, one can construct a « Hamiltonian» # 
which contains these two terms together with some other «non-local » terms, 
and which is equivalent to the original Hamiltonian as far as the scattering 
matrix S is concerned. 

We work in the interaction representation, which has not been used for 
these investigations except by LEPORE (*). The derivation of the new « Hamil- 
tonian » is as follows. | 

In the interaction representation, we have the following equations of motion 
for the state vector Y(t) and the field variables w(x) and ,(x) in the PS-PS 
theory (*), 


a id as 
(1) i =a | dw(pyst-py)x) P(t) , 
(a =t) 

(2) Dit x) pix) — 0, 
(3) (D + 1) 9, (x) = 0 
or 
(4) a i d'v(Pyst:Pp)(x) UE — 00) , 
where 

Vit) = U(t, — )P(— co) 
and 
(5) U(+ 00, — 00)=$, U(— 00, — co) =1. 


Substituting the values of y(x) and y(x) obtained from (2) in (1), we get (5) 


. QU(t, — co) LG À op SA, oy | 
Ge pen e ASA P YS NS) 


œuf 
iG [fore ai: Opal 7 
_— 32M dx ET (pysy"t- py) — Yy:y Tap Bae US) 


(3) J. V. LEPORE: Phys. Rev., 88, 750 (1952). 

(4) We shall use flat space-like surfaces instead of arbitrary ones for the sake of 
simplicity. 

(5) Only the limits of integration of the time component shall be specified, as 
the integration is always over the whole of allowed domain of space-components. 


Luni 
© 
a 


802 V. SINGH 


+ 


| G lo 
Ul, — 00)=1+ saz |? (ren a) METE 
n 
he 
+2, fare core) 86 — 
LA ui ae da: (rs: + 7) Ue 


G A ne QU(t, — co) 
mr: D°Yy5Y PY Gt 


— (es) 


È 
G 
| a af > ATI U(t, — c0)) . 
Substituting the value of 0U(t, — co)/et in (6) from (4), and using the 
commutation relations for y’s, we get 
- 0 
I G = Wa 
(7) U(t pri oo) =1 a 2M 20° Pyh ysT aM ant U(t, co) 


Ses) 


+ 


i È — 2 
= sar | ar Port U(t, — 


+ fal do ee 
os 
pia, vg) (yet ru oe iy ee 
tf 
G 
tam 


00) 


O - 
a cal (Pysyit- py UC co) : 
The last term can be integrated (5) and (7) then becomes 


(8) U(i', — 00) =1~ifare- x(x) 0 Be: 


=i OD) 


G la 
er (pysyt pyp)(x, YU, — co), 


(6) We make use of the adiabatic hypothesis in writing (8), that is we supply con- 
vergence factors (7) exp[—|¢|/7] and take the limit 7’> +co, after the integration 
has been performed. The adiabatic hypothesis is also used in writing equation (15). 

(7) F. J. Dyson: Phys. Rev., 82, 428 (1951); 83, 608, 1207 (1951); Proc. Roy. 

oc. (London), A 207, 395 (1951). 
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where 


iG _ 0 Gy "GA : 
(9) Fox, t) = om Ÿ? very aa + oM PPP? + es FE ,(X 3X1, t) d'A: , 


(10) Fn(x3x1,t ae (vst P(X)) ap Pp (Hr st) Palx) (7 yat P(%1 st) po Yo(r 1) 
Defining | 

di) Pw, — 00) = (1-57; ft 1) UW, — co, 

where 

(12) f(x, t) = (pysy?t-py)(x, 1°). 


The equation for V(#,— co) can be written as 


(13) We ita if ‘rFe(x, t) V(t, — 0) , 
where 

(14) Bex, t) — Hal, À sp [are fe, A} 
Now (cf. footnote (°)) 

(15) Po, — 00) = U(4- co, — co) = S. 


Therefore as long as we are interested only in the case t’—>-+ oo, 1.e. S-matrix, 
it does not make any difference whether we use U(t', — co) or V(t’, — co). 
Hence we can interpret A(x, t) as a « Hamiltonian » which would give the 
same results for scattering problems. 

We can write our new « Hamiltonian » d(x, t) as 


(16) F(x, t) = (x) + (x, t) 
with 

G G2 _ 
(17) Hex) = è PU Tru i. + op PPPs» 
and 
(18) AUX, t) = Hy(x,0)+ Hoo(x, t) +... , 


c 
[=] 
a 
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where 


(19) Fast) far, 


9 G3 
Pre Fe, (x) fx: t) f (x2, 1) a Tees X1,t)f(x2, 7) … etc. . 


Thus we see that Æ,(x, t) involves implicitly the field operators at more 
than two space points, and is in this sense «non-local». The local part 
FH, (x) of F(x, t) consists of two parts which give rise to P- and S-wave meson 
seattering respectively at low energies. The interpretation of the rest of the 
terms (the «non-local part ») is not any more clear here than in the other 
derivations using canonical transformation. 

The new quantity (x, t) which we call « Hamiltonian », has been shown 
to be equivalent to #,, only for scattering states and the two could give 
different results for the bound states. It appears therefore that #(x, t) cannot 
be derived from #,, using a canonical transformation. 


The author is thankful to SHRI B. M. UDGAONKAR, Dr. L. K. PANDIT and 
Dr. S. N. Biswas for helpful discussions. 


RIASSUNTO (*) 


L’equivalenza degli accoppiamenti pseudoscalare e pseudovettoriale dei pioni pseudo 
scalari ai nucleoni è stato qui esaminato in maniera diretta senza far uso di trasfor- 
mazione canonica, come fatto precedentemente (!). Si dimostra che l’accoppiamento PS 
è equivalente, per quanto riguarda la matrice S, agli usuali termini in onde p ed s as- 
sieme ad alcuni termini non locali la cui interpretazione non è ancora chiara. 


(*) Traduzione a cura della Redazione, 
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Energy-Band Structure of Lithium Atoms in the Diamond Lattice. 


F. BASSANI and V. CELLI 


Istituto di Fisica dell Universita - Pavia 


(ricevuto il 31 Dicembre 1958) 


Summary. — A hypothetical crystal of lithium atoms in the diamond 
lattice is considered. The value of the energies at the special points of 
the reduced zone with k= (0, 0,0), k= 27a-1(1, 0, 0), k= 2xa-1(4, 0, 0) 
and k= 2na~1(4,4,4) are calculated by the Orthogonalized Plane Wave 
method using the Wigner-Seitz crystal potential and the same electron 
density as in metallic lithium. The curves of energy versus k are drawn 
in the [1, 0, 0] and [1, 1, 1] directions. The energy bands so obtained 
are compared with the energy bands of diamond evaluated by F. HER- 
MAN. The relative positions of the energy bands appear to be the same 
in the two cases and this points to the possibility that the energy-band 
structure is mainly determined by the symmetry of the lattice. A com- 
parison with metallic lithium gives a higher cohesive energy in the case 
of the pseudolattice, when one neglects electrostatic, exchange and cor- 
relation terms. 


1. Introduction. 


The idea of studying the energy bands of pseudo-crystals, i.e. of lattices 
of atoms which are known to crystallize in some other lattice, was first sug- 
gested by F. Serrz as a tool to study the influence of lattice symmetry and 
atomic potential on the energy-band structure. 

CASELLA (1) estimated the energy bands of carbon atoms in the face cen- 
tered cubic lattice by the Orthogonalized Plane Wave method. The results 
he obtained are considerably different from the results HERMAN (?) obtained 
in diamond; in particular at the point J’ (k= 0) the position of the energy- 


(1) R. CASELLA: Phys. Rev., 109, 54 (1958). 
(2) F. HERMAN: Phys. Rev., 88, 1210 (1952); 93, 1214 (1954). 
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states I,;, and Ij, are reversed. À calculation of the energy bands of sodium 
atoms in the diamond structure was performed by F. BASSANI (3) by the 
O.P.W. method and the results displayed a strict similarity with the energy 
bands of silicon (*). 

A calculation of the energy-band structure of lithium atoms in a diamond 
lattice having the same atomic volume as metallic lithium is presented in this 
paper. The O.P.W. method is adopted because of the advantage of evaluating 
several energy bands and because of its successful application both to metals. 
and to insulators. 

In Section 2 an outline of the method is given in a form particularly suitable 
for the diamond symmetry. The energy states of the empty lattice are related 
to the irreducible representations of the group of the wave vector to which 
they belong. The crystal symmetry combinations of plane waves are also 
given at one of the points of interest. 

The parameters which enter the actual caleulations are discussed in 
Section 3 and numerical results for them are presented. 

In Section 4 the results for the energy states of the lowest bands are given 
at four symmetry points of the reduced zone. A comparison is made with 
the energy bands in diamond and in metallic lithium in Section 5. 


9. General method and symmetry analysis. 


The method of orthogonalized plane waves was introduced by HERRING (5) 
and was recently discussed by many authors (5). It is actually à variational 
procedure in which linear variation trial functions for valence and conduction 
states are formed with plane waves made orthogonal to the Bloch functions 
constructed from the atomic core states. For a certain value of the wave 
vector k, inside the reduced zone, the plane waves to be considered are of 
the type exp[i(k+h,)-r], where h; is any vector of the reciprocal lattice. 
These plane waves are the eigenfunetions of the «empty lattice », 16.101 4 
lattice with a constant potential, as classified in the reduced zone scheme. 

From the set of plane waves belonging to a given value of |(k+h,)|*, 
linear combinations are constructed that transform according to a given co- 
lumn of a certain irreducible representation of the group of the k vector. 
In this way the number of undetermined coefficients in the trial function is 


(3) F. Bassani: Journ. of Chem. and Phys. of Solids (to appear). 

(4) T. 0. WooDRUFF: Phys. Rev., 103, 1159 (1956); F. BASSANI: Phys. Rev., 108, 
263 (1957). i 

(5) C. HERRING: Phys. Rev., 57, 1169 (1940). 

(6) J. CALLAWAY: Phys. Rev., 97, 933 (1955); T. O. Wooprurr: Solid State Physics, 
edited by F. Serrz and D. TURNBULL, vol. 4, 367 (1957). 


806 


ENERGY-BAND STRUCTURE OF LITHIUM ATOMS IN THE DIAMOND LATTICE 807 


greatly reduced and convergence is obtained to energy states that are auto- 
matically classified. 

We shall estimate the energy levels at four points of the reduced zone: 
the point I" with k—0, the point X with k= (2x/a) (1, 0, 0), the point A(4), 
with k= (2x/a) (4, 0, 0) and the point L with k = (2x/Ja) (4, 4, 4). For each of 
these points the lowest eigenvalues of the « empty lattice » are listed in Table I 
with the number of plane waves and the irreducible representations of the 
group of the wave vector belonging to each eigenstate. The symbols used are 
those introduced by BOUCKAERT, SMOLUCHOWSKY and WIGNER (7) and by 
HERRING (8). 


TABLE I. — Olassification of the «empty lattice» eigenvalues. 


Number Empty lattice 
of eigenvalues in units Irreducible 
plane waves (72/2m)(47?/a?) representations 
1 (OMO0)= IPE 
Point 8 (il, TIE TP AP yes Me IP 
1h =) 6 (2, 0, 0)? Los (Lo) Por 
12 (2, 2, 0)? r, 25115 Ip T's; 
Point X 2 (1, 0, 0)? À 
k = 2a/a 4 (Oe 151)? xX, X, 
(Us (Ws Ce) 8 (i, 2 Oe AG, AG AG BG, 
8 (250 IE DICA I 
1 (4, 0, 0)? Al 
4 Gat 21) AREA 
Point A(4) 1 (3, 0, 0)? Albe 
k= 2aja 4 (4, 2, 0)? “Alyy, As Abe 
(4, 0, 0) 4 (eee lew lz Ab Aly Als 
4 (3, 2, 0)? A; A, A; 
1 (Cs OO, We Aly, 
2 $$ nu 
Point L 6 (3, 4, 4)? Ie Ab Ie Mire, 
k = 2a/a 6 G, 8, 2? Ty Tigh Lig Uy 
4,3. D 6 (&, 4, be Ty Ly Ls Ly 
2 (8, 3, 5° Ly Dy 


(7) L. P. BoucxaAERT, R. SmoLucHowski and E. WiGNER: Phys. Rev. 50, 58 
(1936). 
(8) C. HERRING: Journ. Franklin Inst., 233, 525 (1942); see also, W. DORING 


and V. ZEHLER: Ann. der Phys., vi, 13, 214 (1953). 
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Linear combinations of plane waves (L.C.P.W.) belonging to à given row 
of any of the irreducible representations are derived from the sets of plane 
waves indicated in Table I by the standard method of projection operators (*). 
Typical crystal symmetry L.C.P.W. at the points I' and X are given in 
Table II of HERMAN®S paper (2). At the point L the C.S.L.C.P.W. are the 
same as in the face centered cubic lattice; a few of them are reported in 
Table V of CASELLA’s paper (1). The C.S.L.C.P.W. we have used at the 
point A(4) are given in Table II. 


TapLe II. — Orystal symmetry linear combinations of plane waves at the point 
k= (2x/a) (3, 0, 0). 
(2a/a)-(k+h) A, As Ar: A. Ae da 
Me 100) a1 — — = à lei 
St +1 A = +- Real = 
Œ 1, 1) ai n a = tal E 
(3; us 1) +1 — ] — — — ] — ] 
(3, ils 11) + 1 — ] — = +1 +1 
| - 4 
EN) “i “i = 7 = ie 
( 1, 1) +1 +1 — — ni i 
(I ACRI = 2a = = +1 =] 
| (3, 1, 1) — ] + 1 — = — 1 — 1 
22.0) | = +1 +1 - 0 —1 
CASO, = == 1 1 ta 0 x 
| (40,2 — +1 —1 - +1 0 
| 0,2) = "sl sl a seal 0 
(8, 0, 0) — +1 - — — — 
(3, 2, 0) +1 — — +1 +1 0 
(E, 2, 0) a] = = =I Les] 0 
E, 0, 2) "al = = 1 0 zl 
(&, 0, 2) —1 _ _ +1 0 — 1 
In the first column the P.W.’s are indicated by the components of the vector (2a/a)~! (k +h). 
The linear combination from a set of plane waves with given value of |k+h| can be obtained 
from the table above multiplying each plane wave by the coefficient listed on the same row 
and in the column of the irreducible representation of interest. The two basis functions of the 
two-dimensional irreducible representation 4; are indicated by dE and AO 


(9) L. D. Lanpav and E. M. Lirsnirz: Quantum Mechanics (translated) (London,. 


1958), p. 336; G. F. Kosrer: Notes on Group Theory, in Technical Report No. 8-M.I.T 
(1956). ty 
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The problem of finding the energy eigenvalues Æ corresponding to every 
irreducible representation thus reduces to that of diagonalizing a secular deter- 
minant whose elements are sums of matrix elements of the operator (H — E) 
between the plane waves of two C.S.L.C.P.W. orthogonalized to the core 
states. The order of the secular determinant is fixed by the number of 
C.S.L.C.P.W. used. If the crystal one-electron potential is chosen to be a 
sum of suitable atomic potentials V, centered at the atomic sites, the matrix 
element of (H — E) between the O.P.W.’s having wave vectors k-+h,; and 
k-+h, for a crystal having diamond symmetry and only (1s) core states has 
the form: 


(1) (k+h,|H—E|k-+ hj) = [(k+ h;)?— E]ò;+ 
+ cos [(h;—h,)-d] [V(|h;—h,|) — (Bi—E) A; (| + hi |) Ar.(| + Fy). 


The vector d has components a/4 (1,1,1), a being the lattice constant; 
V(|h,— h;|) is the Fourier coefficient of the potential defined by: 


= 
bo 
>, 


V(|h.— h;|)) = 2 (V.0) exp[—i(h, —h,)-r]d°r, 


where ©, is the volume per atom; H,, is the energy of the core state (1s); 
finally the A,, are the orthogonality coefficients defined by: 


(3) Axs(lk + h,|)= 9 fre exp[i(k + h,)-r]d?r, 


where g,, is the atomic wave function. It is assumed that the Bloch sums 
for the core states are eigenfunctions of the same hamiltonian H as the valence 
and conduction states. | 


3. — Calculation of the parameters. 


From Section 2 it can be seen that the parameters needed for the appli- 
cation of the O.P.W. method are: the lattice constant a, the core state energy 
E,,, the orthogonality coefficients given by (3) and the Fourier coefficients 
of the potential given by (2). For the present calculation it is assumed that 
the atomic volume Q, of the hypothetical lithium in the diamond structure 
is the same as the atomic volume of actual body-centered cubic lithium. This 
seemed a reasonable choice also because of the experimental evidence that 
allotropic transformations cause only slight density variations. The lattice 
constant of metallic lithium extrapolated at 0 °K is taken as 6.5183 in atomic 
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units; hence the atomic volume is ©,— 138.47 a.u. and the lattice constant 
in the diamond lattice turns out to be a= 10.347 a.u. For the choice of the 
potential we have followed a recent calculation of energy levels in metallic 
lithium performed by GLASSER and CALLAWAY (1). They have used the semi- 
empirical free atom potential constructed by SEITZ (11) and have found by 
numerical integration a core energy H,,=— 3.765 Ry and its appropriate 
wave function. We have taken the same value for £,, and have approximated 
the wave function by the sum of two exponentials q,,(r)= A exp [— ar] + 
+B exp[— fr]; the parameters A, B, a, B were adjusted so as to obtain all 
the orthogonality coefficients given by GLASSER and CALLAWAY to within 1%. 
The error so introduced is probably less than the error involved in using atomic 
functions in the Bloch sum, neglecting overlap. The satisfactory values of 
the parameters were found to be A = 6.422; B= 2.171; «= 3.019; B — 1.907; 
it was verified that with these values the approximated function y,, was auto- 
matically normalized to unity within 1%. The orthogonality coefficients ob- 
tained from formula (3) by using this analytical approximation for @,,(r) are 
listed in Table III. 


TABLE II. Orthogonality coefficients between the 1s core wave function and the plane waves. 


(a?/420°)(k +h)? A,,(|k +h |) (a?/4r°\(k+h)® | A,,(\|k+h)) 
0 0.32901 4 0.19955 
i 0.31701 w 0.194 41 
3 0.29517 2 0.18478 
1 0.285 21 5 0.18026 
2 0.25059 6 0.16390 
; 0.243 05 25 0.16019 
a 0.229 09 a 0.15318 
3 0.22261 8 0.13771 


The crystal potential was constructed following the Wigner-Seitz method 
as in metallic lithium. In the diamond structure the atomic polyhedra are 
truncated tetrahedra rather different from spheres. In order to obtain the 
Fourier coefficients of the potential, however, the integral (2) has been extended 
over the atomic sphere of radius r, — (3.Q,/42)* = 3.209 4 a.u. Overlapping of 
nearest spheres corrects in the right way the atomic spherically symmetric 


(9) M. L. GLASSER and J. CALLAWAY: Phys. Rev., 109, 1541 (1958). 


(4) F. Serrz: Phys. Rev., 47, 400 (1935). The potential for Lithium is correctly 
tabulated by W. KoHN and N. RoSTOkER: Phys. Rev., 94, 1111 (1954) 
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potential about each nucleus. The Fourier coefficients calculated from eq. (2) 
by numerical integration of the Seitz potential are listed in Table IV. 


potential (in Rydbergs). 


TABLE IV. — Fourier coefficients of the erystal 


4. — Results. 


(a?/42)|h|? V(Ih}) (a?/4r?)| h|?2 V(Ih|) 
0 — 1.0022 19 — 0.07277 
3 — 0.38012 24 — 0.06856 
8 — 0.06216 oF — 0.05799 
11 — 0.04070 32 = 002813 
16 — 0.063 07 


The values given in Section 3 for the various parameters can be substi- 
tuted in the matrix elements of the secular determinant described in Section 2. 
The various roots of each of the resulting secular equations converge to the 
successive eigenvalues of the Hamiltonian as has been proven by MACDONALD (1). 
The determinantal equations of third, fourth and fifth order have been solved 
with the electronic computer of the Engineering School in Milan. The results 


obtained for the lowest energy 
levels at different stages of ap- 
proximation are given in Table V. 

It can be seen from Table V 
that the convergence is not very 
good with the number of plane 
waves we have included. This 
applies particularly to the eigen- 
states are automatically 
orthogonal to the s states by rea- 
sons of symmetry, as it was found 
by F. HERMAN in the case of dia- 
mond. Nevertheless we feel con- 
fident that the present approx- 
imation is sufficiently accurate 
to determine the energy-band 
structure in the reduced zone. 
Adding plane waves of higher 
energy would lower all levels with- 


which 


L, 

à 
22111 =0 27100) K- 28000) 
Kao Se) We Tag Yer a 
Fig. 1. — Profile of some energy bands of 


lithium atoms in the diamond structure in 
the [1, 0, 0] and [1, 1, 1] directions. 


(12) J. K. L. MACDoNnALD: Phys. Rev., 43, 830 (1933). 
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out probably changing their relative positions. This was shown to be the 
case for the J” point in diamond (2), silicon (*) and germanium (1°). 

The energy values given in Table V, the compatibility relations and the 
non-crossing rule allow to draw the profile of the lowest energy bands along the 
symmetry directions [1, 0, 0] and [1, 1, 1], which is presented in Kies Ae Nae 
turally this procedure leaves uncertain some of the detailed features of the 
energy bands such as the curvatures, but we have preferred not to use an 
interpolation scheme, since we are mainly interested in the sequence of the 


energy bands and in their qualitative features. 


TABLE V. — Energies at the points k= 0 k= 2ra-* (1, 0, 0), k= 2a (4, 0, 0), 
k= 2a (4,4,4,) (in Rydbergs). 


(13) F. HERMAN: Physica, 20, 801 (1954). 


le IB ys, ID Wy Fe ER 
Ei — 0.667 0.166 0.166 0.828 — = 
E, — 0.757 (II 0.068 0.825 — — 
Es =O 771! 0206 = == = 
X,(1) XxX, X,(2) Xx, = = 
Ki 20950 — 0.203 — 0.904 — — 
E, — 0.439 0922 0.208 04732 _ _ 
IDA — 0.455 — 0.392 0.186 = = ssa 
E, — 0.503 — 0.121 — — — 
E; = 05534 — 0.111 — _ — 
A,(1) A;(1) A,,(1) A,(2) A,,(2) A;(2) 
E, — 0.591 = 0110 0.053 — = = 
E, — 0,647 — 0.243 — 0.127 0.253 0.375 0.627 
Ey — 0.688 — 0.262 — 0.142 0.138 0.316 0.385 
E, — 0.701 — 0.320 — 0.237 0.099 0.233 0.376 
Es => = — 0.250 = = Pe 
L,(1) L(1) Ls: Ls L,(2) Le) 
E, — 0.370 — 0.510 — 0.166 0.314 — — 
E, = 0,427 — 0.551 — 0.223 0.220 0.477 0.668 
Es — 0.565 — 0.553 — 0.268 0.204 0.476 0.666 
E, — 0.574 — 0.569 — i 0.385 0,662 
B; — 0.598 — 0.589 = — 0.327 0.577 
The subscripts attached to Æ indicate the order of the secular determinant used in ob- 
taining the energies. 
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5. — Conclusions. 


First it seems of interest to compare Fig. 1 with Fig. 1 of HERMAN’s paper 
on the diamond energy-band calculation with the O.P.W. method. There the 
crystal potential is taken as a sum of carbon atomic potentials extending to 
infinity, furthermore the carbon atoms are more closely packed than the lithium 
atoms in the pseudo-crystal, the lattice constant in diamond being 6.73 a.u. 
to be compared with our 10.347 a.u. In spite of these differences the sequence 
of the energy bands at each point is the same in both cases and their relative 
positions at different points of the reduced zone are also preserved. The max- 
imum of the valence bands and the minimum of the conduction bands occur 
at the center of the zone, the energy gap being larger in the case of diamond 
(0.45 Ry to 0.27 Ry). This last result is due to the smaller electron density 
in the pseudo-crystal than in diamond, which results in a smaller width of 
the energy bands consistently with the fact that in the «empty lattice » the 
eigenvalues are closer when the lattice constant is larger. A comparison with 
Herman’s data also shows that in our case the states J and 7”, are shifted 
to higher energies with respect to /\; and J.,;,. This can be attributed to the 
larger effect of the (1s) core state, whose wave function is more spread out in 
lithium than in carbon. 

These results seem to point to the possibility discussed by T.O. WOODRUFF (°) 
and by J. CALLAWAY (4) that the sequence of the energy bands is independent 
of the specific atomic potential and can be predicted by the lattice symmetry 
and the atomic structure. 

A calculation analogous to the one reported here was recently performed 
by F. BASSANI (3) on sodium atoms in the diamond lattice. The energy bands 
which result in that case display a strict similarity with the energy bands of 
silicon. They both differ from the bands in diamond and in our pseudo-crystal 
because the energy at J; lies higher than at X,(2). This may be attributed 
to the presence of the (2p) core states in sodium and silicon, since they raise 
the energies at 7°; and /,,, more than at X,(2). In fact the leading crystal 
symmetry combination of plane waves which corresponds to the eigenstate 
X,(2) is orthogonal to the p states and consequently, in the perturbation ap- 
proach described in Appendix A, the eigenvalue is insensitive to the effect 
of the p core states, up to the second order. 

It might be of interest to compare our results with the energy bands in 
metallic lithium (1%), In that case there is no possibility for an energy gap 
due to the order of appearance of the energies at the points k=O and 


(14) J. CALLAWAY: Phys. Rev., 103, 1219 (1956). 
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k—27at (1, 0, 0) in the b.c.c. lattice. The value of the energy at /7 lies lower 
in our case than in lithium. The value of the effective mass obtained averaging 
the values calculated by fitting parabolas to the points given in Fig. 1 for 
the lower band is about the same as in metallic lithium (53 

For these reasons, if one neglects the contribution to the cohesive energy 
due to electrostatic interactions, exchange and correlation effects, the cohesive 
energy should be greater in our case than in metallic lithium by about 
0.1 Ry. Very likely this means that in the diamond structure electrostatic 
interaction, exchange and correlation effects do not compensate as in the 
body centered cubic structure (1). 


* OK * 


The authors are indebted to Professor F. G. FumI for his interest and 
advices; they wish to thank also Dr. D. LYNCH for a reading of the manuscript. 


APPENDIX 


The determination of the crystal eigenstates by the O.P.W. method may 
be regarded approximately as a perturbation procedure on the «empty lattice » 
eigenstates. 

A crystal eigenfunction y%(k, r) belonging to the « irreducible represen- 
tation of the group of k can be expanded as follows: 


(A.1) vi(k, r) = 8.(k, r) — > <@,|8,> Ok, r) + 
+ > a;S,(k, r) — Y D a;<®,|S,;>@,(k, r) +... 


Here S,(k,r) are the linear combinations of plane waves belonging to the « 
irreducible representation; the unperturbed wave function S,(k, r) belongs to 
an eigenstate of the « empty lattice », whose degeneracy has been fully removed 
by the symmetry analysis. The Bloch wave function of the / core state has 
been indicated by @,(k,r). The expansion coefficients a; and the Square 
moduli of the « orthogonality coefficients » <®,|8;> will be considered small 
quantities, higher order terms having been omitted in expansion (A.1). 


The crystal Hamiltonian H, in atomic units, is: 
(A.2) ETNIE IO Va 


In this formula V(0) is the Fourier coefficient of the crystal potential V for 
h=0, i.e. the space average of V, and V, = V— V(0) shall be considered as 


(5) R. H. PARMENTER: Phys. Rev., 86, 552 (1952). 
(6) F. Serrz: Modern Theory of Solids (New York, 1940), p. 364. 
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a perturbing potential. Let the eigenvalue of H, belonging to the eigen- 
function S; be W,; let the first and second order corrections to the unperturbed 
energy W, be W; and W;. By substituting expansion (A.1) into the equation: 


(A.3) Hyi(k,r)=(W+Wi+ Wi+...)yi(k, r) , 
and using the techniques of perturbation theory one obtains to first order: 


(A4) <S;|V,|S5 — DD, |S;> <S;|H|@> + a;W, 
L 
SS > <B,|8,> (8; |D) +a,W;. 
I 


As usual, one can assume for the core states: 
(A.5) HD(k,r)= ED{(k,r), 
where £, is the atomic value. Then from (A.4): 


(A.6) Wi = 8. |Vil|So + > (W:— E)|kK(D.| Sp]. 
(A 


The values of the coefficients a; can be obtained from formula (A.4) for è # j. 
The value of the coefficient a; is found to be + 3 DD | |S;>| by requiring nor- 
malization of w;(k, r) to first order. 

Proceeding to second order one finds: 


(A.7) Wi = > (W; EE = W,)- | <8; | VANO, — (E, e Wi) {D, |S;> {S;| Py fe 
i) 1 
+ X|KD:]S) Wi. 


It can be seen from formulae (A.6) and (A.7) that the core wave function 
@,(k, rr) has no effect on the energy when <@,|S,> = 0. 


RIASSUNTO 


Si considera un ipotetico cristallo costituito di atomi di litio nella struttura del 
diamante e avente la stessa densità di elettroni del litio metallico. Usando un potenziale 
del tipo di Wigner e Seitz, si calcolano col metodo delle onde piane ortogonalizzate 
i valori delle energie nei punti speciali della zona ridotta con K= (0, 0, 0), k—2xa-t (1, 0,0), 
k=2aa-1(k, 0,0) e kK=27a-(}, 3, 3) e si tracciano le curve dell’energia in funzione 
di k nelle direzioni [1, 0, 0] e [1, 1,1]. Le bande di energia così ottenute sono para- 
gonate con le bande di energia del diamante valutate da F. HERMAN. Il fatto che le 
posizioni relative delle bande risultino le stesse nei due casi indica la possibilità che 
la struttura delle bande di energia sia principalmente determinata dalla simmetria 
del reticolo. Dal confronto con il litio metallico, si trova che l'energia coesiva è mag- 
giore nel caso dello pseudoreticolo, quando si trascurino i termini elettrostatici, di 
scambio e di correlazione. 
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Quadratic Lagrangians and General Relativity. 


LAR\VARHIGGS 
Department of Mathematics, Imperial College - London (*) 


(ricevuto il 31 Dicembre 1958) 


Summary. — Two out of the three sets of field equations associated with 
Lagrangians quadratic in the Riemann-Christoffel tensor, which were 
recently discussed by Stephenson, are shown to be equivalent to the 
source-free Einstein equations with an arbitrary « cosmological » term. 


1. — Introduction. 


In à recent paper STEPHENSON (1) has discussed the field equations which 
one obtains by applying the Palatini method (independent variation of the 
metric tensor g,, and the affine connection 1,5) to action integrals constructed 
from three Lagrangian densities quadratic in the Riemann-Christoffel tensor, 
it being assumed that both metric and connection are symmetric. The fol- 
lowing are the action integrals and associated field equations. 


Case (A). 12 = fata V— y; 
(1) R(Ru— ERg,,) = 0, 
(2) (Rg*? /— 9), = 0. 
Case (B). I; = [star heyy, 
(3) RR lg + odo Sally 
(4) VD, =0. 


(*) Now at the Department of Mathematics, University College, London. 
(1) G. STEPHENSON: Nuovo Cimento, 9, 263 (1958). 
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Case (C). Io = fata Re mae 
o) i Te ET + or Bey at 7 ou > Lai tee Om =0, 
apy 
(6) (RB, V—9)y=0. 


‘ 


In these equations 


# ss x # x œ x œ 
R Ai te Ay nig Li Avy de uP li An — te li Au ? 
— v 
R,, = È Aur 


and the remaining tensors are constructed from these with the aid of g,, in 
the usual way. A comma indicates ordinary differentiation, a semicolon co- 
variant differentiation with respect to /”,,. A bar placed under tensor indices 
denotes symmetrization. 

STEPHENSON (2) has further shown that all these Lagrangian densities are 


invariant under the Weyl gauge transformations 


Guy Fs PI ny ? 


x # 
LR", 


(7) 


where g(x) is an arbitrary scalar function of the co-ordinates. A direct con- 
sequence of this invariance is the vanishing of the trace of one set of equa- 
tions (namely (1), (3), (5)) in each case. 

The purpose of this note is to show that in cases (A) and (B) one set of 
equations may be transformed into field equations of the Einstein type with 
an arbitrary « cosmological constant » in terms of a new gauge-invariant metric. 
The otber set provides algebraic relations between the old metric and the new, 
the solution of which contains the arbitrary gauge function. No such trans- 
formation has yet been found in case (C). 


2. — Case (A). 


Assuming that R 0, we first take the trivial step of rewriting (1), (2) 
in the form 


(1a) Ruy = Ay" Inv ’ 


(2a) (prig'V_g),=0, 


(2) G. STEPHENSON: Quart. Journ. Math. (Oxford), to be published. 
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where y(x) is an arbitrary dimensionless scalar function and / is an arbitrary 
constant of dimension (length)-2. If we now introduce a new metric tensor 4,,, 


defined by 
(8) I uv <= Yaw ? 


together with its contravariant inverse a” and its determinant a, these equa- 


tions become 


(9) Rw ca My ? 


(10) (a? V/— a), = 0. 


Clearly we may regard a,, as a gauge scalar and write the gauge transfor- 


lod 


mation (7) as y — py. 
The most general solution of (10) is 


x pei il xO = 
(11) uf: GA 24 (Cr = Wo Bo. 46,0) ; 
Consequently È, is symmetric and equations (9) are of the Einstein type with 
an arbitrary cosmological constant 2. They are to be regarded as second-order 
differential equations in the gauge invariant metric a,,. From any solution 


of these equations we may recover a family of metrics In» by using equation (8), 
which contains the arbitrary gauge function y. 


3. — Case (B). 
This time we use as our definition of the new metric b,, the equations 
(12) DEV —D = DIR gg —g , 


where again À is an arbitrary constant of dimension (length)-2. Equations (4) 
then become 


(13) (b%-/— 5), = 0, 

of which the general solution is 

(14) 182 ca 30°°(b,.0.8 Da De bp.) 2 

As a consequence, R,, is again symmetric and we may rewrite (12) as 


(15) Buy = AV BG Guy Iyp 0” . 
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On substituting this expression into (3) we obtain 


Lo Ce a 
PRO po 0° Oey i= Poe I pod TRE 


which may be rewritten as 


‘ 


(16) Grab gb = pd e; 


where (x) is an arbitrary scalar function. By combining (15) and (16) we 
finally obtain once more 


(17) Re =D 


which is to be regarded as a second-order differential equation in b,,. From 
any solution of (17) a family of metrics g,, may be constructed by solving the 
algebraic equations (16). These have amongst their solutions the simple re- 
lation 


(18) Iu» = yo, : 
However, (16) may be written in matrix notation as 
(gb!) = y?1, 
from which it is clear that besides (18) there are other classes of solutions 
corresponding to various square roots of the 4X4 unit matrix. In all the 


is to be regarded as a gauge scalar and y as an arbitrary 
gauge function transforming under (7) as y qu. 


elasses of solutions b, 


uv 


8. — Case (C). 


In both cases (A) and (B) we have been enabled to define the new metric 
by the fact that the second set of equations, (2), (4), has the form 
da Fa 0 ? 
where T* is a symmetric tensor density. This is not true in case (C), so it 
does not appear that equations (5) and (6) can be transformed into the 
Einstein equations. 


a 
Lon] 
ao 
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4. — Remarks. 


One comment should perhaps be made on the interpretation of the metric 
tensors in these gauge invariant theories. Presumably one wants the concept 
of length in such theories to be gauge invariant. Therefore it seems likely 
that more direct physical significance may be attached to a,, and b,, than to 
the original dynamic variables g,, which enter into the action principle. This 
is analogous to the situation in electrodynamics. 


* * % 


I am grateful to Dr. G. STEPHENSON for arousing my interest in this problem 
and for various helpful comments. The work was performed during the tenure 
of an I.C.I. Research Fellowship. 


RIASSUNTO () 


Si dimostra che due dei tre sistemi di equazioni di campo associati a lagrangiani 
quadratici nel tensore di Riemann-Christoffel, recentemente discussi da STEPHENSON, 
sono equivalenti alle equazioni di Einstein prive di sorgente con un termine « cosmo- 
logico » arbitrario. | 


(*) Traduzione a cura della Redazione, 
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Spectral Functions in the Static Theory. 


J. D. WALECKA (*) 
CERN - Geneva 


(ricevuto il 31 Dicembre 1958) 


Summary. — The spectral functions of the electromagnetic form factors 
of the nucleon are discussed within the framework of the static theory. 
1/M recoil corrections are also included and it is shown that these cor- 
rections are close to those given by dispersion theory and that they pro- 
foundly alter the spectral functions indicating that the application of the 
static model to this problem is of doubtful validity. 


1. — Introduction. 


The spectral functions G, ,,(0°) describing the electromagnetic properties of 
nucleons are related to the form factors F’,,,(q2) by 


œ 


n 9903) 
(1) ea a) =i af ’ 


au? 


where qg? is the four-momentum transfer occurring at the electromagnetic vertex 
and u= m_. The spectral functions can be regarded as the quantities of fun- 
damental theoretical interest and there have been two recent attempts at the 
calculation of these quantities using dispersion theory, one by CHEW, KARPLUS, 
GASIOROWICZ, and ZACHARIASEN and the other by FEDERBUSH, GOLDBERGER, 
and TREIMAN (1°). 


(‘) National Science Foundation Postdoctoral Fellow. 

(1) G. F. Cnew, R. KaRPLUS, 8. Gasiorowicz and F. ZACHIARIASEN: Phys. 
Rev., 110, 265 (1958). 

(2) P. FEDERBUSCH, M. GOLDBERGER and 8. TREIMAN: Phys. Rev., 112, 642 (1958). 
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The purpose of this paper is two-fold. First it will be shown how one can 
calculate the spectral functions within the framework of the static theory and 
that the functions thus calculated agree with the static limit of the dispersion 
theory result. The validity of equation (1) has not yet been proved in the 
framework of relativistic field theory. Moreover, the computations of the 
spectral functions mentioned above are based on the analytic continuation of 
the pion-nucleon scattering amplitude for unphysical values of the angle. 
Therefore, it is thought useful to associate with these calculations a straight- 
forward computation based on a simple model. 

The advantage of dealing directly with the spectral functions instead of 
the form factors in the static theory is that in the no cut-off limit the former 
are finite whereas the latter are not. Now it might be hoped that the static 
theory would at least explain the low o? part of the spectral functions, but 
the work of references (1) and (?) has shown that the discrepancy between 
dispersion theory results and their static limits are great and that they persist 
even down to o? = 4u? (*) Therefore, the second problem dealt with here will 
be the result of including the 1/M nucleon recoil corrections into the frame- 
work of the static model and showing, that 


a) the presence of these purely kinematical corrections profoundly alters 
the spectral function; 


b) the kinematical 1/M corrections are in good agreement with the 1/M 
corrections calculated by dispersion theory techniques indicating that 
this is the main source of difference between the two theories. 


The first conclusion indicates that one can have little faith in the pre- 
dictions of the static theory in the electromagnetic structure problem. Also, 
since the partial waves of higher ? value (the scattering in these states is known 
neither experimentally, nor theoretically), become very important in the cal- 
culation of G, ,(0°), there may be some doubt as to the validity of using 
only S and P waves in the analytic continuation of the meson-nucleon scat- 
tering amplitude as is done in the dispersion theory approach. 

The calculation of the electromagnetic properties of the nucleon in the 
static theory has been discussed by several authors (4). In this work the nota- 
tion of the last author (F) shall be used. The introduction of nucleon recoil 
is based on the work of Bosco, FUBINI and STANGHELLINI (5), who use the 


(3) S. DRELL: Proceedings of the 8-th Annual Conference on High Energy Physics 
(Geneva, 1958). 

(*) G. SALZMAN: Phys. Rev., 99, 973 (1955); 105, 1076 (1957); H. Miyazawa: Phys. 
Rev., 101, 1564 (1956); S. TREIMAN and R. G. Sacus: Phys. Rev., 103, 435 (1956); 
F. ZACHIARIASEN: Phys. Rev., 102, 295 (1956); S. FUBINI: Nuovo Cimento, 3, 1425 (1956). 

(5) B. Bosco, S. FUBINI and A. STANGHELLINI: to be published. 
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procedure to reproduce, up through D waves, the meson-nucleon scattering 
phase shifts calculated from dispersion theory. 

In Section 2, the working equations are derived. In Section 8 the spectral 
functions are calculated and in Section 4 they are related to the corresponding 
quantities from the relativistic theory. 


2. — Derivation of equations. 


In the following, the Breit frame will be used, that is, the initial nucleon 
momentum is — q/2 and the final nucleon momentum is +q/2, where q is 
the momentum of the photon from the external field. In this system there 
is no energy transferred to the nucleon and hence q? = q? and also there is 
only momentum vector, q for constructing invariants. Now 


A|+q/2> = Baus| + q/2> = (M? + (q/2)*)*| + 4/2) , 


where H is the total meson-nucleon Hamiltonian and | + g/2) is a physical 
nucleon. One thus includes the proper kinematics by replacing H by E,,, when 
acting on a physical nucleon state but then uses the static model for the 
dynamics and sums overs states. The reason for this is that the dynamics of 
the static model are relatively simple. To calculate corrections rigorously to 
1/M one must also, of course, include 1/M modifications of the dynamics or H, 
but this complicates matters and it will be shown that the kinematical cor- 
rections alone are sufficient to give 1/M corrections similar to those of disper- 
sion theory. 
One starts from (Ff) 


ie exp [i(k + kai 


oe 20 (cop)? [(@e — x) (Agg dgr — À pra DE 11) — 
And k,k' xp )® 
— (Wz + Wy) (AT pra ge — A ra Aer)] » 
; — ie x exp[i(k + k)-a] : 
j(æ) PATTES 20 si (0x0p HE a CA, k,}(4% we (07251 ae 


K ye kK 
+ at pa pre + pe LE er + 02411) » 


: + 
which are the expressions for the charge and current densities in the meson 
field. The subscripts 1 and 2 refer to the isotopic spin of the Hermitian charged 
fields. What one wants to calculate is the transition probability 


wo = <q/2 fatelo (0) 0 (0) — j(x): A%(x)) | —q/2>, 


DI 
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where one takes 
@™ (x) = y exp [iq: x] 
A” (@) = a exp lig: x] 
and defines (using invariance properties). 
We f,(q°)p UE (g/\a)fy(q°) È 


The fs are the form factors in the static theory. The procedure is to relate 
the expectation value of pairs of creation and annihilation operators to the 
meson nucleon scattering amplitudes (F). 47,4, Shall be used as an example 
in the following and then the results given for all pairs, starting from 


H=H,+H,, 
[Ho, Ogg] = Opty > 
[Fa Gel = Vis 
{q[2|[H, aka da ]|— q/2> = 0, 
one finds 


ES co 


Op — On % 


nlare|— g/2> — <q/2 | ak |n> <n | Vitw |— 9/2) - 


q/2 | ia Ga’ |—q/2> = 


Taking matrix elements of the following relations (p is the total momentum 
operator) 


LP; Ai, = kay, ’ 
Lp, Dig Oy | =. (k DE CE et , 
one gets 
(Pn Pa — k)<m | ax, |n> = 0 ) 
(Pm— Pn + k)<m|a,,.|n> = 0 ’ 


(Pn — pr k + ')<m| ap, aya |) = 0 


expressing the conservation of momentum. 
Also from 


[H, ax, | | a q/2> = (H — E 12) Uy | + q/2> 
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one has 
Cy, | ab q/2) = eee Lies Vi | le q/2) 
Sad PME 7 x | Ay 
H+ — Egg © 
Therefore 
<q? lat, aug | — le) | = 2 y <al2 Pr in) Gil Vibe | > | 
\ petal | Yard CA) APN SE a 33 E = 
VO eee hoa oat Wa + @y — Lys y 
(q=k-k’) (ky, = —q/2—k’) 
IT, | + ony 
<q/2| y ka | > <n| Via | == q/2> | 
ar, le CE Se BEA | 
| 9 
a RO E are | 


(ky, = q/2—k) 


W, is the energy of the intermediate state n. Using 


W, = (M?+ KR) 
and 


(= Nh, AVE 
one finds.to order 1/M (i.e. if w,/M<1) 


? k- k' 

di ma Hare = y» + oM . 

Until now the derivation has been rigorous. 

One now uses the static model to evaluate the numerators. 

Replacing | + q/2>—|R>, where | R> denotes one of the four physical nucleon 
spin states one has 

; ‘ 
rie 1 tt |v ka | UD, «n | Vis | R> 


cal? |ats axe'|—42) == 
È di ! 4 Ri 4 On Oe wae OO, =F Ox = (k- k')/2M 


CR | Vinw| n> <1 | Via |) 
On ae Op! ae (k- k')/2M 


1 6B | Vra|n) <1 | Vive | BD 


Ox 5 Ox n On = On (k- k')/2M KE 
ER | Vine | n) «n | Vee | È) 


On + Op — (k- k')2M ; 


‘ 


f Guests | ‘= 
<q/2 | diate’ | — q/2> = 


= 


(42 | Ghee axy |— 9/2) = - LL CR) Vision] Va E) 
| | x % i O% + Ox n On == Om (k- k')/2M | 


(B| Vity In) <n | Vie LR) 
On+ @y —(k-k')2M ? 


| 


a 53 - Il Nuovo Cimento. 
fos) 
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fe | 
Fa i(4) (sax) (+ h)r,0(h) 


Notice that the only difference from (F) is the extra term in the denominators. 
Proceeding as in (F), one writes the unitarity principle 


where 


VERRE Ge aes 
pvp) Kia; r| DT! TT \ jo; ry =(F Ab: D b(n — 0,)(R|o:t, 


n —)(n—|o;ty|R) , 


where 7 is defined between free meson and bare spinor states and n — is an 
incoming scattering state. Multiplying this equation by 1/(@,+@,+(k:k'/2M)) 
integrating on ©, from w to co and noting that the sum one needs is a sum 
over all states and not just scattering states one finds (*) 


(AE (R|o.t,.|n)(n|0;ty |R) 


u varese] k-k)2M 


Tor (di (leo: dest) do, Im T(@,) | 
Re at oi uf ox+(k-k)2M p?v2(p) ox + @,+ (k-k')/2M 


CETTE 


Vad 


where f is now the renormalized coupling constant. One now takes into ac- 
count the experimental evidence of the dominance of the 3-3 scattering reso- 
nance. Approximating the contribution of the resonance by a 6-function 
one has 


(cel 


il do, Im T(,) _ fa\? (2 + €-4)(2 + 1-6) 
J v?(p)p° Ope CU k')/2M  \u) w+ 0x4 (k-k')/2M' 


li 


with a? the effective strength (a? = (4/9)f? theoretically) and ©, the resonance 
energy (W, 2u), one finds 


TE 
(g/2|ax,a_xy |— q/2) = > al : =) k,kj<iar | R(w,) — Rwy) \ja'r> , 
ig Or — 
, i ; , C(k, k') i cee STAR 
(q/2 | ak 4% |— q/2) = — 2 ore kik; tar | R(x) + R(0w)|jatr), 
c(h, k') So = Ce 
(q/2 rx Ure! og | DZ q/2) = — - nu k;k;Ciar | R(0x) al R( wx) | joe rv) , 


(*) Note (6): 
Gia; r|1 — bot |jals 1) = 8X7 | tate |) = Cho's r|1+ Et lin; n), 
ius 7 |1— Lo jet: 7) = dgr |0:0; |") = Cjo’s r|1+ L-o]ia; 79 . 
(5) M. Cini and 8. FuUBINI: Nuovo Cimento. 3, 764 (1956). 
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CO 
bo 
SI 


with 


i 
È a (Mx Ox! )E 


were 


R(w) = R;(w) + R.(0)(0-14 tt) + R;(o)(o-l)(E-T) 


Cicer | (aw) |jx'r > = <ja'r| R(w) |iar’d , 


Paap, k- k' a\? il k- k' 
Ro) = |. — A TREES 
n) \ox 2 Mo A) \0xr + Oe 2M (ox + ©)? 


ae kek me “\"( 1 k-k' 
u} \@x OM UL) \Ox + Op | 2M (oz + wr)? 3 


R(x) a | 
Di leo a kik' 
Fis(ox) = ¢ ox 2Moë) | u Wy + We | 2 M (co, + = ? 


R(©,)is obtained from R(©,) by interchanging k and k'. Substituting in the 
expressions for c(x) and j(x) one has 


= ak d3k'v k-k\dk+ k' + 
FAC be) = oa = [TS v(k)v( oe = ( = q) [R.(0,) — Ri(ox)] , 


> 


—io-(q/\a)f,(q?) = E: 


m2 


|| dk d*k'o(k)v(k)(K- a)(o oa: (k/\ k')(0(k+k'+q)[R (x) R,(ow ) 


k? — (Ve OX Ox! 


3. — The spectral function. 


The problem now is to transform f,(g?) and /,,(q?) to a spectral represen- 
tation. To do this it is useful to have the following integral representations 
of the terms appearing in À;(w,). These representations are derived by con- 
sidering the Cauchy integral in the complex k? plane with cuts from iu to 
ico and — iu to — ico and deforming the contours to run along the cuts paying 
attention to the proper sign of the square root along each side of the cuts. 


‘0,(a?) da? 
h;(@x) -( ’ 


~ 
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a h(@x) 0: (a) ’ 
1 1 Vo?— n? 
Gi E À 
n On Op mau + WF 
i 20 p Vo — Ha 
2 ETE LAG SATO 
(0x4 Or)? Tui On)? 
3 oli 1 Il Op 
[ay a n ——— Toe) ; 5 
Ox (Or + Op) a Vo pu? On + a? — y? 
ih al ai OR + u?— oe 
4 — A = = rear 2 = Xe AUD 
Or (Ox + Op)? Te W/o ee ie (02 uw)? 


The contributions from the resonance will be calculated first and the con- 
tributions from the Born terms can easily be found by the substitutions 


Mega: Ds 10 — 2a? > fr, 


Ia (di): On > 0 a? > fr. 


One may note in passing that since the contribution to f,(q?) from the 
Born term and resonance enter as f? and — 2a?, and since a? = (4/9)f?, this 
term will be quite sensitive to assumptions about the resonance whereas f,,(q?) 
will not be. 

One now proceeds as follows: 


a) Notice first that 


LI EN il DI — 1 
ke k'°|a® +k? +R] (02 + 2) (a+ he!’ 


b) It is convenient to forget about the cut-off at this time. It will be 
needed further on to justify a partial integration but the final spectral functions, 
since they give the low mass behaviour, do not depend on the cut-off. 


e) Introduce a change of variables 


d°kd3k'= d'u dv. 


This gives a ò(u — q) and the w integration can be performed. 


828 


SPECTRAL FUNCTIONS IN THE STATIC THEORY 


co 
bo 
© 


d) Introduce the Feynman parametrization 


1 


il il 
af [aa + b(1— 2)?’ 


0 


and make a second change of variables to t= v — (q/2)(1 — 2x). The deno- 
minators become independent of angle and the angular integrations in the 
numerators can be carried out. Note that 


Jac, (terms odd in t)=0, 
dor 


fron = alias b). 


The results of these four steps are 


SA GELO ds Le) 
ian, ja fe Jre + ae 


oy On a ee 


: fat TI, 
ACRI he = (=) fu ‘fv 425 TAN ! 203(02) — ON t2(t2 — q?x(1 — 2))0,(2?) 


wheres 4 —='q'o(l ==) a 


e) The q? dependence must now be removed from the numerators. To 
do this use 


When this is done, one finds one term in f7(q?) which has no q? dependence. 
Thus one must write 


f2(q*) = (00) + fa(q? 


and including the Born term one has 


lr) = 


nm 
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where an upper limit À has been put in the ¢ integration. Since this term cor- 
responds to f,(q°? = co), one can say nothing about it on the basis of the static 
theory. It can be removed by relaxing the condition that equation (1) also 
gives the correct charge and performing one subtraction (!) 


aye rag 2 d 2 
Fa) #,(0) = — | ARE 


MIA) 
ap 
Therefore it is the function FE (q?) which shall be investigated further. 
The remaining terms in 1/[&+4?| are removed by a partial integration 
setting 


er _tdt pate i à s =[what’s left]. 
[EP AP 


dr 


Note two things: 
1) Since there is a cut-off sr NM = 


2) Since one is interested in the spectral function for low mass behaviour, 
the derivative of the cut-off funetion may be set equal to zero. This results in 


a co (co) 
2etz (a\? fi di 1 | 
PES 3 i. ; ) ‘ a 
fu(q?) = ra (5) Joo] def [se 03(a?) — OM (8t4 + 32a?) o,(a?)} , 
0 (Tia 0 
1 (ce) 
n on —4e:,;fa\? di 
e(Q°) = se (5) faefar PL Aa 
0 tt 0 
Pe (ee ee Cr 
Jar a OL") 05 at) aM ee ain) OCS ance? Nie ae) CAI) 
gq) Consider 
+ A? = + «4+ q'a(l — 0) = 
4 ao? 
== Gode Se ea 
| È da æ(1 — i 
One can write 
i n do°d(0° — (t+ a2)/e(1— a)) 
q? + (+ «?)/a(1— x) JC = 


4 
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where the lower limit comes from the fact that the minimum value of 


ORS Sico 
1? + 2 
— for 2 <a? < 00 is 4u?. 
r(1 — 2) a a 4 : & 
<p Al 


Putting this in and interchanging limits of integration one is left with an 
integrand which is just the spectral function (*) 


GG = OE fa a fa (+ a? — 0°x(1— 2))- 


1 
sr CAC AT (St + 3t?a?) o4(a?)] , 
gi(o2) = — 4et, (“) fur de [ara + a? — o?x(1— 2))- 
2 _| 2 2 Ziel 2 OP + oc?) (1 20)? 2 
Jia + aqua) gg À + at) (6 + ot) — IE AED oa) 


h) Now note the following: for the 6-fanction to vanish one must have 
o?æ@(1 — x) — a? >0 


1) since o?a(1 — x) < o?/4 this implies a? < o?/4; 


2) for given a, æ=%[1 + (1 — (422/0?))*] 


is the region of vanishing of the argument of the 6-function. 
The ranges of integration thus reduce to 


07/4 3[1+(1-(4a?/0°))!/2] co 


foce da Joi 2(0 ak : 6(t —[o%a(1 — x) — a? ])[...]. 


2g(1 a) a*)# 
se 3[1-(1-@a2/03)!/2] 0 
The x integrations are easy to carry out since the square root in the integrand 
vanishes at the end points. This results in 


uè 
0° /4 
27073 (A lola 1/0: 
gz (02) = — ea(?) Jove (5-2) 01 (x?) x C - >) Où(œ?) 
————_——_—— le 
ERI die ola) dor 
(*) d(ar) = — d(a) an ER ({— a) for A 


a 
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The remaining integrations are elementary and yield, after the Born terms 


Ì 2 
PET 
“ tg 2 
( (ee kh 20% + à (1 gr sell 
= Le, Ce a ? 
a 3 DI fe ix Ix On) | 
ayers Heya) eee) 
Al CONTE 3 \u 3 \u Lt ca 


eue (AREE (0 SH RETI 
OUEN CAIN a 3dr \ to 


(2g2 + u? + we) (2q% + uw? + Sez) 
92(Q%, + w2) 


are included 


quo) er, (È Mat 


bo 


_ 


2 5 
la Sa 


F 2 
Git On, dn On | 


where g° = (o?/4) — u?. 


4. — Relation to relativistic form factors. 


The S matrix for a single scattering in an external electromagnetic field 
is given by (*) 


® x == / a | ex 
<p'|S|p> = — o u(p') L(g?) u(p) Aq) ; 


where p’=p-+q, 2 is the volume of the normalization box, and A%*(q) is 
the Fourier transform of dr (x). 
Now 
INI) E (EH 


Making the non-relativistic reduction to order 1/M by taking 


(2) Vu Vy FF Vy Vu = 26 ? 


1 
Ow cs 9j (Yury a VV y) ? 
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one finds (in the Breit frame-note again q? = q’) 


Lae + — (27) ; 
NEMI re NI 


GE f 
[= (2) — <a F,(¢)} y —ia:(q/\a) A (4%) + - OM ) , 


where 
AM (x) = a, exp Liga] , 


a, = (a, ig) . 


The above is the full S matrix to first order in the external field. Going now 
to an interaction representation where the time dependence is given by the 


full meson nucleon Hamiltonian and H, = H , one has 


electromagnet'e 


q/2|S|— q/2> = — Kq/2 (o ) di&|— g/2) = 


— = Inid(B, —- E; )cal2] [a %(x) dx | |— g/2) = 


= — 2ni0(E; — E;)o . 
Therefore, since 


<5 (p'—p—4q) = 


= (GE 2)g —io- (qd\a)f(4°) ; 


: q? 
ANT (9°) — sa F2) ; 
F,(g 
(BOIANO DA ‘ 


For the second equation, one can add the corresponding spectral functions 
however, for the first equation one must write 


Go?) qa G,, (6?) : 
f,(q?) fac: Sup OM 0° + @? |= =) 
dr? 
i [ 2 2 
MR) = DEMO G,(0 )do?, 
Ape 
j A 2M\G 
(4?) = fan (0? ) pees Fy | eo) 
Cie di 
4u 
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Thus if one uses a superscript (1) to indicate terms up through 1/M and a (0) 
to indicate the static limit one finds for the predictions of the static theory 


for G®(02) and Gy (0°) 


GB(6")sx, = gurl?) — Sa 


o? 
y(1) Lee NILO 
te (Oar = Gr (0°) gm 2% (0°); 


m(f\? (a\?(—@, , & + OR ao ts) 
AE ARIETE 


Solera ace eter | 
e 3 (EVE ceo 


We) Pa) (aa) a È 


CU) r=2e7, (“) | 
(Op 
(ge a Ge) Fa) 
: la (60) 2 | I 


k 


5 (2) (ee + wn) (297 + u?- Boz) — (q2 = yg 2). 
n ts FAURE + Wi) 


(0 etna) 20,1 à du + Bo 
qe è OR Gn) 34 GE il I | | 


The corresponding equations of dispersion theory obtained by expanding equa- 
tions (15.1) to (15.7) of reference one to order 1/M are 


E GIE (4 an Op dr at, 2 da) ce 
Hi, DATE i op) M 


.Jo (1) Pere) x È 2 (o Tipe + On ( Op tos de 
es Me VE x Or 


GD (0°?)pr. = 267; (È) (A) AS ) 2 (= ae fox dr ni aE | 


see tee E eye 295 ae Lo (a +4 
Oh On 2M | Lage 


a: (4) (€ 29m + On + 4°) — (gz + A a JT 
u (AU + à) 


(5 == On te-1 Wa Da) 4 207, Gn ale we 
: la À On Ia = Mi i 
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Fig. 1. — Spectral functions for the magnetic form factors calculated from dispersion 
theory (D.T.) and the static theory (S.T.). (0) indicates the static limit and (1) indi- 
cates first order in 1/41. 


Fig. 2. — Spectral functions for the charge form factor calculated from dispersion theory 
(D.T.) and the static theory (S.T.). (0) indicates the static limit and (1) indicates 
first order in 1/M. 


ns ¢ . amr Y(0) 2 y(1) (1) 7 Alici 
Rio lean dea replovie Ge (ERG) om and “also 


4 = À i 1 4 a a ry 
G,.5(0°)p obtained from the equations of reference (!). 
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5. — Summary and discussion of results. 


To summarize the results. 

a) One can calculate spectral functions within the framework of the 
static theory by a direct transformation of the form factors. 

b) The results obtained agree with the static limit of the results obtained 
from dispersion theory (except for the term proportional to 4(f/u)? — 3(a/w)? 
in g,(o*) in the static theory which must come from the neglect of pair effects ; 
it happens to be almost zero). 

c) The 1/M corrections coming solely from the addition of the kinematics 
of nucleon recoil are very nearly equal (*) to the corresponding corrections in 
the dispersion theory and the fact that they are so large raises serious doubts 
as to the applicability of the static model to the study of the electromagnetic 
structure of the nucleon. 


The physics of the situation is this, the nucleon likes to emit and absorb 
P-wave mesons, and indeed, if it is infinitely heavy can only do so. After the 
meson picks up a photon from the external fields it is no longer a P state 
relative to the nucleon, but has a great probability of being in a higher angular 
momentum state. If now one allows the nucleon to move and hence interact 
with mesons of higher angular momentum, even though the interaction is 
small, the amplitude of such a wave may be large enough so that the result 
may compete with P wave absorption. This is what was found above. 


The author would like to express his great indebtness to Prof. S. FUBINI 
for suggesting the problem and for his assistance throughout. He would also 
like to acknowledge a stimulating correspondence with Prof. 8. DRELL and 
to thank Prof. B. FERRETTI and the theoretical division at CERN for the 
hospitality extended to him during his stay. 


() Notice that the 1/17 Born terms are given exactly. 


RIASSUNTO) 


Le funzioni spettrali dei fattori di forma del nucleone si discutono nel quadro 
della teoria statica. Si includono anche correzioni per il rinculo 1/M e si mostra che 
tali correzioni sono prossime a quelle date dalla teoria dispersiv 
fondamente le funzioni spettrali indicando che l’applie 
questo problema è di dubbia validità. 


a e che alterano pro- 
azione del modello statico a 


(*) Traduzione a cura della Redazione. 
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On the Radiative Decay of Charged 7-Mesons. 


N. CABIBBO 


Istituto di Fisica dell Universita - Roma 
Istituto Nazionale di Fisica Nucleare - Sezione di Roma 


(ricevuto il 7 Gennaio 1959) 


Summary. The radiative decay m+—+p.*++y-+y is examined in the frame 
of the A—V theory of weak interactions. The modifications of the spectrum 
of u* emitted due to the structure of the u-meson are in part related to 
the process 79-y+y. These modifications turn out to be negligible. 


A recent experiment by C. CASTAGNOLI and M. MUCHNIK (1) has shown the 
possibility of measuring the spectrum of u* coming from the radiative pion 
decay: 


+ N Se ay 
TG se wae . 


In the present paper we shall attempt a detailed analysis of this process 
in the A—V theory suggested by FEYNMANN, GELL-MANN and MARSHAK- 
SUDARSHAN. 

If one deals with e.m. and weak interactions only to the first order, Lorentz 
invariance and gauge invariance impose strong limits on the form of the 
relevant S-matrix elements (2). Assuming universal Fermi interaction in the 
form: 

(1) H,, = V8G(Hy,av) (yap) ; a=}(1+y3) ( 


(1) C. CasragnoLI and M. Mucunik: to be published in Physical Review. 

(2) S. KamerucHI and 8. ONEDA: Nuclear Physics, 6, 114 (1958). 

(8) According to the hypothesis of non renormalization of the vector part in weak 
interactions, this part shall be written: G(u(x)y,av(x))J, (x) with the condition: 


Ia) = 0. 


cd 
co 


838 N. CABIBBO 


the matrix element will be a combination with coefficients depending on (pk) 


of the following terms 


A (à) 
M,= (qh) (u{ek + 2(qe)}au) , 
(2) WM, = (pk) (wean) , 


Gi (wy au”) c's, Dy Ke ee (u{kep Le e(pk)}au") 3 


k and e are momentum and polarization vector of the photon emitted, p 1S 
the four-momentum of the 7-meson, 4 and h the momenta respectively of the 
u-meson and of the neutrino, M is the mass of the 7, m the mass of the u; 
(pk) is the scalar product p,k,, e —e,y,. We use units in which #—= e 

Neglecting effects due to the internal structure of the z, the actual (weak 
and strong) interactions can be substituted by a simple derivative coupling 
between mesons and leptons (4). In this simple theory the only term which 
appears in the S-matrix is 977,. 


e Il 
(3) Sa (= Se ATC 
V 4a AV ky Zi 


the constant A is defined in terms of the matrix element for the normal 
mode rt—ut+ v 
(4) SO AG’ paw”) . 

In order to evaluate the structure effects on the + spectrum, we made a 
perturbation calculation of the S-matrix to the lowest order, making use of 


standard Feynmann-Dyson techniques, with the following interaction Hamil- 
tonian : 


(5) Jel’ = Ha + Ha "a I(MYP)P, + Ta 


a gere a ! O . 
The result is an addition of 977, and 977, terms to the principal 977, term: 


(6) S= a Eas A{M, + R[ WU, + WM} . 
AT 29 ky 


(4) G. E. Frarno and J. Tiomno: Anais Acad. Brasil. Cienc., 24, 245 (1952); B. 
Iorre and A. RUDIK: Dokl. Akad. Nauk SSSR, 87, 359 (1952). S.A. BLUDMAN and 
M.A. RuDERMAN: Phys. Rev., 101, 910 (1956). 


838 


ON THE RADIATIVE DECAY OF CHARGED 7-MESONS 839 


The constant A is the same that appears in eq. (3) and (4), while À is de- 
fined as 


— 
= 
= 


ay a (eens 
oi, = Fain wana 
y 


R= he a | Ie dt Dp, 
Jeemecmein/lemer"; ant 


where N is the nucleon mass. If we do not use a cut-off in the evaluation of 
the integrals, R turns out to be zero, and we fall back to the case of equa- 
tions (3) and (4). 

From equations (6) and (2) we get the following expression for the diffe- 
rential probability, integrated to give the spectrum of u* 


(8) W(E)dE = ia 


Haute AD ae ee AR 
Te Fe + 7) - 


+ RM(E,— E)[2(M — E)f(E) — Mg(E)] + 


4 3 
+ R sua P(H; E) ue, E° 


Î(E) and g(E) are the following functions: 


(E + P)M — m° 
= eyo eee LE 
age o8 (E — P)M — m2’ 


Here W° is the probability of the normal decay m+ > ut+v; Py = (M? — m?)/2M 
and H, = (M?-+m?)/2M, are the maximum values of the momentum P, and 
of the total energy Æ of the p. The coefficient of 977, in equation (6) is equal 
to the coefficient of 277,; this result rests only on the perturbation treatment. 
The coefficient of 277,, on the contrary, can be related by the following equa- 
tion (5) to the lifetime of the neutral r°-meson under the hypothesis of non- 
renormalization of the vector part 


Pym 


AT peri [+ 
IAT AEREO 10 17 1% 
(9) | Fe | ana)! a 0.3-10 = MeV, 


(5) V. G. VAKS and B. L. Iorre: Nuovo Cimento, 10, 342 (1958). 

In the Appendix we give a new derivation of a result due to VAKS and Iorre, based 
on straightforward application of reduction formulae, in which the role played by the. 
different approximations is clearly exhibited. 


Gel 
co 
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t and 7° are in the lifetimes respectively of z+ and r°-mesons. The sign of È 
cannot be determined in this way, but we assume it to be positive, as in à cut-off 
evaluation of eq. (7). At present we have only an upper experimental limit 
for 7°: 10-15 s (6), while the most reliable theoretical evaluation gives 5-107! $ (7). 
The correspondent values of R are R—0.14:10- MeV? and E=0.64-10-° MeV=2; 
from eq. (7), with a cut-off equal to the nucleon mass N: R—0.38-107$ MeV-?. 

An experimental determination of the u* spectrum could give some infor- 
mation on the value of À, and hence, through eq. (9), on 7°. However, as- 
suming k =10-* MeV-?, the terms in eq. (8) proportional to R and R? turn 
out to be negligible: the contribution of the term proportional to À? to the 
total probability of radiative decay is of order of 10-1W°, If we consider the 
coefficients of f(E) and g(£) in the R term, they turn out to be smaller by a 
factor 10-? and 105 respectively than the coefficients of the same functions 
in the first term, even for the more favorable values of E. Of course the 
coefficient of 977, may actually be different from that of 977,, but we think 
that our estimate will still give the right order of magnitude. 

Moreover the coefficients of f(#) in the first term and in the À term have 
nearly the same energy dependence. 


I am grateful to Prof. R. GATTO who suggested this problem and has given 
constant advice, to Prof. M. CINI for useful suggestions, and to Prof. C. CASTA- 
GNOLI for his interest in this work. 

I am also grateful to Dr. FRANCESCO CALOGERO for his active collaboration 
in the early stages of the work. 


APPENDIX 


The vector part of the matrix element YutVont | Tin > is related to the matrix 
/ 0. R à + = SO . 
element <yyo | rh) (5): from invariance arguments these matrix elements can 
be written in similar form: 


A Ge 
(10) NI = Cri B6(p — q — h — Baden el ; 
(10') Vout | Tin = 2B 0(p — k — Pe Mele 40 ,. ol =a ant ee 


(°) G. Harris, J. OREAR and S. TAYLOR: Phys. Rev., 106, 327 (1958). 
(?) M. L. GOLDBERGER and 8. B. TREIMAN: Nuovo Cimento, 9, 451 (1958). 
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In the following we will show that the factors B! and B? are the same func- 
tion B(pk). 

In (10): k! is the total momentum q+h of the leptons emitted. 

In (10’): e! and k! are the polarization vector and the momentum of one 
of the photons, e and k those of the other. In both (10) and (10'), A=4(k — #1). 
All other notations have already been used in the text. 

The tensor form of (10) is equivalent to that of 977, in equation (2), as is 
easily seen using the energy-momentum conservation: k + k!= p. 

B will be a function of the number & = (pk)/A?, where A has the dimension 
of a mass, and plays the role of a cut-off parameter. If we assume A = N, 
the values of € are in both cases very small, and we can assume B(É)=B(0)=B. 

In the Heisenberg representation the u+-meson field (7) and the e.m. 
potentials A, satisfy in our approximation: 


(11) (0 + mu(x) = 4(y,av(x)) I (2), 


(11) DrA,(0) = Ji) = e(J (x) + J'(x)) . 

The e.m. current is divided in an isovector part J, and an isoscalar part J®. 
In the derivation we use the following equations, which are true only in the 
lowest order of the constants e and @: 
(12) BE Je ay fe bee 
(13) v(x) = Vout(®) * 


By repeated application of reduction formulae to the matrix element 
uva, we get: 


:G | 
(14) (You | Tin > = aa np fato your | (vs@v(@)) JI) | at) exp [— iqa] = 
= na I fase Ire (x) |a+> exp[— i(q + h)x] = 


mina Che, | at dy CO Fay) TV (0) |) exp [— tht” — iky] = 
AT )2 


= = Gi mÔ(p — k — k1)l, aa QPAE I (- Le exp[— iad]. 


Here equation (13) was used to take the neutrino out of the state, and equa- 
tion (12) allows us to reduce a 7 product to a simple product. Moreover 
equation (14) is still satisfied if we commute J;(x/2) and J{(— a/2): 


G p 
(14') CUP ae) Mi 3 2m 0(p —k—k)le; 


EE 
fato 


a 54 - Il Nuovo Cimento. 
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The same procedure can be applied to the matrix element <yYou |Ttin}} and 
we get: 
al 


(15) Your | Tin = — e2n0(p — k — kh) ce, 
[are <O|J, (3) J, (= |x» exp[— dad]. 


Now, as in equation (11), we divide the current J, in an isoscalar current Je 
and the third component of an isovector current J®, of which J‘ is the 


(—) component. 

Under the generalized charge conjugation G = 0 exp [ixT,/2] these currents, 
and the vacuum and one meson states transform as follows: 

GIOIE JO, G\x> = —|n), 

(16) 
GIA Go = 102 
Therefore it is easily seen that the only non vanishing terms are those of the 
form <0|J°J®|z> and <0|J°J*|z>. From charge independence we get the 
following identities: 


TEE 


Comparing equations (14) and (14'), with equation (10), we get: 


EX 
— 
= 

~w 

ET 
wis 
ee 
= 
zx 
sg 
© | 
SS e 
à 
[=] 
N 
I 


18) —2a | di (0 2) (- 5) lat) exp[— iv4] = 


0 4 A À wv . 
a | dix (0|J; (- J; (5) |a) exp[— iw] = ep, A,B) . 


From (15) and (17) we get: 
(19) Dern AB? = fern A, + Zia AS) BHE 
Therefore bo —= 62" 


In order to derive equation (9) we need the lifetime of 7*- and 7°-mesons 
in terms of A and B: 


(20) ct = (W)-? =2M(mP,A)?, = 16M (ee BM). 


RIASSUNTO 


SI esamina il decadimento radiativo mt—ut+v—+ nello schema della teoria 
A—V per le interazioni deboli. Le modificazioni dello spettro dei ut emessi, dovute alla 
struttura del mesone x sono in parte messe in relazione col dat TOS vt 
Queste modificazioni risultano trascurabili. | bit 
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L/K-Capture Ratio in !°°Au. 


A. BISI, E. GERMAGNOLI (*) and L. ZAPPA 


Istituto di Fisica del Politecnico - Milano 


(ricevuto il 14 Gennaio 1959) 


Summary. — The L/K-capture ratio in orbital electron capture decay 
of Au was investigated by using a measurement procedure based on 
coincidence techniques. It was found for the transition to the 128 keV 
level of 15Pt: P,/P, = 5.5 + 0.9 which corresponds to a transition energy 
of 106 keV. In addition a determination of the half-life of 199Au and of 
the Æ conversion coefficient of the 97 keV transition was made. It was 
obtained 7,= (192-+5) d and xx = 8.4 + 0.5. The proposed decay scheme 
agrees quite well with previous results. 


1. — Introduction. 


The decay of !9Au to the excited levels of Pt has been investigated by 
several authors (+), and so the y-ray energies and the orbitals of the excited 
levels (Fig. 3) are well known. The obtained level pattern is in agreement with 
that deduced from the decay of Pt”, while there is only a partial agreement 
with the results of Coulomb excitation of Pt, which suggest that the 97 and 
30 keV transitions are interchanged. 

A few uncertainties still remain on the conversion coefficients and relative 
transition rates of the excited levels of Pt. Moreover the currently listed 
value of the half-life of Au (7, = 185 d) appears to be in contrast with the 
value (7, = 90+ 10d) recently obtained by CRESSMAN and WILKINSON (?). 

As far as concerns the L/K-capture ratio and the disintegration energy 
of Au no further progress has been made after the investigation carried out 


(*) Laboratori C.I.S.E. e Istituto di Fisica del Politecnico di Milano. 

) D. STROMINGER, J. M. HOLLANDER and G. T. SEABORG: Rev. Mod. Phys., 30, 
1958) 
) 


D. J. CRESSMAN and R. G. Witkinson: Phys. Rev., 109, 872 (1958). 
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by some of us (*) a few years ago. It is worthy of mention that in the above 
investigation the L/K-capture ratio was deduced from the measurement of 
L and K X-rays intensity. 

Errors in the result may have arisen from the fact that the L fluorescence 
yield values are rather uncertain, that in the decay scheme are present low 
energy y transitions which are highly converted in the Z-shell and whose con- 
version coefficients are not well known, and that the procedure was rather 
indirect. 

The present work was intended to determinate once again the L/K capture 
ratio by using a direct measurement procedure based on coincidence techniques 
already employed in analogous investigations (7). As a by-product of this 
research some data on the decay of Au and its half-life were obtained. 


ean KX-rays 2. — Experimental results. 

15Au was obtained after 
irradiation of Pt with deu- 
terons in the Amsterdam 
cyclotron and supplied as 
carrier-free AuCl,. It can be 
mentioned that the decay of 
Au was observed over a 
period of about 250 days and 
that its half-life resulted to 
be (192 + 5) d; moreover it 
was observed that the rela- 
tive intensities of K X and y 
peaks did not noticeably 
change in the same time in- 
terval. This fact points out 
that no other important acti- 
vities are present within the 
examined active specimens. 


6000 + 


4000 


Counting rate 


2000 
escape peak 


Channel tila Fig. 1. = KX.- “and. y-ray 


Is se 
0 10 20 30 40 50 60 spectrum from 1°Au. 


(8) A. Bist and L. ZAPPA: Nuovo Cimento, 12, 539 (1954). 


(4) A. Bist, E. GERMAGNOLI and L. ZAPPA: Nuclear Pia 
dora si yystes, 1, 593 (1956); N 
Cimento, 4, 764 (1956); 6, 299 (1957) y ( ); Nuovo 
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21. — Fig. 1 shows the y and AX spectrum, which was obtained with the 
help of a single crystal spectrometer (NaI(T1) crystal, of 3.75 em diameter, 
2.50 em high, optically coupled to a Dumont 6292 photomultiplier): during 
the present measurements the spectra were examined by means of a ten channel 
fast pulse analyzer. An intense peak whose energy (about 67 keV) corresponds 
to characteristic K X-rays of Pt and a y line at 97 keV can be easily noticed: 
the weaker peak at about 130 keV is a composite one and is mainly due to 
the fact that two K X quanta, which are emitted in the same disintegration, 
can both be detected by the crystal spectrometer, the resulting peak actually 
being at about 130 keV. Another consequence of the integrating features of 
the y spectrometer is the 
presence of a very weak A Is 
peak noticed at about Tor 
165 keV: this is due to 10 
coincidence between a A X 
quantum and a 97 keV 
y-ray. However the inten- 
sity of this peak is so low  05F 
that no quantitative con- 
clusions could be drawn 


T 


L 
in this case. 0 : is, 
> 
The ratio between the 0 05 1.0 
intensities of the summa- Fig. 2. — Determination of the J,,,/I,, ratio by the 


summation method. 


tion peak at 130 keV and 
the y line at 97 keV was 
measured after locating a point source of Æ5Au at different distances from 
the surface of the detector. The obtained results are given in Fig. 2. The 
ratio decreases linearly if the source-detector solid angle is decreased and be- 
comes constant at distances greater than 15 cm. This procedure allows the 
conclusions that a weak y line at about 128 keV is present within the sum- 
mation peak. 

Table I gives the results concerning the relative intensities of K X and 
y lines which where observed: the efficiencies of the detector and the escape 
probabilities of I K X-rays from the crystal were taken into account. 


TABLE I. 
Energy Relative intensity 
K X-rays 100 
97 keV 9.61 
128 keV 0.71 
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9°2, — All the peaks which are listed in Table I were used in coincidence 
measurements. In any case, only pulses belonging to one definite photopeak 
were allowed to trigger the coincidence circuit (2 us resolving time). In coin- 
cidence work with the 97 keV y-ray the triggering peak was filtered by means 
of Ag absorbers to eliminate the contribution of the tail of the intense K X 
peak: an analogous method was used when the triggering pulses belonged to 
the 128 keV line: in this way the above mentioned summation peak was com- 
pletely cut down. Table IT summarizes the results which were obtained in 
coincidence measurements; corrections for the probability of escape of I K X- 
rays from the crystal are included. 


TABLE II. 
anno Coincidence intensity/intensity of triggering peak 
peak si pa 
K X-rays 97 keV 128 keV 
128 keV (6.800110) 10 = a 
97 keV 0.141 + 0.010 _ — 
K X-rays 0.227 + 0.018 (1.42 + 0.12)-10- weak 


According to the decay scheme given in Fig. 3 the L/K-capture ratio for 
the orbital electron capture transition between the ground level of “Au and 
the 128 keV excited level of Pt can be deduced from the coincidence measu- 

rements with the 128 keV line used as 
triggering peak: 
195 195 
7ePt 79 AU 
ea 
VEC Px 


(1) E Ep 
1128 fo Sees 


Here P, and P, are the electron capture 
probabilities for the Z- and K-shells, w, is 
the fluorescence yield for the Æ-shell in 


o/2- 128 keV : 

je rey Pt (w, = 0.95) and e, is the counter effi- 
ciency for K X-rays (e, — 0.461 for our 
experimental arrangement). From the data 
of Table II we obtain 

D - 0 

1 
Fig. 3. — Decay scheme of 195Au. FT DEEE 
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which according to the calculations by Brysk and Rose () corresponds to 
a transition energy of 106 keV, if the transition itself is either allowed or first 
forbidden (AI < 2). 

As far as the 97 keV y transition is concerned, the Æ-shell internal con- 
version coefficient «, can be readily obtained from the coincidence measu- 
rements with the A X peak used as triggering peak. Actually we get 


where ey, is the counter efficiency for 97 keV y-rays. For a thick Nal crystal, 
€ can -be regarded as equal to e, and consequently from 7,,/1,,, — 16 +1 
we obtain 

C= 84 0D 


Finally, if the 97 keV y line is used to trigger the coincidence circuit, 


Thos SII | IP I a ¥ ) IE, 
lee po NGS ID AI 


where w is the branching probability towards the 120 keV excited state of 
195Pt, (1— &) is the branching probability towards the 97 keV excited state 
of 1%Pt and the indexes 1 and 2 refer to the two considered transitions. The 
relationship (3) shows that the L/K-capture ratio for the orbital electron capture 
transition between the ground level of Au and the 97 keV excited level of 
195Pt cannot be directly deduced. Nevertheless if the experimental result 
JT, xlTor = 0.141 is inserted in (3) we get 


(3) 


Pe Pr 
099220287 L (1 di 
: (Fx) r { n) (a) 


Taking (P,/(P,+P,)),=0.154, which was deduced above, and (P,/(P,+P,)).= 
= 0.49, which was calculated from the results by Brysk and RosE for an 
allowed or first forbidden transition (AI <2) whose energy is about 
(106 +128 — 97) keV = 137 keV, we find 


[EOS 


It may be concluded that the two considered electron capture transitions take 
place with approximately the same probabiliy. It has been assumed that 


(5) H. Brysk and M. E. Rose: Rev. Mod. Phys., 30, 1169 (1958). 
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there is no appreciable branching between the ground states of Au and !95Pt: 
this was suggested by the fact that the observed intensity of K X-rays is 
completely accounted for by the contributions due to the K internal con- 
version of the 97 keV y-ray and to the 106 keV and 137 keV K electron cap- 


ture processes. 


9°3. — The internal consistency of the above conclusions was checked by 
the results obtained from the observed intensity of K X-K X coincidence 
pulses. If the probability of the 128 keV y-transition is supposed to be so 
small, that its contribution to K X emission can be neglected, we get 


Mier 9 XK ali 
—AO0rErx 


Tx dea re) CEA a) [u(Pxl(Pi+Px)}+ (À — fi) (Pri(P Po) i 


(4) 


The total conversion coefficient « of the 97 keV y transition was considered 
equal to 9.9, in accordance with the results by DE SHALIT et al. (5); if this 
value of « and the above experimental results are used, we have 


Tak - 0.20 À 


K 


which is in satisfactory agreement with the value 1,,/1,= 0.227 + 0.018, 
obtained from coincidence measurements. 


3. — Conclusion. 


Fig. 3 shows the decay scheme of Au as suggested by the results of 
Section 2; it appears to agree completely with previous results. The disin- 
tegration energy to the 128 keV level is slightly different from that given in 
ref. (*), although the new P,/P, value (5.5) is much greater than the old 
one (0.58). This is because in the low energy region the transition energy 
varies very slowly against P,/P,. The disagreement between the new P,/P, 
value and the old one can be explained by recalling that measurements in- 
volving L X-rays are made difficult from the necessity of estimating the counter 
efficiency and L fluorescence yields; moreover the old value was obtained 
from the difference of two near terms. 

By assembling the obtained information on disintegration energy, branching 
ratio and half-life, the log ft values for the transitions to the 128 keV and 
97 keV levels were found respectively equal to 5.7 and 6.8. 


(6) A. pe SHaLiT, O. HuBer and H. SCHNEIDER: Helv. Phys. Acta, 25, 279 (1952). 
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These values are typical of allowed (A7 — 0.1; no) or first forbidden 
(AJ = 0.1; yes) transitions, and together with the X conversion coefficient 
value for the 97 keV transition obtained in Section 2, represent a set of results 
consistent with the orbital assignement to the levels involved in the decay, 
made by previous authors and included in Fig. 3. 

An additional information can be deduced after combination of data given 
in the previous paper (*) and of the present results. The intensity of uncon- 
verted 30 keV y-rays, measured by means of a proportional A--CH, counter, 
was equal to 1.42, taking the AK X intensity equal to 100. If the probabilities 
of the two electron capture transitions towards the 97 and 128 keV levels 
are supposed to be equal, as suggested from the present data, from 7,;,/1:5—6.8 
one gets that the total conversion coefficient of the 30 keV y transition is 
= 30. This is consistent with the assumption of a M1 transition. 


RIASSUNTO 


Viene studiato il rapporto tra cattura L e cattura X nel decadimento dell’!95 Au 
usando un procedimento di misura basato su tecniche di coincidenza. Si è ottenuto: 
P,/Px= 5.5 +0.9 per la transizione che conduce al livello di 128 keV del #5Pt. A 
questo valore corrisponde una energia di transizione pari a 106 keV. Sono stati, inoltre, 
misurati il periodo di dimezzamento dell’ Au (7, = (192 + 5) d) ed il coefficiente di 
conversione A della transizione da 97 keV (xx = 8.5 +0.5). Lo schema di decadi- 


x 


mento proposto è in accordo con i risultati di precedenti ricerche. 


NOTE TECNICHE 


Performance of a Large Area non Focusing Cerenkov Counter 
and Absolute Yield of Cerenkov Light. 


G. CAGLIOTI (*), R. CERVELLATI (*) and L. MEZZETTI 


Istituto di Fisica del Universita - Roma 
Istituto Nazionale di Fisica Nucleare - Sezione di Roma 


(ricevuto il 23 Febbraio 1959) 


Summary. — The design and the performance of a large area, compa- 
ratively thin toluene Cerenkov counter are described. The velocity 
resolution of non focusing Cerenkov counters and its basic limitation, 
due to the Landau fluctuations in the emission of fast knock-on electrons, 
are discussed and compared with those of scintillation counters. The 
absolute yield of the Cerenkov light is evaluated from the measurements 
and found to be in good agreement with the theoretical yield. 


1. — Introduction. 


In the last few years focusing and non-focusing Cerenkov counters have 
been extensively used for the detection of fast charged particles (*). 

The measurement of the velocity f of a charged particle can be performed 
with great precision with a focusing counter, when the incident beam has à 
very small cross-sectional area and angular divergence (?). 

Non focusing counters can be used also when these conditions are not ful- 
filled, using the total amount of light produced as a measurement of f (35). 
A further advantage of a counter of this type is that it can be converted into 


(*) Now at the Comitato Nazionale Ricerche Nucleari, Roma. 

(1) J. V. Jeviny: Cerenkov radiation and its applications (London, 1958). 
(2) R. L. MATHER: Phys. Rev., 84, 181 (1951). 
(8) P. Bassi, A. M. BranCHI and C. ManpucHI: Nuovo Cimento, 9, 861 (1952). 
(4) M. Manpò: Nuovo Cimento, 12, 5 (1954). 

(5) J. R. WincKLER, E. N. MircHELL, K. A. ANDERSON and L. PETERSON: Phys. 
Rev., 98, 1411 (1955). 
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50 mm 


“i 


50 mm 
Fig. 1. — Schematic arrangement of the Cerenkov counter. Some of the materials are 


indicated: Gl = Glass, Al = Aluminum, Re = Reflector (polished Aluminum or white 
coating), Po = Polyethylene, Br = Brass, Sr = Soft rubber, Lu — Lucite. 
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a multichannel velocity selector by applying conventional pulse height analysis. 
Serious drawbacks, however, are represented by the relatively small amount 
of light usually available, which gives rise to poor photomultiplier statistics, 
and the fundamental limitation set by the Landau fluctuations. “an 
In this paper we want to discuss quantitatively such possibilities and limi- 
tations on the basis of a few results obtained with a comparatively thin, large 
area Cerenkov counter built in this laboratory for measurements on Cosmic 
Ray particles. The counter was originally designed to select fast Cosmic Ray 
secondaries travelling upwards, from a very large cross-sectional area. 


2. — Experimental apparatus. 


The Cerenkov counter (0) is shown in Fig. 1; it is enclosed in two cylin- 
drical Aluminium boxes 47 em in diameter, separated by a 5 mm thick circular 
glass window. 

One of the two boxes hides seven Dumont 6364 photomultipliers (PM). 
The PM are mounted in a close-packed hexagonal arrangement; their cathodes 
are kept in optical contact with the glass window by the spring action of the 
silvered steel contacts providing the electrode voltages, mounted on a lucite 
frame. 

The other box is completely filled, to a thickness of 10 cm, with about 
12 liters of pure toluene. 

The choice of the Cerenkov material was based on the requirement of max- 
imum light output for a given energy loss. For ultrarelativistic singly charged 
particles one can use as figure of merit M the quantity: 


436 
M = 
oZ/A ( 


= photons g-' em? , 


which gives the number of photons between 3000 and 7000 À generated in 
a material of density o, index of refraction n, and effective mass to atomic 
number ratio A/Z. For lower energies, the calculated light yields for a few 
transparent materials are plotted in Fig. 2. Good light transmission and weak 
scintillation (5) were obviously also required and checked. 

The sensitive region of the liquid is limited laterally by a hexagonal re- 
flector whose cross-section matches the array of the seven photocathodes. 

The seven PM anodes are connected to an open ring of heavy wire which 
forms a line of approximately 200 Q impedance correctly matched at one end 
(« fast » output). A «slow» output, more suitable for accurate pulse height 
analysis, is obtained by mixing the dynode signals at the common plate of 
seven preamplifiers. 

The anode supply voltages of the seven PM were chosen so as to give equal 
pulse for equal amount of light collected by the photocathode. Standard light 
pulses were obtained from a small plastic scintillator with a weak Polonium 
a-source, placed in turn near the center of each photocathode. 


(°) R. Mapry and L. LEIPUNER: Nucleonies, 14, no. 4, 51 (1956). 
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The «slow » signal was fed either to a threshold discriminator or to a nine 
channel differential discriminator through a linear amplifier (7, = 0.08 us, 
maximum gain = 80 db). 


T 


6:10 


410 + 


Number of photons emitted by a 
Cerenkov threshold 


2:10 


y meson slowed down fromy fol 


15 20 25 30 35 


Fig. 2. — Calculated number of photons (between 3000 and 7000 A) emitted by a 
u-meson slowed down, in various transparent materials, from initial energy H = m,c?y 
to Cerenkov threshold, vs. y- 


In both cases the discriminator was gated by the twofeld (or threefold) 
coincidence of a Cosmic Ray u-meson vertical telescope AB (or ABD). 

Two kinds of measurements have been performed. 

In the first series of measurements we tested the counter for its « anti- 
directionality » (+): to this end, the bottom plate was covered with black velvet, 
and the lateral reflector consisted of a mirror polished sheet of aluminum alloy. 

In the second series of measurements we tested the counter for optimum 
pulse height resolution: to this end, the lateral walls and the bottom plate 
were covered with a white coating of TiO, dispersed in water glass on a white 
enamel background. The Cerenkov counter was usually operated with the 
photomultipliers looking down (Own). This arrangement, in which the light 
undergoes at least one diffusion before reaching the photocathodes, was found 
to be the most convenient, because of the greatly improved uniformity of light 
collection, in spite of the loss of light involved. 


8. — Results. 


31. Antidirectionality tests. — The results of the antidirectionality test are 
summarized in the plots of Figs. 3 and 4. 

The main features of the curves in Fig. 3 are in agreement with what one 
expects. We would only like to emphasize that a high efficiency and a good 
antidirectionality can both be obtained with a comparatively thin radiator 
(8.7 g/cm? in our case). 
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To investigate the origin of distributions ITI and IV we performed the measu- 
rement indicated in Fig. 4 (curve V). 


À / 
\ 
| 
1.00 
n 
G 
0.80 Pb 
UA) 
060} | -—-with lead 
= — without lead 
IS — 20cm 
040- = 
tu 
020 
ae 
0 20 40 60 80 
Fig. 3. — Efficiencies of the antidirectional Cerenkov counter vs. minimum pulse 


height N accepted by the integral discriminator, for the telescope shown in the inset. 


A glance at curves IV and V shows that approximately equal amounts of 
light going upwards are somehow produced in the toluene and in the glass 


window. 
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N, Dr a are 
Fig. 4. — Efficiency vs. minimum pulse height N. Curve V was obtained with the 
arrangement shown in the inset, i.e. with counter emptied and with a black lucite 
plate in optical contact with the glass window so as to absorb completely the Cerenkov 
light produced in the glass. Curve IV of Fig. 3 is reproduced for comparison. 


Actually it is a known fact (*) that most of the substances used for Ce- 
renkov counters do scintillate when they are crossed by ionizing particles. 


+ 
10 
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D.C. measurements performed in this laboratory indicate that our toluene and 
our glass have a scintillation efficiency, for 5 MeV «-particles, about 10-2 and 
5:10 times that of a standard p-terphenyl in toluene solution. This explains 
the relative positions of curve IV and V. The long tails towards the larger 
pulse heights can be qualitatively accounted for in terms of small knock- 
on showers. 


3'2. Pulse height resolution measurements. — Fig. 5 shows an overall Ciown 
distribution obtained by combining measurements performed with a narrow 
beam of u-mesons over a period of several weeks (*). The individual distri- 
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Fig. 5. — Overall differential pulse height distribution obtained with the « white » 
Cerenkov counter and the telescope shown in the inset. 


butions were not exactly reproducible because of a slight long-range instability 
of the electronics, including the calibration equipment. For this reason the 
spread of the combined distribution is somewhat larger than that of the indi- 
vidual runs. 

In Fig. 6 are plotted two typical distributions obtained with the white 
Cerenkov counter looking up and down and a lead filter 20 cm thick. 

As expected, the most probable pulse height is larger for the disposition 
C. than for Cawn, While the relative width of the Ci distribution is smaller 
than that of ©. From the ratio of the most probable pulse heights we 
obtain a value of about 0.7 for the average light diffusion coefficient, in agree- 
ment with direct spectrophotometric measurements. 


(*) The lead filter used was sometimes 10 and sometimes 20 em thick, but its 
thickness did not appreciably affect the results. 
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The counter was tested also for uniformity of response over its large area, 
by shifting it horizontally up to about 10 em with respect to the telescope. 
The Crown response was found to be uniform within the limits set by the men- 
tioned instability of the electronies, i.e. better than + DE 
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Fig. 6. — Differential pulse height distributions obtained with the « white » Cerenkov 
counter and the telescope shown in the inset. 


4. Landau type fluctuations and resolution. 


The velocity resolution of a non focusing Cerenkov counter is limited by 
many effects, and ultimately by the Landau fluctuations in the emission of 
fast knock-on electrons. 

The latter have been discussed qualitatively by LINSLEY and HORWITZ (°) 
who showed that for a sufficiently thin Cerenkov counter the « theoretical » 
resolution is better than that of a scintillator of the same thickness, because 
of the existence of a threshold in the Cerenkov emission. 

According to LINSLEY and HORWITZS criterion our counter cannot be con- 
sidered «thin » and its behaviour is therefore expected to be similar to that 
of a scintillator; this is in agreement with the shape of our experimental distri- 
butions, showing the typical Landau asymmetry. The contribution of this 
type of fluctuations to the fractional total width at half maximum can there- 
fore be estimated directly from Landau’s treatment (8), giving, for the average 


(7) J. LinsLey and N. Horwirz: Rev. Sci. Instr., 26, 557 (1955). 
(8) L. LAnDAU: Journ. Phys. (USSR), 8, 201 (1944). 
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u-meson crossing our telescope (see eq. (4)): 


A 
(1) Ga = 0.185 (*). 


One may notice that the relation between the error in the measurement, 
for example, of the quantity y = H,./moc? and the fractional error AH/H in 
the measurement of the light yield is in favour of the Cerenkov counter, in 
a range of y-values depending on the index of refraction, as indicated in Fig. 7. 


A 
10} 2Y 
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BL 
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=== Cerenkovini= Nes 
A 
0 27 Pash, 1 LE 1 1 Y 
1 2 3 4 5 6 
Fig. 7. — Calculated fractional error Ay/y due to. a fractional error AH/H in pulse 


height, vs. y= (1—B?)-? for a scintillator and for typical Cerenkov counters. The 
scintillation pulse is assumed to be proportional to the specific energy loss; the slowing 
down of the particle within the counters is neglected. 


To exploit this fact it is of course necessary to reduce the fluctuations due 
to photoelectron statistics to a sufficiently small amount, by improving the 
collection and photomultiplier conversion efficiency. Such improvements do 
not seem beyond the reach of foreseeable technological developments. 


(*) A calculation using the gaussian approximation (°) leads to 


A 1 Ay = 
= = 024, — — 0:26: 
H Cerenkov H Fcintillator 


(9) B. Rossi: High Energy Particles (New York, 1952), p. 31. 
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5. — Absolute yield of Cerenkov light. 


It is interesting to try to evaluate from our results the absolute yield of 
the Cerenkov light per particle, in order to compare it with the theoretical 
yield given by the Frank and Tamm formula. Estimates of this kind were 
performed firstly by Cerenkov himself (1°) (luminescence of a strong f-source 
in benzene), and later, through a somewhat indirect method, by BELCHER (1!) 
(again with B-rays in water), and by ASCOLI-BALZANELLI and ASCOLI (2?) (with 
Cosmic Ray u-mesons in gases and vapours). 

To do this we have estimated the absolute average number g of photo- 
electrons per particle in two ways: 


A) From the measured size of the most probable voltage pulse height. 
B) From the relative width of the pulse height distribution. 
A) The first method consists obviously in applying the equation 


CH 


= Haven à 
eG, G, 


(2) q 


where C is the average output capacitance, G, is the PM gain, G, is the total 
gain of the amplification chain, and H is the most probable pulse height as 
measured (in volts) by the discriminator. 

For G, we took the average gain of the seven PM computed from the 
nominal tube characteristics; G, was calibrated by means of a mercury switch 
pulse generator. 

Inserting the numerical values we obtained: 


4 = 135 photoelectrons/partiele . 


To obtain the number of Cerenkov photoelectrons we must subtract the 
contribution due to the scintillation of our radiator, evaluated from the curves 
of Fig. 3. The final result is 


(3) {oerenkor = 106 photoelectrons/particle . 


B) The second method requires that the contribution a) to the measured 
total spread due to the PM statistics be isolated from the other contributions, 
such as: b) energy and e) angular spreads of the impinging particles; d) Landau- 
type fluctuations of Cerenkov light; e) non-uniformity of the light collection 
efficiency and photocathode sensitivity; f) statistical fluctuations in the number 
of photons collected. 

We assume that the relative width of the pulse height distribution (full 
width at half-maximum 4, divided by the most probable pulse height H) is 


(CO) RIPA CERENKOV: Compt. Rend. Acad. Sci. USSR, 21, 116 (1938). 
(7) E. M. BeLcnER: Proc. Roy. Soc., A 216, 90 (1953). 
(12) A. AscoLI-BALZANELLI and R. AscoLi: Nuovo Cimento, 6, 1392 (1957). 
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made up by partial widths due to the above effects : 


ee SO Oey 


The most important partial width in the sum at the right hand side of 
eq. (4) is that due to the Landau-type fluctuations, and has been discussed. 
in the preceding section. The second important is the contribution f) which 
can be calculated as a correction with a successive approximation procedure 
assuming an average cathode efficiency of 11%. Taking into account also 
the other small contributions we obtained 


2 
(2) = (0.34)? — (0.21)? = 0.072. 


The effect of photomultiplier statistics can be represented with the for- 


mula (*) 
ASEM 2:30) 
(5) (7) ( > ) 


a Oro —1)| 7 
where q is the number of useful photoelectrons, f is the ratio of q to the 
number of impinging photons, 6 is the average dynode multiplication factor 
times the electron collection efficiency, 6; is the analogous quantity for the 
first dynode (which is usually operated at a higher voltage than the others), 
and e is a parameter connected with the supernormality of the statistics of 
the multiplication process, ranging, according to the various authors, from 
1.0 to 1.5; for the latter we take the value 1.3 (*). 

Inserting the appropriate numerical values and correcting as before for 
the scintillation photoelectrons, we obtained: 


(6) Qrerenkov = 79 photoelectrons/particle . 


It should be emphasized that this result is affected by a rather large un- 
certainty (and is probably too low) because of its critical dependence on the 
estimates of the individual partial widths and specially on the value of the 
measured total spread. For example, a change of the latter from 0.34 to 0.33 
would increase the value of g by about 10%; furthermore the spreads due to 
instability of the electronics and to non-uniformity of light collection and 
photomultiplier cathode efficiency, which have been disregarded, could easily 
account for an increase of several percent each in the final result. 

On the other hand, the number of Cerenkov photoelectrons to be expected 
can be calculated starting from the Frank and Tamm formula. The proce- 
dure is straightforward, but requires the knowledge of the optical properties 
of the whole apparatus used as a function of wavelength. We write approx- 
imately: 


270 1 he 
(7) Iverenkor = Faz (1 si Lg | A t(2)R(2) da , 


(*) For a complete discussion see U. AMALDI, jr.: Thesis (Rome, 1957). 
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where B is the velocity of the impinging particle, n is the (constant) index 
of refraction of our radiator, L is the radiator thickness in em, and g = 0.53 
is the geometrical collection efficiency for the Cerenkov photons. In the in- 
tegral {(2) is the attenuation suffered by light of wavelength / along its average 
path through the toluene and the glass up to the photocathodes, R(A) is the 
absolute spectral sensitivity of the PM used (S4) as given by the manu- 
facturer (*), and he//® is the energy spectrum of the Cerenkov light. #(4) has 
been obtained by spectrophotometric measurements; g has been calculated 
taking into account in an approximate way the multiple diffusion of the light 
on the white walls of the counter. To this end, the diffusion coefficient was 
estimated as mentioned in Section 3. 
In this way we obtained 


(8) Querenkov = 102 photoelectrons/particle . 


Disregarding the factor {(2) one would obtain, instead, 265 photoelectrons/par- 
ticle. Considering the uncertainties involved in the evaluations Of Qéerenkoy the 
agreement between the three values obtained is remarkably good, and may 
be considered as a further check of the absolute intensity of the Cerenkov 
light as given by FRANK and TAMM. 


6. — Conclusions. 


In conclusion, we would like to stress the following points: 


1) The Cerenkov counter described in this paper can be used as a high 
efficiency antidirectional detector, or as a good resolution non focusing device; 
in both cases good uniformity of response over its large area was obtained. 


2) The contribution of the Landau-type fluctuations to the overall spread 
of the pulse height distribution is found to be comparable with that of photo- 
electron statistics. 


3) The absolute yield of the Cerenkov light is experimentally found to 
be in good agreement with the theoretical one. 


() Differences between individual phototubes are at least partially compensated 
by the fact that the light is distributed among seven photocathodes. 


RIASSUNTO 


Vengono illustrati la costruzione e le prestazioni di un contatore di Cerenkov a 
toluolo di spessore relativamente piccolo. Viene discussa la risoluzione in velocità dei 
contatori di Cerenkov non focheggiati e la sua limitazione essenziale dovuta alle flut- 
tuazioni di Landau nell’emissione di elettroni knock-on veloci. Viene effettuato un 
confronto fra l’intensità assoluta della luce di Cerenkov ottenuta sperimentalmente 
e quella calcolata a partire dalla formula di Frank e Tamm. L'accordo fra i v 
teorici e sperimentali è buono. 


alori 


860 


IL NUOVO CIMENTO Vor. ly NG 16 Marzo 1959 


Experiments on a Ion Source 
for an Electromagnetic Isotope Separator. 


F. Bist and B. DE MICHELIS 


Istituto di Fisica Sperimentale del Politecnico - Milano 


(ricevuto il 24 Febbraio 1959) 


Summary. — The paper describes a hot cathode, magnetic ion source 
for an electromagnetic isotope separator. The source, useful both for 
gaseous and solid charge materials, gives a total ion output of about 
20 mA. Measurements on discharge characteristics, extracted currents 
and source efficiency are reported and the effects of introducing reflecting 
electrodes in the discharge chamber are discussed. 


1. — Introduction. 


We have developed a ion source to operate in an electromagnetic isotope 
separator. The choice of the type of source was determined by the usual 
requirements: 

1) maximum ion output not lower than 10 mA; 
2) low operating pressure; 
3) good performance both for gaseous and solid charge materials. 

To satisfy these requirements we studied a hot cathode, magnetic ion 
source with radial ion extraction (+3). The behaviour of our source is satis- 
factory for good efficiency and high ion output: in particular, a considerable 


increase in efficiency has been obtained by the introduction of two reflecting 
electrodes in the discharge chamber, as is described in Sect. 6. 


() H. Herr: Zeits. Phys., 120, 212 (1943). 

(2) R. BERNAS: Journ. Phys. et Rad., 14, 34 (1953). 

(3) J. KIstEMAKER, P. K. Rot, J. SCHUTTEN and C. DE VRIES, in M. L. SMITH: 
Electromagnetically Enriched Isotopes and Mass Spectrometry (London, 1955), p. 10. 
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2. — The source. 


The discharge chamber is cut from a block of graphite and supported by 
a stainless steel tube (Fig. 1). Two copper rods support the tungsten filament 
(diameter — 1 mm) and serve the current supply. The filament is connected 
to the conical extremities of the copper rods by two invar nuts: the negligible 
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Fig. 1. — Ion source: a) graphite box; b) copper rods; c) invar nut; d) filament; e) gra- 
phite furnace; f) extraction electrode; g) ion collector; h) magnet; à, 1) reflecting elec- 
trodes at cathode potential. 


thermal dilatation coefficient of invar assures a good electrical contact between 
rods and filament, even at high temperatures. The exit slit on the front plate 
of the chamber was shaped as a Pierce electrode (ey 

The graphite furnace for evaporation of solid charge materials is electrically 
heated and can be easily replaced by a system for introducing gas in the 


(4) J. R. Pierce: Journ. Appl. Phys., 11, 548 (1940). 
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discharge chamber. When solid charge materials are used, it is convenient 
to maintain the discharge chamber at a higher temperature than that of the 
furnace, in order to prevent material from being deposited in the source: to 
obtain this condition an electric resistance (200 W) is inserted in the walls 
of the chamber. Some parts of the system are water-cooled to maintain the 
O ring seals at low temperature, while the source is operating. 

The graphite extraction electrode is shaped as a Pierce surface for a col- 
limated ion beam. The dimensions of the slits are respectively 2x50 mm 
for the exit slit on the front plate and 455 mm for the extraction slit. The 
distance between the two slits can be varied: most measurements were taken 
at a distance of 8mm. The ions are collected in a Faraday cage which is 
always at a positive potential with respect to the extraction electrode, so as 
to avoid errors in current measurements due to secondary electrons from the 
collector. The magnet for the alignement of electrons emitted from the fila- 
ment gives a maximum induction of 800 G; the gap is 170 mm. 

The operating conditions used during the measurements are reported in 
the following Table: 


TABLE I. 
Are voltage (50200) V 
Are current (0.1+5) A 
Filament heating current (40 —60) A 
Magnetic field (100 = 800) G 
Extraction potential (2—7) kV | 
Furnace heating power | 500 W 


The study of the characteristics of the ion source was made with argon. 
We were interested to know the pressure in the discharge chamber; therefore 
the system for gas introduction was supplied with a Pirani gauge head or 
a ionization gauge head for different ranges of pressure. In experiments 
with argon the source ran continuously for several hours with great stability, 
good efficiency and high ion output. Periodical dismounting and cleaning of 
the source was necessary; the operation is very easy due to the simplicity 
of assembly. Some preliminary experiments with zine gave satisfactory results. 


3. — Measurements on the discharge. 


In Fig. 2 are reported some current-voltage characteristics of a discharge 
in argon at different pressures. At low discharge currents the are voltage is 
generally less than 40 V and the are current is mostly carried by primary 
electrons emitted from the filament. When the arc current exceeds the satu- 
ration current from the cathode, higher are potentials are necessary to increase 
the current, to which the secondary electrons in plasma contribute. The curves 
show a decrease of are voltage at low currents when pressure is increased; 
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this fact may be due to the higher collision probability that allows the same 
values of are current with lower accelerating potential for electrons being 
reached. 


Vv 
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Fig. 2. — Current-voltage characteristics of a discharge in argon for different pressures 
(magnetic induction = 630 G; filament heating current = 58 A). 


Fig. 3 shows the effect of increasing pressure on are current: the curves, 
at a constant arc voltage and different values of filament heating power, show 
the existence of a minimum pressure (5) required to maintain the discharge; 
that for our source is about 10-* mm Hg. 
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Fig. 3. — Are current as a function of discharge pressure for different filament heating 
currents (magnetic induction = 630 G; anode voltage = 170 V). 


4. — Measurements on ion current. 


Due to the low extraction potentials, the ion currents collected on the 
Faraday cage are low. However the total ion current, including the current 
collected on the extraction electrode, can exceed 10 mA. Fig. 4 shows 


(5) W. WaLcHEer in J. Kocn: Electromagnetic isotope separators and applications 
of electromagnetically enriched isotopes (Amsterdam, 1958), ch. XIII. 
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that the collected saturation current increases at increasing extraction po- 
tentials: this means that measured currents are space charge limited. With 
an extraction potential of 6kV at an extraction distance of 8 mm we obtained 
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Fig. 4. — Ion currents at collector as a function of total extracted current for different 
extraction voltages (magnetic induction —630 G; distance between the electrodes=8 mm). 


a ion current of about 8 mA with a ratio 1:1 between the collimated current 
and the current collected on the extraction electrode. Obviously the value of 
collimated current and the extraction ratio can be greatly improved with 
higher accelerating potentials valuable in the electromagnetic separator 
operation. 

In Fig. 5 the saturation values of collected currents are plotted as function 
of extraction potential; the data fit well the Langmuir formula (f) 


. Ju: Al 
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where U is the extraction potential, A,, the area of the emitting electrode, d the 
distance between emitting and accelerating electrode and M the mass number 
of ions. G is a constant depending on the geometric configuration of the 
extraction system. 


(6) I. Lanemuir e K. T. Comeron: Rev. Mod. Phys., 3, 327 (1937). 
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In measurements of collected currents as function of the distance d between 
the electrodes, the Langmuir formula was well verified except In case ote 
distance, probably due to the deformation of the plasma meniscus at strong 
electric fields. 


e È 
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Fig. 5. — Collected current at saturation as a function of extraction potential (magnetic 
induction = 630 G; distance between electrodes = 8 mm); % experimental data; 
— — — theoretical curve. 


5. — Measurements on efficiency of the source. 


It is of interest to operate a high efficiency source, chiefly in separating 
radioactive or rare substances. We have therefore measured the efficiency 
of the source, usually defined as 


7] 
a ; 
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where % is the ratio of the number of ions per second emitted from the source 
to the number of neutral atoms per second emerging through the exit hole. 
This last term can be calculated when pressure p and temperature 7, in 
the discharge chamber, and the area A, of the exit slit are known (5). Then 
for a gas of atomic weight M and an extracted current i, we have: 


ML 90 VI M , 
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From the data for our source, when operating with argon, we obtain: 


where è is expressed in ampere and p in um Hg. In the calculation we have 
used a discharge temperature 7, of 700 To 
With regard to the total current è emitted from the source slit, the theory (5) 
enables us to calculate the total 
efficiency of the source as: 


à M - E%o 
Ni L Je O az ’ 
DL 
where j, is the electron current En 
density, o the ionization cross Ir 
section for the gas, V the volume 
of plasma, A the total limiting | | ‘ Hg 
area of plasma, and A,, the area of 0 1 2 3 4 5 Jo € 
the ion extraction meniscus. This | 1 
relation is independent, at first ap- Fig. 6.— Total source efficiency e, and source 
proximation, of pressure in the efficiency for collimated beam e, as a function 
source. If we consider only the of pressure (magnetic induction 800 G; extrac- 
current i. of the collimated beam, tion potential = 5 kV; distance between elec- 
D j € 
we have: trodes= 8 mm). 
MALTA, 
Vi — = 
7» © par oe ? 


the efficiency of the source then decreases with pressure. 

The experimental curves for e, and e, in our sources are reported in Fig. 6 
as a function of pressure. The total efficiency of the source can reach the 
value of about 4.5%. 


6. — Effect of reflecting electrodes. 


To increase the efficiency we introduced two reflecting electrodes i and 1 
(Fig. 1) at the same potential as the cathode (?*); the electrode just behind 
the ane increases the primary electron current that enters the main 
discharge chamber, while the other electrode compels the electrons to oscillate, 
and so lengthening their path in gas and enhancing collision probability. 
With such an arrangement we obtained an overall efficiency higher than 
25%: the total extracted currents can reach (20-25) mA. 

Total efficiency curves e, for different operating conditions are reported 
in Fig. 7: the effect of introducing reflecting electrodes is clearly visible. When 
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the source is operating with reflecting electrodes, efficiency over 25 % can be 
obtained by increasing the filament heating power: however the curves show 
a marked dependence on pressure, 
which can probably be attributed 
«x to the variation of plasma meniscus 
x area A,,, at increasing are current, 
and to the variation of electron pri- 
mary current, due to the ion bom- 
bardment of filament. 
10 x At high arc currents an intense 
effect of filament sputtering takes 
place, which shortens the filament’s 


i S i life and causes the metal vapours to 
2 O be deposited on insulators. To ob- 
a Hg viate this inconvenience we have 
eee Tin SG Pk! used a different type of filament, 


whose coils have a diameter just 
greater than the diameter of the 
electron collimating hole. Such a di- 
sposition greatly reduces the sput- 
tering as well as the efficiency, 
though the latter has still a tole- 
rable value. 

Naturally an effective evaluation of the optimum operating conditions of 
the source is related to a study of the beams at the collector of the separator. 


Fig. 7. — Effect of reflecting electrodes on 
total source efficiency: X both electrodes 
at cathode potential; e without reflecting 
electrodes; 0 electrode J at cathode poten- 
tial; electrode i at ground potential. 
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Thanks are due to prof. G. BoLLA for his kind interest in this work. 


RIASSUNTO 


Viene descritta una sorgente di ioni a campo magnetico e catodo caldo per un 
separatore elettromagnetico di isotopi. La sorgente, progettata sia per sostanze solide 
che sostanze gassose, può fornire una corrente totale di ioni di circa 20 mA. Vengono 
riportate misure sulle caratteristiche di scarica, sulle correnti estratte e sul rendimento 
della sorgente e viene discusso l’effetto di elettrodi riflettori, posti nella camera di 
scarica, sul funzionamento della sorgente. 


LETTERE ALLA REDAZIONE 


(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lasciata 
dalla Direzione del periodico ai singoli autori). 


The = —p + y Decay Induced by Heavy Vector Mesons (*) (*). 


P. PRAKASH 


Centro Brasileiro de Pesquisas Fisicas - Rio de Janeiro 


A. H. ZIMERMAN 


Instituto de Fisica Teorica - Sao Paulo 


(ricevuto il 27 Dicembre 1958) 


A Universal Vector minus Axial Vector coupling for the Fermi interaction among 
four fermions was proposed by several authors (1). It was suggested by Feynman 
and Gell-Mann that such Universal interaction may arise from the exchange of 
« electrically charged » vector mesons of very high mass between the two currents 
formed by the two pairs of fermions. LeITE LoPES (?) suggested the possibility 
of including the neutral counterpart of the vector meson in the theory also, provided 
that we assume that the neutral meson enters the coupling with the fermion field 
in such a way that the current which creates the neutral meson is conserved. The 
neutral meson then leads to an effective non-electromagnetic potential between 
neutron (or neutrino) and electron. 

In a theory which allows only direct four fermion interaction the radiative 
decays of the type u ~e+y, X>p-+y¥ etc. involve at least the square of the Fermi 
coupling constant G in the matrix element and consequently are very slow. How- 
ever, if the interaction is due to the existence of vector mesons (?) these decays are 
possible with only the first power of the coupling constant (or rather eG) in the 
matrix element. These decays may not be then very small compared to the cor- 
responding other modes of decays known experimentally and involving the first 
power of the coupling constant G. FEINBERG (+) has shown that, through the vector 


(*) Supported in part by the Conselho Nacional de Pesquisas of Brasil. 

(**) Presentato al XLIV Congresso della Società Italiana di Fisica, Palermo, 6-11 Novem- 
bre 1958. 

() R. P. FEYNMAN and M. GELL-MANN: Phys. Rev., 109, 193 (1958); E. C. G. SUDARSHAN 
and R. E. MARSHAK: Padua-Venice Conference on Mesons and Recently Discovered Particles 
(September 1957); J. J. SAKURAI: Nuovo Cimento, 8, 649 (1958). 

(2) J. Leite Lopes: Nucl. Phys., 8, 234 (1958). 

() The vector meson must be of very high mass in order to be unstable for the decay into 
known particles [reference (?)]. This is also implied by the close agreement of the experiments 
with the predictions of the local Fermi interaction theory. 

(4) G. FEINBERG: Phys. Rev., 110, 1482 1958). 
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meson, the branching ratio of pu e+Y to u>e+ v+v is — 10-4 while that for 
ut-et+et+e- arising from the internal conversion of the photon is ~ 10°. 

In this note we have calculated the branching ratios for the decay 2* =p+y 
which is also possible through the intermediate vector meson. The 2+ decays vir- 
tually into a charged vector meson and a neutral baryon which then recombines 
to create a proton. The photon may be emitted due to the interaction of the electro- 
magnetic field with the charged meson or with the ingoing or the outgoing charged 
particles (5). 

The interaction Hamiltonian for the coupling of fermion fields with the vector 
meson field is given by (°) 


FE (x) = JPaYo(1 un Y5)V5Po an GPVo(1 + Y5)VnPo + h.c.. 


The quadratically divergent terms as well as the non-gauge invariant terms 
disappear when the contributions from the three diagrams are summed and the final 
result is only logarithmically divergent. The matrix element Ql is given by 


Geena a de 
WM =e ya (Ps) = D[(1— y5)mMe + (1+ 5) 4] Op (E phy — yy) + | U=(p>)- 


Here p, and p, are the four momenta of the proton and the = of masses m, and my 
respectively and k is the four momentum of the photon of the polarization four vector e. 
D = (1 fis?) {d4q/(q? + M3)? is the logarithmically divergent integral and F is the 
finite term. When the mass of the vector meson is assumed bigh (say ten proton 
masses) F can be shown to be proportional to the inverse of the square of the vector 
meson mass and can be neglected. 

For the branching ratios we obtain (summing over the various intermediate 
baryons), for a cut-off of the order of the mass of the vector meson: 


Pe Y) theo 
a DDT Jin A (A 10 
PACE pe AT) 


and (7) 


FRE > P+ Y) theo. 
P(X 1+ u++V)theo. 


431035: 


(5) For a similar decay of A° into a neutron and a photon, the photon may be emitted by 
the intermediate charged meson or the baryon only. The decay XZ°— A°++y should be treated 
differently since it conserves strangeness and is known to be very fast. 

() The coupling constant g is related to the coupling constant G involved in the direct 
Fermi interaction written in the form 


¥ 


(CRS ne 
FH (x) = Ave [Par +75) Pp] [vere(1 +75)¥al +h.c. 


by g?2/M2= G/4/2. Here M, is the mass of the vector meson [see reference (1)]. 

È no) The leptonic decay of A° (A°— p+e-+v) has recently been confirmed experimentally. 
Similar leptonic modes of decay are expected for 2 with the branching ratios to the known 
pionic modes of the order 107?, [J. LErrE LOPES: Ann. Acad. Bras. Ci., 29, 521 (1957); see also 
FEYNMAN and GELL-MANN: reference (1)]. i 
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It may be remarked that for the muon decay (u >+e+v+¥) the hypothesis of the 
non-locality in (V—A) Fermi interaction (with G,=G,=G) from the exchange of 
vector meson leads to an effective increase in the Michel parameter o while the 
life time of the decay is decreased. In fact we have 


m 


and 


à 3m, mi 
FOE 10 Ve ee | 
5 M? My 


Hence it is suggested that the vector meson should have a mass higher than about 
ten proton masses (see also reference (?)). 


* OK OK 


The authors are grateful to Professor J. Leire Lopes for the suggestion of 
the problem and the continued aid in its completion. One of us (A.H.Z.) wishes 
to express his gratitude to the Instituto de Fisica Teòrica for the grant which made 
this work possible at the Centro Brasileiro de Pesquisas Fisicas. We are grateful 
to the Centro Brasileiro de Pesquisas Fisicas for the hospitality. 


x» OK A 


Note. - An event which, as suggested by M. GOLDHABER, could be interpreted 
as the decay 2 >p+y was found by E. P. GEORGE, A. J. Herz, J. H. Noon and 
N. SoLNTSEFF: Nuovo Cimento, 3, 94 (1955); see R. E. BEHRENDS: Phy. Rev., 111, 
1691 (1958), footnote (4) 
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For nuclear disintegrations initiated 
by cosmic ray particles of great energy 
(E per nucleon > 10" eV), there exists 
at present no satisfactory method of 
estimating the relative proportion of x- 
mesons, K-mesons, hyperons, and nu- 
cleon-antinucleon pairs produced amongst 
shower particles. Since direct mass meas- 
urements of secondary particles are im- 
possible at such high energies, several 
attempts by indirect methods have been 
made. 

A method due to BRISBOUT et al (1), 
and later used by other workers (?), of 
determining the proportion of neutral 
rather than charged K-mesons, hype- 
rons, and nucleon-antinucleon pairs (here- 
after called x-particles) amongst shower 
particles, involves the following assump- 
tions: 


a) because of their very short life- 
time, r°-mesons cannot contribute to 
nuclear interactions; 


(*) Formerly known as P. KUMAR. 

(1) F. A. BRISBOUT, C. DAHANAYAKE, 
A. ENGLER, Y. FUJIMOTO and D. H. PERKINS: 
Phil. Mag., 1, 605 (1956). 

(2) B. EDWARDS, J. Losty, D. H. PERKINS, 
K. PINKAU and J. REYNOLDS: Phil. Mag., 3, 
237 (1958). 


b) secondary interactions produced 
by neutral particles are due to x°-par- 
ticles or ejected nucleons; 


c) the mean free path of x-par- 
ticles, nucleons, and z-mesons is the 
same; 


d) nucleons and (or) heavy me- 
(K*, 6°, 6°) with isotopic spin, 
I= 4, are produced in a charge-inde- 
pendent way, so that if one assumes that 
it is true also for the hyperons, then un- 
like pions for which n,+=2n.0, NnE=Vp. 


sons 


Using N°, the number of secondary 
interactions produced by neutral par- 
ticles and N=, the corresponding number 
of interactions produced by charged par- 
ticles all observed within the same an- 
gular interval (defined by a strip of width 
720 um and length — 10 em) within the 
narrow cone of high-energy showers ini- 
tiated by singly and doubly charged par- 
ticles of energy = 10!° eV, BRISBOUT 
et al. (1) obtained the ratio N°/N* as 
O = 0.25. Because of assumption (c) 
above, Q= nan or Nyof(nE+7*), 
where n_, ny and n, refer to the numbers 
of x-mesons, x-particles, and the total 
number of charged shower particles, 
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respectively. While considering the pos- 
sibility that the large proportion of ob- 
served non-pions could be due to nucleon 
contamination, they (1) discussed the pro- 
bability that nucleons produced from the 
target nucleus as a result of multiple 
collisions within it, would be found in 
the region of the narrow forward cone, 
but in view of the high multiplicity in- 
volved at such high energies, estimated 
that such a contribution from the nuc- 
leons would be negligible. No account 
seems to have been taken of the nucleons 
of the primary, which as we shall see 
would continue undeviated or only 
slightly affected in a collision even with 
a single nucleon of the target. It is our 
aim to consider presently the processes 
by which high-energy neutrons of the 
projectile enter the forward cone and to 
see whether such considerations affect 
the result. Since charged secondary inter- 
actions are produced by both x* and 
x+ while the neutral interactions are due 
only to the x°-particles (r°-mesons being 
excluded because of assumption a) above) 
the relative proportion of x-particles is 
much more sensitive to the presence of 
a few neutrons than it would be if an 
equal number of protons were present. 

There are two ways by which neu- 
trons can enter the forward cone of high- 
energy interactions. These are: 


i) If the primary is a complex 
nucleus, an «-particle or a heavy singly 
charged particle (deuteron or triton), 
only a part of the constituent nucleons 
may participate in the first interaction 
and the remainder proceed in the shower 
with the same energy per nucleon as 
that of the primary. 


ii) The primary nucleon may par- 
ticipate in the collision and emerge from 
it with an energy still higher than that 
of any of the nucleons from the target. 


We are at present not particularly 
interested in heavy primary interactions 
for which it follows naturally that during 


ne 56 - JI Nuovo Cimento. 
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collisions with nuclei only a fraction of 
the incident nucleons take part in the 
interaction. For «-particle interactions 
within a wide range of primary ener- 
gies (°°), it is known that, during col- 
lisions with emulsion nuclei, on the aver- 
age only half of the incident nucleons 
contribute to meson production. In col- 
lisions with lighter elements or peripheral 
collisions with heavy targets (i.e., events 
of smaller number of heavy prongs), one 
would expect a still smaller fraction to 
take part in the collision. It is reasonable 
to assume that all such nucleons which 
do not take part in the first interaction 
would appear within an angle « to the 
primary direction where « ~ (4)p,/po 
(p; is the internal momentum of the 
nucleons within the nucleus and py is 
the external momentum per nucleon of 
the primary). For p; corresponding to 
the average binding energy per nucleon 
(— 8 MeV) and p, corresponding to the 
energy per nucleon of the incident «-par- 
ticle, E (in eV), the angle « (in radians) 
is given by the relation «~108/H. For 
«-particles of energy from 300 MeV to 
> 10 GeV per nucleon, several workers (°°) 
have experimentally obtained the angular 
distribution of singly charged particles. 
From these data it is possible to find the 
energy dependence of the maximum angle 
within which are found most of the nu- 
cleons of an «-particle which presumably 
do not take part in the first event. We 
find this energy dependence to be in 
agreement with the relation for « derived 
above. Consequently, this relation may 
be used also for the energy region we 
are interested in at present, and one 
concludes that the nucleons of the 
projectile, which do not take part in the 
primary event, would be able to contri- 


(3) M. CECCARELLI, G. QUARENI and G. T. 
ZORN: Nuovo Cimento, 1, 669 (1955). 


(4) G. QUARENI and G. T. ZORN: Nuovo 
Cimento, 1, 1282 (1955). 
CANNES RS APPAMRAO RP DANIEL 


and K. A. NEELAKANTAN: Proc. Ind. Acad. Sci., 
43, 181 (1956). 
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bute to secondary neutral interactions 
in the region under investigation. On 
account of the secondary processes 0c- 
curring in the uppermost layers of the 
atmosphere, it seems reasonable to as- 
sume that some fraction of the high- 
energy interactions initiated by singly 
charged particles in the stratosphere may 
in fact be due to deuterons or tritons (*). 
Similarly as for multiply charged pri- 
mary events, the nucleon which does 
not participate in the primary collision 
is likely to produce secondary inter- 
actions. The calculation of this contri- 
bution is not straightforward and so 
we shall not enter into more detail 
here. 

As to process ii), the nucleons that 
have participated in the primary colli- 
sion can produce secondary interactions 
only if they appear in the angular in- 
terval defined by the scanned volume. 
A number of factors, such as the nature 
of the primary and target nuclei involved 
in the collision, inelasticity of the inter- 
action, and primary energy determine 
this probability. In the first approxi- 
mation, let us consider for example an 
elastic nucleon-nucleon collision. At re- 
lativistic velocities, the angle between 
the two departing nucleons is given by 
2y./(y2 — 1) < tg 8 < co, where y, is re- 
lated to the velocity f, of the c.m. system 
by the relation y,=[1—Bî]-?. (It may be 
seen that the factor 2y,/(y5—1) reduces in 
the non-relativistic case to & 4/W, where 
W is the kinetic energy of the nucleon 
in terms of its rest mass, so that the 
angle 0 is always about 90°, as is well 
known). The minimum angle between 
the two nucleons is ~ 2/y,, so that both 


(*) In such a case the stripped off nucleon 
would proceed without deviation or loss of 
energy, so that if it interacted, the second inter- 
action would have similar characteristics as the 
first event. An event of this type observed in 
an earlier work (°) led us to some interesting 
conslusions, which we shall describe separately 
while considering the probability of the presence 
of deuterons and tritons in the stratosphere. 


ADITYA 


the nucleons will only rarely appear to- 
gether within the narrow core. It also 
means that, if one of the nucleons makes 
an angle between 1/y, and zero, the other 
nucleon will appear correspondingly be- 
tween angles from 1/y, and 90°. In 
actual practice, the collisions are no lon- 
ger elastic, but since at very high ener- 
gies the primary particle does not lose 
more than a small fraction of its energy 
in meson production, the likelihood of 
the primary nucleon appearing in the 
narrow cone of half-angle 1/y, is very 
high. On the same grounds, the struck 
nucleon of the target is not likely to 
appear in the scanned volume unless 
there are multiple collisions. If one 
makes further the reasonable assump- 
tion that in half of the cases the inter- 
acting nucleon would suffer a charge 
exchange, then both, protons undergoing 
charge exchange and neutrons which 
leave the interaction after suffering only 
energy loss, would contribute to secon- 
dary neutral interactions. 

Concluding that the nucleons from 
the projectile (whether or not they have 
participated in the primary collision) can 
contribute to secondary interactions, let 
us look at the data of BRISBOUT et al. (À) 
for which some experimental details are 
available... 9 neutral interactions were 
observed while following 166.2 cm of 
axial length along showers produced by 
singly charged particles and 82.1 cm of 
the shower axis for «-particle stars. In 
the same region, 36 charged interactions 
were also observed. If it is assumed 
a) that one neutron is present for a 
fraction F of the length 166.2 cm, b) that 
two neutrons are present along 82.1 em, 
and c) that the interaction mean free 
path of the neutron is the same as that 
of the proton, i.e., A, = 34 cm, one finds 
that even for F= 0.5 (which corre- 
sponds to the presence of an exchange 
neutron in half of the cases and takes 
no account of the possible contribution 
of neutrons resulting from heavier singly 
charged primaries) 7.3 interactions out 
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of the 9 can be attributed to neutrons, 
thereby reducing the proportion of x-par- 
ticles to about 5% of the shower particles. 
This value may further be modified if 
account is taken of the possibility of 
deuteron and triton primaries. In view 
of the limited statistics such detailed 
considerations may not be meaningful, 
but it appears that the proportion of 
x-particles determined from the ratio of 
neutral to charged secondary interac- 
tions is highly sensitive to the presence 
of nucleons amongst the shower par- 
ticles. 


* OK OK 


Some of the arguments presented here 
are based upon an earlier work (*) done 
at the Panjab University, Chandigarh, 
India. 

Thanks are due Professor B. M. 
ANAND for encouragement and to the 
Government of India, Department of 
Atomic Energy, for financial assistance. 


(5) P. KUMAR: Proc. 45-th Ind. Sci. 
III-3, Abs. n. 9, Jan. 1958. 


Cong., 
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The author is indebted to Professor B. 
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Bombay, and for many stimulating dis- 
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to thank Professor NIeLs BOHR for the 
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Note added in proof. 


Since this paper was written, the 
results of an investigation by LOHRMANN 
and TEUCHER (Phys. Rev., 112, 587 
(1958)) have come to our attention. 
Amongst secondary particles produced 
at an average energy of ~ 101° eV, they 
obtained the ratio of interactions pro- 
duced by neutral to charged particles 
as (16+6)%. Their procedure of cor- 
recting for the presence of nucleons has 
been different from ours, and has yielded 
a value of (9*$)% for the proportion of 
strange particles and anti-particles among 
charged secondary products. This result 
is in agreement with our present consi- 
derations. 
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Zur leitfähigen Haut des Titanoxyds. 


S. YAMAGUCHI 


The Institute of Physical and Chemical Research - Tokyo 


(ricevuto il 2 Febbraio 1959) 


Auf eine Glasfliiche wurde Titan durch der Haut befindlichen TiO,-Kristalle auf 
Vakuumverdampfung niedergeschlagen. dem Substrat orientieren. 
Der erhaltene Titanfilm lie sich mit der 
ozonhaltigen Luft auf Zimmertempe- 
ratur oxydieren. Die auf diese Weise 
hergestellte Haut (Dicke: ca. 300 AE) 
lieferte ein in Fig. 1 wiedergegebenes 
Beugungsbild. Sowohl die Beugungs- 
flecken als auch die Ringe in diesem Bild 
entsprechen dem Titan-Oxyd (Ti0,) des 
Rutiltyps. Die den Flecken entspre- 
chenden Kristalle sind regelmäfig orien- 
tiert. indem die (110)-Ebenen derselben 
der Substratfläche parallel hegen. Dieser 
Oxydfilm war elektrisch leitfàhig. In der 
Tat erkannte man keine Elektronen- 
aufladung beim Beugungsversuch. Die 
betreffende Oxydhaut wies eine bemer- 
kenswerte Adhäsion für das Substrat 


Fig. 1. — Reflektionsbild der leitfähigen TiO,- 
È i n x pers Haut. Wellenlänge: 0.0323 AE. Kameralänge: 
auf. Dies spricht dafür, dafi sich die in 495 mm. Positiv in 2.3 facher VergròBerung. 


n 


IL NUOVO CIMENTO MOLAEXTIPENEEC 16 Marzo 1959 


A Remark on One-Dimensional Quantum Electrodynamics. 


V. GLASER and B. JAKSIG (*) 


Theoretical Study Division CERN - Geneva 


(ricevuto l’11 Febbraio 1959) 


In a recent note, I. BrALyNICKI-BIRULA (!) has shown that the one-dimensional 
quantum electrodynamics with zero mass electrons can be solved exactly. In what 
follows we would like to show that this one-dimensional model is essentially equi- 
valent to the THIRRING model (?). 

The field-equations in this model are as follows: 


(la) PY ey", A, ; 
(1) (= p)An = —j = —epyny, 
where p = — iy"Ou3, y = io,, y? = 02 (7), and y is the mass of the « photons ». 


Tt will turn out that the case «x = 0 has no physical solution due to the fact that 
the Coulomb force in the one-dimensional case does not vanish at infinity. (1 a, db) 
should be supplemented with the usual (non-renormalized) canonical commutation 
relations and the subsidiary condition %“(«)|a> = 6, y = 0,4%, |a> being a physi- 
cally allowed state. 

From (la) there follow the conservation laws: 


CARE 0,0, 


where 


1 
sh = 25 = py'y = (01— 22) A+ Oz) > 


du = py" yy = (01+ 02: O1 — 02) > PESO 


(*) On leave of absence from the Institute « Ruder Boëkovié », Zagreb, and the University 
of Zagreb. 

(3) We are indebted to the author for sending us a preprint of his letter to the Nuovo Cimento. 

(2) W. E. THIRRING: Ann, Phys., 9, 91 (1958); Nuovo Cimento, 9, 1007 (1958); V. GLASER: 
Nuovo Cimento, 9, 990 (1958). 
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With the introduction of new variables w= t+%, Vv = tx, the conservation 
laws lead to 


f @ = alr) = vile —t, O pile — à 0), 
9 
se | = ou) = ya(x +t, 0) per + t, 0). 


0: 


This means Lls¢ = Lite = 0, and 
(3) [j“ (x), j’(x')] = 0 for all x and 2. 
The last fact follows from (2) and the canonical commutation relations. 
From this basic result we see that (1/u?)j“(x) is a solution of (16). The most 


general solution of (1b) will be 


1 
(4) AM) = A+ A 


where /4“() satisfies the homogeneous equation (L]— y?)4%(x) = 0. From (3) and 
the canonical commutation relations it follows furthermore 


(5) AG) (Dl CAS), 0 
and 
(6) LA# (©), A°’(x')] = [Ab(a), 46(0)] = tg” Ala — x!, 3) . 


This all implies 44 — 45 = 45w, so that the model does not give rise to any 
creation or scattering of « photons ». 
Inserting (4) into (la) we obtain 


e? 
(7) PY = Aou?" Y + ie (Pye ya Y - 
We introduce now the « Ansatz » 
(8) pla ) = exp |i À AG Of ale 
where 
A(æ) = OF AS 3 ve” Oy Ao, Ré ipy Ao, ? 
sn e, Mo +1, . [A(m), Aw)]=0, 


for all + and a’ by virtue of (6). One sees immediately that py“y = yryl{yr and 
the equation (7) leads to 


e? e2 


(9) PYr = PALE, Yr)yV" Wo = 2 o (WrWr)Wr - 


This is exactly the Thirring equation for y, with the coupling constant gr = — e?/u?, 
compare (?) . 
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The «gauge transformation) (8) leading to the Thirring model means that the 
S-matrix, defined by ww = S1wnS, Aout = S14,,8 = Ain, is identical with the 


S-matrix of the Thirring model: S = exp {igrQ,0,}, compare (2). It is a unit matrix 
with respect to the « photons ». 


In the case 7 = 0 equation (1b) can be written in the form (JA, = — eg, with 
A, = 4(A? + A1). The formal solution 
eu 
AE Sp [eset ae’ A (uw, Vp) > 


contains a factor which increases linearly with time, so that the asymptotic condi- 
tions cannot be satisfied. The physical reason for this behaviour has already been 
mentioned before. 

The study of the Green’s functions and the Wightman functions in this model 
does not seem to reveal anything essentially new as compared to the Thirring model. 


In conclusion, we would like to thank Dr. W. THIRRING and Dr. F. L. SCARF 
for discussions, which led to the clarification of the physical content of this model. 
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The Attractive K--Meson Proton Interaction (*). 


E. M. FERREIRA (**) 


Department of Mathematics, Imperial College - London 


(ricevuto il 12 Febbraio 1959) 


Experiments (1?) seem to indicate 
the K~-proton interaction to be attrac- 
tive. However it is known that the lowest 
(2nd) order perturbation theory matrix 
element K p+Y—K p, where Y represent 
A or X (virtual) hyperons, has a sign that 
corresponds to a repulsive potential be- 
tween K° and proton for either parity 
of A and © hyperons. The question 
then arises whether or not the higher 
order matrix elements give contributions 
of the correct sign and large enough to 
make agreement with experiment pos- 
sible. 

We have calculated one of these 
higher order diagrams, that obtained 
from the lowest order one by adding the 
exchange of a m-meson between the inci- 
dent and the outgoing proton. We spe- 
cialized to the case of forward and thre- 
shold (zero kinetic energy) scattering. 
The ratios M®/M® of this fourth order 
matrix element to the second order one 


(*) Supported by the National Research 
Council of Brazil. 

(**) On leave of absence from the Centro 
Brasileiro de Pesquisas Fisicas, Rio de Janeiro. 

(1) W. ALLES, N. N. Biswas, M. CEcca- 
RELLI and J. CRUSSARD: Nuovo Cimento, 6 
DU (MOTO 

(2) R. H. DALITZ: Proc. of the Geneva Con- 
ference (1958), p. 187. 


> 


for the several possible parities of the 
A and & hyperons are indicated in the 
Table below. 


Trsn|Trgen MO/MP 
1 1 Gi Gian <= 5GixN 
4nhc Gian + Gion 
aca ee Gian + 4.24Gieon 
° | 4rhe GNT 0.125625. 
ed Fe q.41 Fran + 5-06kex 
Ms 4ithe i Ca + 8 cx 
oe hi Gr 0.4) ras + 56kon 
ù i 4ahe Gian a Grox 


We note that the ratio is always ne- 
gative for any parities and values of the 
coupling constants and that for G2/4aie= 
— 15 its extreme values are — 0.9 (for 
Gon = 0, Tgan=1) and — 30 (for Gran=0 
Tken=%). 

The matrix element for the graph 
considered presents thus the desired 
features, having the convenient sign and 
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being large enough to compensate the 
second order term. 

The K -proton system is a mixture 
of T—1 and T— 0 isotopic spin states. 
It is interesting to note that if we assume 
that only one of these states contribute 
to the interaction, the potential derived 
from our graph is still attractive for all 
possible parities and values of coupling 
constants. That is, no minus sign arises 
in the splitting into iso-spin states, For 
example, in the case of equal parities 
attributed to the A and & hyperons, the 


Le] 
Le] 


T=1 state will give a matrix element 
proportional to Gz,y+2Gfryx, and the 
T= 0 one will be proportional to 3Gisy: 
This equality of signs of the 7=1 and 
T= 0 interactions is in agreement with 
the result of the phenomenological ana- 
lysis presented by R. H. DALITZ (see 
reference (?), p. 189). 


*o** 
We are greatly indebted to Professor 


A. SALAM and Dr. P. T. MarrHEews for 
guidance and encouragement. 
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Quadratic Lagrangians in Relativistic Hydrodynamics. 


F. HALBWACHS, P. HILLION and J.-P. ViGIER 


Institut Henri Poincaré - Paris 


(ricevuto il 12 Febbraio 1959) 


In a preceding paper (1) we proposed a lagrangian formalism for relativistic fluids 
made of extended rotating dropplets, defined within time-like tubes along the lines 
of flow. Because the canonical tensors deduced from our lagrangians and satisfying 
the well known conservation relations O,t,,=0 and fu = tt ty,) take the 
particular forms 


py 
= == er 
(1) Luv si July and uv ane 2 Mur > 


the vector g, (linear momentum density) and m,, (angular momentum density) 
obey the equations: 


(2) In = 0 and Mur = Gud JyUn » 


as a necessary consequence. If we suppose the internal tensions between the indi- 
vidual droplets to be negligible, the latter satisfy consequently Weyssenhoft’s 
equations of motion: 

G,=0, M y= Gu, — Gu 


wy DO 
However the proposed lagrangians with a proper rotation energy term 3.M abi ba 
rested on the special assumption that M u did not depend on the derivation of 
the EX variables be, so that the non relativistic approximation furnished as rotation 
energy a term depending linearly of the angular velocity 0955 =b Pb — HOH), 
thus different from the usual classical expression 


(3) T= 310,0, 
where @;= €;;,@;, and I,, denotes the usual Symmetrical tensor of inertia. The 


quantization of such an expression would naturally raise special difficulties. 


() F. HALBWACHS, P. HILLION and J.-P. VIGIER: Nuovo Cimento, 10, 817 (1958). 
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We shall now show that the preceding assumption is not necessary and that 

more general Lagrangians quadratic with respect to ®,, give correctly, not only 

relations (2), but also the usual non relativistic rotation energy (3). We first define 
the relativistic generalization of the tensor of inertia through the relation 


Myg= T8 760,5 (AE 


where M,, is the particle angular momentum and the tensor Z,4.,5 is antisymmetric 
ha Re 


with respect to « and f, and to y and è, and symmetrie with respect to substitution 
of the pair of indices (#5) with (yò). We take as rotation energy density 


= j — GAC Xia) FAC AGS) 
T= $1485 @ap%y5 = 4MapVag » Ogg = 3(ba di — bg da ) - 


And then the non relativistie approximation gives evidently the above classical 
quadratic expression. 

On the other hand, we make the supplementary assumption that this tensor 
I,pys does not depend on the derivatives 9,0 of the EX variables. 

Let us start from the lagrangian 


ZO SALO (Op) 
L= OMC? LE icob, 0,5 = +elxByò WaBAyd ar PAU by ie Suv) x 


Because 2=0 (as a consequence of 22/29 —0 and LOK) —5,,) the two canonical 
tensors ¢,, and f,,,, depend only on the derivatives a £ with respect to the deri- 
vatives } (and also 3,8) namely: 


E É . 3 5 . 
Îuva = Bo) 3 oie bY Mb 08 È 3 (au By DE Ò pn Oxy) ’ 


that is, because by be = dg 
fav = 30M, ZE 
tuy = i00bV 0,8 + Hoicb YP, M,, 0,00 = 


ay 


= ou,(0,S + $M, 20,0509) . 


Thus we obtain the correct form for these two tensors and our formalism 
satisfies the equation of Weyssenhoff with 


Iu = 0(0,8 + 4,00 bY) , My = OM y: 


@) The symbols + and — dencte skew symmetric and symmetric tensors. 
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Note on Possible Parity Non-Conservation in Strong Interactions. 


L. M. Brown () 


Istituto di Scienze Fisiche del Universita - Milano 
Istituto Nazionale di Fisica Nucleare - Milano 


(ricevuto il 23 Febbraio 1959) 


As experimental tests of the conser- 
vation of parity in strong interactions 
have not been completed, it may be of 
interest to observe that a principle which 
has already been successfully applied to 
the theory of the weak interactions may 
allow a distinction to be made between 
the strong interaction terms of the 
Lagrangian which do, and those which 
do not conserve parity. We refer to 
the principle of chirality invariance (1), 
or its equivalent expression in terms of 
the two-component fermion theory (2). 
Whether the presence of such parity 
non-conserving terms in the strong inter- 
action Lagrangian would necessarily 
imply prominent parity non-conserving 
effects in pion-nucleon or nucleon-nucleon 
interaction at moderate energies is an 


(*) On leave from Northwestern University, 
U.S.A., for the academic year 1958-59 as a 
Fulbright Research Scholar in Italy. 

(*) J. TIOMNO: Nuovo Cimento, 1, 226 (1955); 
S. WATANABE: Progr. Theor. Phys., 15, 81 (1956); 
E. G. C. SUDARSHAN and R. E. MARSHAK: Phys. 
Rev., 109, 1860 (1958); J. J. SAKURAI: Nuovo 
Cimento, 7, 649 (1958). 

@) R. P. FEYNMAN and M. GELL-MANN: 
Phys. Rev., 109, 193 (1958); L. M. Brown: 
Phys. Rev., 111, 957 (1958). 


open question related to the role which 
virtual strange particles play, in general, 
in such processes (8). 

In order to apply the chirality inva- 
riance principle, we may suppose that 
the strong interactions, like the weak, 
can in some sense be represented as 
direct four-fermion interactions and adopt 
the model of the pion as a nucleon- 
antinucleon compound, following FERMI 
and YANG (4). Similarly we represent 
the K-meson as a compound of nucleon 
and lambda (5). 

A consideration of this type on the 
part of OKUN (°), but assuming parity 
conservation in the strong interactions, 
led him to propose that not all the 
known baryon fields are fundamental, 
but that the Y and the E are compounds 


(5) A. PAIS: Phys. Rev. Lett., 1, 418 (1958). 
Pais has discussed how parity nonconserving 
effects can be limited by imposing other exact 
or approximate symmetry principles. 

(*) E. FERMI and C. N. YANG: Phys. Rev., 
76, 1739 (1949). 

(5) S. SAKATA: 
686 (1956). 

(°) L. B. OKUN’: Proc. Annual International 
Conjerence on High Energy Physics at CERN 
(L953) ED 2225 


Progr. Theor. Phys., 16, 
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of the others. Here, instead, we assume 
that each baryon field is fundamental 
and see what the postulate of chirality 
invariance tells us about parity con- 
servation. The isotopic spin is omitted 
in the writing as it is irrelevant to this 
question (7). We treat each isotopic 
multiplet given by Gell-Mann-Nishijima 
scheme as a single field. 

The chirality transformation y —7;#, 
y —>— py; (with y =1) is to be applied 
to each fundamental baryon field sepa- 
rately. When the same fermion field 
occurs twice in a bracket, for example, 


(1a) (WO py) (PPO), for 0;=V, A; 


(1b) += (PPO yp), 


TODO == Snel ele, 
independent of the transformation cha- 
racter of the field y under y;. But 
when two different fields occur, 

(2) GO y) + (PMO VV) = 
= + (FOO), 


the equality holding only when y® or y® 
are chosen to be eigenfunctions of y;. 
In order to form a chirality invariant 
four-fermion interaction with other fields 
from expressions like (2), it is clear that 
yD and y® must belong to the same 
eigenvalue of y;. (Note that y;y= +y 
implies — yy;=+y). For (2) to be non- 
vanishing it follows then that O, is 
either V or A. 

Thus while chirality invariance im- 
poses severe restrictions on the inter- 
action form when two different fields 
are in the same bracket, it gives much 
weaker results when the fields are the 
same. We shall minimize the restriction 


(7) Of course, if one requires charge inde- 
pendence, then the isotopic spin is relevant. 
In fact SOLOVIEV has shown that in renormali- 
zable pseudoscalar meson theory, imposing 
CP invariance and charge independence requires 
parity conservation. See V. G. SOLOVIEV: 
Nucl. Phys., 6, 618 (1958). 


10 
el 
Lee) 


(0 0) 
Qt 


in writing the interactions by putting, 
whenever possible, the same fields in one 
bracket. 

Adopting the model of Fermi-Yang, 
we see that the pseudoscalar pion field 
can be represented by an expression of 
type (1b). The « vertex » interaction 


(3) (NysN)-(Ny;N) ; 


where the second bracket represents the 
pion, is thus automatically parity and 
chirality invariant. Pseudoscalar inter- 
actions between this pion and the other 
baryons (X=), (EEr) cannot be written 
in chirality invariant form, but only 
as V or A direct interactions, according 
to Eq. (1). With this restriction they 
can still conserve parity. When three 
fundamental fields are involved, however, 
this is not the case. 

Consider, for example, 1+ N > A-+N, 


(4) (AO;2)-(NO,N) . 


Here the nucleon bracket is chirality 
invariant if 0; is V or A, but it is nec- 
essary for X and A to be positive eigen- 
functions of the chirality operator. In 
this case we get maximum parity non- 
conservation. Similarly, considering the 
K-meson as a nucleon-A compound we 
conclude that (KNA) can be parity con- 
serving, while (KN) is parity violating. 

In summary, we can classify the 
strong interactions as follows, letting 
a, b, c, d represent the four fundamental 
baryon fields, and letting a*, a- be the 
positive and negative chirality parts 
of a, ete. 


Class I: requiring V, A interaction 
and parity violation. For (a, b, c, d) 
we must have a,, b,, ci, d.; for (aacd) 
we must have ¢,, d,. This class includes 
the interactions (NN)(AX) or TAX or 
KN; (NE)(AA) or KEA; (NE)(AX) or 
KEX; (NE)(2Y). 


Class II: requiring V, A interaction 
but not requiring parity violation, form 


886 Tey Mo 


(a a b b). 
or KNA; (NN)(2Z) or x22; (NN) 


OT THe. 


Class IIT: no restriction, form (aaaa). 
This class includes (NN)(NN) or zNN. 


With regard to the question of 
renormalizability it may be remarked 


BROWN 


that, in principle, this question has 
nothing to do with the strength of 
coupling. The objection on ground of 
mathematical consistency is the same 
here as for the usual formulation of the 
weak interactions. If this is resolved, 
for example by the attribution of an 
intrinsic non-locality of some type, the 
same solution could apply here. 
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On Ionization Measurements in a Hydrogen Bubble Chamber. 


I. DERADO and N. SCHMITZ 


Max-Planck-Institut für Physik und Astrophysik - München 


(ricevuto il 9 


Using bubble chamber photographs 
from the 10 in. hydrogen chamber of 
the Alvarez Group (Berkeley) we have 
examined the inelastic = -p-scatterings 
resulting from an exposure of the chamber 
to à m--beam of 1.14 GeV/e from the 
bevatron. Mainly, in addition to mul- 
tiple production, two reactions occur 


(a) cain a e io pil a A 
(b) Ta4-) Tarn ET. 


For most of the events we were able by 
using the conservation laws for energy 
and momentum to decide whether a 
reaction of type (a) or type (b) had 
occurred. However, in some cases it 
was impossible by this method to dis- 
tinguish between the two types of reac- 
tions. (The reasons for this will be 
discussed in a later paper.) 

In order to decide also in these cases 
whether the positive particle was a 
proton or a z+ we have examined the 
momentum dependence of the bubble 
density o both for proton and pion 
tracks, using elastic events 
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where the particles can be identified 
and where the momenta follow exactly 
from the scattering angles. In all cases 
the measured values for o lay between 
20 and 25 cm (Gayerage= (22.5 +2.5) cm) 
regardless of the momentum and of 
whether the track was that of a proton 
or of a pion. This means that no mo- 
mentum dependence for o could be 
established. The reason for this apparent 
independence lies probably in the fact 
that the bubbles had rather large radii 
in our films and that even at momenta 
for which theoretically minimum ioniza- 
tion was to be expected the bubble 
density was so high that the bubbles 
almost bordered on one another thus 
preventing a resolution of higher den- 
sity. Besides, in the case of small 
momenta bubbles had often 
grown together to form one lengthy 
blob or a larger bubble thus making 
the measurement of the bubble density 
unsuitable. 

In the course of these examinations 
we noticed, however, that at momenta 
lower than 650 MeV/c a proton track 
presents a visual form quite different 
from that of a pion track of the same 
momentum. In the case of the pion 


several 


16 Marzo 1959 


888 I. DERADO and N. SCHMITZ 


track the bubbles are comparatively together and are not all lying on a line; 
small and seem to be arranged exactly instead, many bubbles lie outside the 
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Fig. 1. 


on a line whereas in the case of the proper track (8-rays) thus causing a 
proton track the bubbles are larger as quite different visual impression of the 
a result of several bubbles growing track. This difference ceases to exist at 


+ 


5 E + Proton 
(©) 
= o Pion 
À 
k te 
+ 
10+ 
TANT & 
+ 
05 
— I 
0 0.25 


ima A 


888 


ON IONIZATION MEASUREMENTS IN 


higher momenta where approximately 
minimum ionization for both particles 
is to be expected. In a quantitative 
examination we have measured the 5-ray 
density o for many tracks of elastic events 
having a dip angle <50°. A « S-ray » was 
counted if one or more bubbles lay 
clearly outside the track. The results 
of these measurements for 29 tracks are 
shown in the two figures. Fig. 1 gives 
the « d-ray» density o as a function of 
the momentum p, Fig. 2 as a function 
of the velocity £. (Crosses correspond 
to identified protons, circles to pions. 
The curves are drawn as a best fit 
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through the experimental points.) From 
Fig. 1 it can be seen clearly that for 
momenta smaller than 650 MeV/c the 
« §-ray» density of a proton track is 
considerably larger than for a pion 
track so that at these momenta a dis- 
tinction between a proton and a pion 
is possible by the measurement of the 
« §-ray » density. 


* % * 


We are grateful to Dr. K. GOTTSTEIN 
for a valuable discussion and to Dr. L. 
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LIBRI RICEVUTI E RECENSIONI 


G. TORALDO DI FRANCIA — La dif- 
frazione della luce; pagg. XI-+-530; 
figg. 159. Edizioni Scientifiche Ei- 
naudi, L. 8000. 


Il presente testo costituisce quanto di 
più moderno sia disponibile nella biblio- 
grafia internazionale in merito ai feno- 
meni di diffrazione della luce. 

Lo studio teorico di questi fenomeni è 
ricondotto alla ricerca di autofunzioni e di 
autovalori metodo già 
usato in vari articoli originali ma mai 
fino ad ora nei libri di ottica, dove si era 
preferita la più intuitiva via tradizionale. 

Questa nuova trattazione è il pregio 
fondamentale del libro sia perchè rende 
possibile una notevole semplificazione 
della trattazione matematica che risulta 
così più coerente e compatta, sia perchè 
rende più accessibile al fisico moderno, 
abituato dalla meccanica quantistica a 
ragionare in termini di autovalori e auto- 
funzioni, una materia molto specializzata. 

Il libro è diviso in cinque parti delle 
quali la prima è solo una premessa ma- 
tematica per rendere intellegibile il se- 
guito anche ad uno studente con la pre- 
parazione del biennio di ingegneria. 

La seconda, che è forse la parte più 
bella, è dedicata ad un vasto ed accurato 
studio dei fenomeni di Fraunhofer. Vi 
si trova una precisa trattazione di tutti 
i vari tipi di reticoli, dal primo reticolo 
a fili dello stesso Fraunhofer ai reticoli 
di fase; dai reticoli finiti a quelli infini- 
tamente estesi, dai reticoli a trasparenza 


utilizzando un 


a quelli a riflessione, con la giusta pre- 
cisazione che i reticoli incisi sono reti- 
coli a trasparenza e di fase. Non manca 
un paragrafo sulla verifica sperimentale 
dell’esistenza delle onde evanescenti, 
verifica eseguita dall’autore sia con la 
luce sia con le microonde. Completa 
queste parti la teoria dei diaframmi non 
periodici e del potere risolutivo di molti 
strumenti ottici. 

La teoria diffrazionale delle immagini 
sia in luce coerente sia in luce incoerente 
è trattata nella III parte. Si trova qui 
una interessante e approfondita discus- 
sione sul potere risolutivo del microscopio 
e sulle informazioni contenute in una 
immagine in luce coerente. Non ci sembra 
invece sufficiente il rilievo dato alla micro- 
scopia in contrasto di fase. 

Un brillante capitolo sui paradossi 
ottici ai quali può dar luogo la teoria 
della diffrazione è inserito nella quarta 
parte dedicata ai fenomeni di Fresnel. 
Son qui descritti tra l’altro i paradossi 
ottici di Poisson, di Lord Rayleigh, di 
Gouy e la propagazione anomala delle 
onde sferiche con una ingegnosa verifica 
sperimentale. 

Lo studio diffrazionale delle aber- 
razioni si trova infine nell’ultima parte 
che è molto utile per le applicazioni 
pratiche. 

Notiamo uno strano errore tipografico, 
più volte ripetuto che fa diventare infi- 
nitesimali i reticoli infinitamente estesi. 


E. PoLacco 
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